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Branching ratios for the K~ capture by the nucleon are studied to investigate the validity of 
Gell-Mann’s global symmetry. The following assumptions are made: 1) Gell-Mann’s global symmetry 
is valid, 2) K-meson is pseudoscalar, 3) mass differences among baryons can be neglected and 4) 
processes involving K-mesons are described by treating the M. S coupling as perturbations in the 
lowest order that gives the effect. Under these assumptions, the matrix element for the above men- 
tioned processes are expressed in terms of the s—WN scattering matrix elements for which experimental 
data can be used. The branching ratios, then, are compared with experimental results. 


§ 1. Introduction 


Recently Gell-Mann proposed an interesting theory concerning the connection between 
symmetries and strengths of the interactions.” The essential points of his arguments are 
as follows. 

1. All the baryons are symmetrically coupled to the pion field with a coupling 
constant ~15—V. S coupling. 

2. The K-meson couplings to the baryons are weaker and less symmetric although 
still charge independent—M. S coupling. 

We examined how his theory predicts experimental results involving K-meson, pion 
and baryons under the following assumptions. 

1. Gell-Mann’s global symmetry is valid. 

2. K-meson is pseudoscalar. 

3. Mass differences between the nucleon, A’ and »' can be neglected. 

4. Processes involving K-mesons are described by treating the M. S coupling as 
perturbations in the lowest order that gives the effect. 

Under these assumptions, the matrix elements for various reactions involving K- 
meson, pion and baryons can be approximately expressed in terms of the 7—WN scattering 
amplitudes. But unfortunately, besides the —WN scattering matrix elements, there appears 
a matrix element for an unphysical process, which corresponds to the absorption of a 
fictitious pion with [=0. However the phase of this matrix element is expressed in 
terms of the z—WN scattering phase shift if we assume the invariance under the Wigner 
time reversal. We take the magnitude of this matrix element as an adjustable parameter. 

In the present paper, we shall apply this technique to the case of the K~ capture 


by the nucleon. The branching ratios may be compared with experimental results by 
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making use of the experimental data for the =—N scattering. In § 2 expressions are 
derived to be compared with experiment, in § 3 comparison is made with experiment and 
finally in § 4 we discuss the result thus obtained. 


§2. Evaluation of the matrix elements for the processes 
K-+N>Y+42z 


It is experimentally known that there are the following modes of reactions of the 
K--meson absorption by the nucleon at rest.” 


K-+p> 24+" 
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First we evaluate the matrix elements of these modes by treating the M. S coupling 
as perturbations in the lowest order that gives the effect. 


process K~ + p— A’+7°. 


As an example, we take a 


The K-meson interaction Hamiltonian for this case, according to the assumption (2), 


is given by 
Hy =ig ar D7 ypR-+iG a, °7s pK Pigex/ 2 D-7,aK- (2-3) 
which may be rewritten as follows. 
Hy =i(Yax+Gux) /V 2 275 pK- +i(Gax—Jox) /V 2 Y°yspR- +igeav 2 S774nK- 


(2-4) 
where 


M=(Z°4Y)A/2 and M=(Z°—Y)//7. 


Notations used here are the same as those in Gell-Mann’s paper.” 


The above mentioned process is given by the matrix element 


(A x°|S|  py=(Pa |S) 0" foo) dtedea PLC), mes, H(t) ]\K-p) 


n=O 
ty 4 fo) \n ps 
= (—i) (ix! | SR yine 5) Ol... dt,:+-dt, 
JV 20, n=0 n! 


x PLH (1), H(4,)0@) |p) (2-5) 


Global Symmetry and the Branching Ratios for the K-~ Capture by the Nucleon 119 
where 
H() = | (He) +H (2) Jd» 
and 


O=i(9ax+Ixzx) /[V227sp+ i(Jax—Jux) /V 2 Y° eB 12 Gux27Ton- 


This expression can be transformed to the Heisenberg representation as” 


(Mx 


py=(—) (Aa Hf ee Zeon ““O\p) (2-6) 


where | 4°) or |p) represents the physical particle with complete K-meson and pion clouds 
around it, O(x) is the operator in the Heisenberg representation and k, is the momentum 
of the incident K~-meson. The above expression is exact. We next rewrite (2:6) as 
follows, using the assumptions (1), (3) and (4), 


(A'2"|$| K-p)=—i//2 | re 


e* | (Z°"|O(x)|pyt+¢Y° z°|O (x) | p>] 


ax 


=—ih72 | e*| ¢ pz*|O! (x) | p) + <n 2"|O" (x) | p>] (2-7) 


where O’(x) corresponds to the operator O(x) in which Z°, Y° and 2 are replaced 
by p, n and n respectively. Or, explicitly we may write 


O' (x) =1(YaxtIsn) /V2N ¥5(1+t,)/2N+i(Gar— Iux)/V 2 Ny; c_N 
FIV 2 Geax NG —) /2N (2-8) 
where the operator NV represents the nucleon wave operator in the Heisenberg representation. 


. < . Y rae 
Now, let us examine one of these terms in detail, say (—)<pa"||d*x/v 2wx 


e** Ny,t,N|p). As was shown by Low,” this quantity, on the energy shell, is es- 
sentially the —WN scattering matrix element (7°p|S|z°p), that is to say, 


(—)(pa| | 7c N7,2,N|p)=— Caps =p) (2-9) 


where k,?+/2=0, #4: pion mass and y,, is the coupling constant of the —WN interac- 
tion. In our case, however, k, is the momentum of the incident K--meson, so the 
matrix element (2-9) describes the off-the-energy-shell process. Nevertheless we can make 
the connection between this and the actual process, since baryons are heavy compared 
with mesons and we may neglect the spatial momentum dependence as far as the s-wave 
is concerned. Then eq. (2:9) is valid approximately even for this case. As for the 
treatment of the p-wave, we shall state briefly in § 3. 


Another term concerned here is 


Ea te Sa 
(—)(pa" | $e Wy, 1p). 
JV 20, 
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This corresponds to the absorption of a fictitious pion with I=0. (We shall call 


this fictitious particle as o-particle hereafter.) So we may write this as 
d‘x dec AT ) 1 

— Po hee OY BEN = 
( dP lean Sask Ge 


where YJ, is the o-particle coupling constant, which we assume to be so small that g,° 
can be neglected. As this matrix element does not correspond to the experimentally 
measurable quantity, we take this term as an adjustable parameter. In general this term 
will be complex, so there are two independent parameters. However, as we shall see 
below, the phase of this matrix element is expressed in terms of the z—WN scattering 
phase shift if we assume the invariance under the Wigner time reversal. Thus, two 
independent parameters could be reduced to one. This follows from the fact that the 


reaction matrix K is real and symmetric if K matrix is invariant under the Wigner time 


n° |S| po) (2-10) 


reversal.” 
For the s-wave, K matrix has two channels for [=1/2 submatrix and is expressed 


as 

-oN 

K=(, 
tan 0, hoaziv 
where @ is real. Then the T-matrix is given by 
0 ce cos 0,-a 

tore an 

‘cos0,-@  e'®* sind, 


Thus, we see the phase of the scattering matrix for the processes ¢+p—> p+ is 
given by 0,. 
Finally, using eqs. (2-9) and (2-10), eq. (2-6) may be expressed as 


(AP? |S| Ko py =G,/4 {¢ pa 


WO) +V 2 <nz|S| px-)} +G,/4 ¢ px |S| por) 


(2-12) 
where 


G,= (Gax—Iszn) / In 
Gy= (Jax+39zx) /In- 


Similarly, the matrix elements for other processes are expressed in terms of the 7—N 
scattering amplitudes as follows : 


(208° |S| Kp) =G,/4 {C pa? |S] pa) — 2 (nz |S| pz-)} +G,/4 ( px |S| po), 

(at 2" |5| Kp) =G,/2< pam |S| pro), (2-13) 

(2° a" |S] Kop) =G,/22 (nm |S| px) + G,/2y2 (nz* |S por), 

(Mar |8] Kn) =G,/4 {C pa- |S] pz-) +2 (nz- |S| nz-y} +G,/4 ( pr |S| no), 

(202 |S] Kn) =G,/4 {¢ pa |S| nz) 2 (ne o/4 (px |S| no, 
(2-14) 
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(37 29|S| Kn) =G,/2y/2 (nz"|S| nt) + G_/20/2 (pre |S| no). 


Introducing the eigenamplitudes of the isospin for the 7—WN scattering, we may 
write these as 


(Ar? |S| K- p) =G,/12(4M,—M,) +// 3 G,/12 M,, 
(3°n° |S| K- p)=G,/4 M,+/3 G,/12 M,, 
(3+ 27 |$| K- p)=G,/6(2M,+M,), (2-15) 
»'-x* |S| K- p)=G,/6(M,—M,) +3 G,/6 M,, 
{A x- |S| K-n)=G,/61/ 2 (4M,—M,) +G,/2/6 M.,, 
(3°x-|$| K-n)=—G,/6/ 2 (M,+2M,) —G,/2V/6 M,, (2-16) 
(¥- x°|S| K-n)=G,/6)/ 2 (M,+2M,) +G,/2/ 6 M, 
where 
M,=(I=3|S| I=) 
M,= (I= $|5|I=3) 
M,=(1=3 (N, 7) |5| 1=3 (N, ¢) >. 


Note that between these matrix elements there are the following relations owing to 


the charge independence. 
3” x* |S| Kp) + (2 a [S| K pp =2 (2°? |S] Kp), 
(23°n- |S|-Kon)+ CS~ 2% |S| K> ny =0, (2-17) 
V7 2607 \S| K- py = M2 |S| K-n). 


§ 3. Comparison with experiment 


Experiments on the processes (2:1) and (2-2) have been done by means of the 
hydrogen bubble chamber and the emulsion, respectively.” 

As the processes (2-2) are studied using bound neutrons, their measurement is 
more difficult and uncertain. So, we first examine the branching ratios for the processes 
(2-1) :i.e., K-+p—>Y+7. In these processes, K~-mesons are absorbed at rest, and we 
may expect that the K~-mesons are absorbed from the s-state. 


In this case, we may write eq. (2-15) as 


M, =-1-e** sind, ; M. matt gtts sin 0s ; 
1 13 Sait 
M,=——e cos0,:@ (3-1) 


; pane Ml eats, 
where 0, and 0, are the s-wave 7—N scattering phase shifts in the states of isospin } 


and 3, respectively. 
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As we have neglected the mass difference between A® and se, it - not certain which 
energy of the pion we should take to evaluate the value of 0, and os: For a 
the outgoing pion energy is about 80 Mev for the processes K-4+N 3247 and i about 
150 Mev for the processes K-+N>A+7. If the incident K~-meson energy 1s trans- 
ferred to the corresponding pion energy, the incident pion kinetic energy would be about 
350 Mev. So, we varied the energy from 100 Mev to 400 Mev, and at each —— 
examined the branching ratios. To evaluate the value of 6, and 0;, we used Orear’s 


experimental formula” 


0,= (0.167 + 0.012) 7, 
0,= (—0.105 0.010) 7 (3-2) 


which seem to hold up to fairly high energy.* 
Then, from eq. (2-15), we have 


[( Ax? |S| K- p)|?=C[Int+ 6+, C7], 
| <2 2° |S| K- p) P= C[Txo+ zo + Mz00*], 
| <3'+ 2- |S| K- p)]/?=C[Ig++ y+ (4+ e+ 0], 


= 


| (2 2* [S| K- p) ?=C[Ie-+ Uy-C 4 Ilp- 2] G3? 
where 


Gy oa Ge : 
ay Sn Kans Ak aa UG aD 
G (= ay G. a (G, ) 


If Jax=Jxx, as is easily seen from eq. (2:15), the process K-+p »+* +7 is forbid- 
den. It will be noticed that this is independent of the assumption (4). As we shall 
see later, Tiomno’s theory of Hyperons and K-mesons'’) is not satisfactory, because the 
experimentally established reaction K~+p—>+*+77~ is also forbidden for the same reason 
mentioned above. 

Now, we see from eq. (3-3) that there remains essentially one parameter ©, as far 
as we examine the branching ratios. This seems satisfactory for testing the global sym- 
metry. We first adjust the parameter £ to minimize the process K~+p—1°+7°, which 
is experimentally known to be very rare. Results are summarized in Table I. 

It is known from experiment that the branching ratios are about 2:1:1:0.25.% 
From the Table I, we see that apart from the qualitative feature that the process K~ + 
p—>+~+7" is dominant, the results are not very good. Next, we try to adjust £ so 
as to give an overall fit for the whole processes. Main features are reproduced in Table 
II, including the processes (2-2). 


* E. g. Proceedings of the 7th Rochester Conference. Introductory talk by Cassels. An experimental 
analysis by Russian work shows that at least 03 seems linear up to 300~400 Mev. This is reasonable since 
the virtual states involved have nucleon pairs with high virtual energies and the range of interaction seems 
extremely short. For the same reason, we assume the linear character for the 3, phase shift. 
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Table I 


150 0.177 —OUE14 es 0567. 1.76 1.40 0.402 — 1.080 1.40 
200 0.237 —0.149 0.693 1.58 0.993 0.490 Gis) 0.993 
300 0.356 = O12 22. 0.897 1.30 0.562 0.702 eZ 0.562 


yee | Tig | IIIs+ | jes ties | siilee e | yerst: xo; a 
0.286 Oatled MMe. S38 = 3,59... 8.95 0.441 14 dS Glos 
0.355 Puna’ nda bono 476 ~3.25 5.62, -| 0.610 (A ROT 0.6 0,05 
0.437 Ob LS ay 0! oa 0577 | 7.08.1 par O799 | 12 al 25 01s 
0.955 0 0 0.716 | —2.58 Dee cdtG 4h ORS t IA 3 ee OS 
R925 p08 le 0 <I t.ao7 ~2.08 1.41 sini aie Boe care a 
ei rc ean at ne A AP a Lene 3 ROR Laake Neel ec ee Ie eS 


Table II 

E (Mev) c Digs Se 50 Ao ss x0 Ao 
100 0.15 4 1 2.0 0.6 0.85 0.85 1D 
150 | 0.20 4 i DD 0.8 1 1 1.6 
200 0.30 | 42 1 2.8 0.7 0.57 0.57 1.4 
300 0.40 222 il 1.4 0.5 0.25 0.25 0.96 
400 0.40 1.9 1 1.5 0.8 0.26 0.26 1.6 

experiment”) De ee lige nae 02 O Sy Wh WS 6, 


From the Table II, it is seen that the branching ratios for the process K~ +N-—-> 
2'+7 are not inconsistent with the experimental results. However, the processes K~ + 
N-—A°+7 is more than several times larger than the experimental results. 

So far, we have confined ourselves to the s-wave, but for the pion energy E=100 
Mev~200 Mev, the p-wave will contribute considerably owing to the (3-3)-state resonance. 
So we must estimate the contribution from the p-wave part. In this case, eq. (3-1) 


has an additional contribution 


M,,=0 


ie : {2e%*5 sin 0, cos #—e* sin 033 sin @ e*?} (3-4) 

where k,?+//=0 and potm=0; m: K-meson mass. In transferring from the 7z- 
scattering to the K~-absorption, we have assumed a linear dependence on the spatial 
momentum since p-wave is involved. Using eq. (3-4), we have examined the contribu- 
tion from the p-wave and it turned out that the result became still worse than that with 
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the s-wave only, ratios being 1:1:0:4. From this fact, we may conclude that the 
K--meson is absorbed mainly from the s-state. 


§4. Concluding Remarks 


In §2 and § 3, we have examined the branching ratios for the processes (2-1) 
and (2-2), considering that Gell-Mann’s global symmetry may predicts at least qualitative 
features of the experiment. The result shows that the branching ratios for the processes 
K-+N—>3'+7 are fairly well fitted to the experiment for rather high pion energies 
300 Mev~400 Mev. It is known from the phenomenological analysis that the phase 
difference between matrix elements of the [=O and J=1 must be large, i.e. 70°, to be 
fitted to the branching ratios for the processes K- +N-—>+ +7," which corresponds to 
the fact that for rather high pion energies we could fit them well. The situation is not 


so well for the processes K- + N— A°+7 as the processes K-+N— ’+7. In general 


it seems difficult to forbid the process K~+p— A°+z=", allowing the process K-+p— 
a°4 7°, 

In this respect, Tiomno’s theory of Hyperons and K-mesons is also unsatisfactory, 
as already mentioned. In his theory, both J, and J,’ (see his paper”) are conserved in 
the processes (2:1) and (2-2). But if this is the case, the process K- +p-+2* +77 
is forbidden, because J;’ is not conserved. And his theory predicts the branching ratios 
S7~:3*:23°: A as 2:0:1:1, which are not consistent with experimental results. 

As a final remark, let us estimate the magnitude of the K-meson coupling constant. 
For this purpose, we must know the absolute value of the cross section for each process. 
Experiments show that for the processes K~-+p-—>+*+7* the energy dependence of the 
cross section obeys the 1/v law." At E=40 Mev, o* is about 10 mb. From this 
and Table II, taking the value of the outgoing pion energy E as 100 Mev, we obtain 
the result 


(Jax—Jsx) 47z7~4, 


It will be noticed that the magnitudes of the K-couplings are not separately obtained 


from these treatments. If g/,x/47 is of the order of about unity, as Gell-Mann pointed 
out, this result shows that 


192/42 < 10. 


We could not, of course, decide only from these results whether Gell-Mann’s global 
symmetry is valid or not. Many attempts will further be needed to decide the validity 
of the idea of global symmetry from various points of view. 


The author wishes to express his sincere thanks to Dr. H. Miyazawa for his con- 
tinuous guidance and encouragement. 
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Many authors have shown that the universal V-A coupling of the Fermi interaction is not incon- 
sistent with various experimental results. Their assumptions are more or less ad hoc. In this paper, it 
will be shown phenomenologically that the universal V-A coupling is a natural and almost unique conse- 
quence of our previously mentioned standpoint ; especially, a possibility of S-T-P coupling of §-decay 
interaction is definitely excluded. To lead to this end, we use one additional assumption setting the 
relation between a weak boson-fermion interaction and our fundamental interaction. 


$1. Introduction 


Our basic standpoint in I” was to admit the common character of the interactions 
belonging to each stratum of interaction strength. This standpoint was expressed by four 
assumptions : 

T) Species v, e, 4, n, p and A are assumed to be elementary or different particles 
and other stronglies are composed of n, p, A and their antiparticles. 

II) The original interactions are of the Fermi type without derivative. 

III) All interactions are invariant under the individual /’, transformation.*** 

IV) The interaction 9((, 0 ¢s) (, 0’ ¢’4) is strong if a=b and c=d, and weak 
if otherwise. Such combination as a=b and c=d, if they contain leptons, are missing 
or much weaker than the weak interactions for some unknown reasons. 

The most noticeable implication of above four assumptions is the parity nonconser- 
vation in all weak interactions. As is well known, it was established that the A-decay 
displays the parity nonconservation.. Because this process contains no neutrino, the 


parity nonconservation is not the characteristic property of the process with neutrino but 


* The preceding article, published in Prog. Theor. Phys. 19 (1958), 475, is quoted as I. 
** The preliminary account of this paper was published in Prog. Theor. Phys. 19 (1958), 589 (L). 


*** The individual I’; transformation of “a”-field is defined as 
Pa>Na® 5 Yar oo" das Po >e™ Yey** 


mq >— May; mn-->mp, Mo—>Mey*** 


where 7's are phase factors: | |2=1. 


Phenomenological Model on the Interactions of Elementary Particles, II 127 


is common in all weak processes. The result of this experiment ensures the validity of 
our standpoint. 

Recent experiments give some information on the asymmetry and polarization of 
particles emitted in decay processes. Many authors obtained the consistency of the assum p- 
tion of universal V—A coupling of the Fermi interaction with these experimental results.°” 
But their assumption is ad hoc. In this note, we show phenomenologically that the 
universal V—A coupling* of the Fermi interaction containing leptons, assumed by many 
authors, is necessary from our standpoint. 

From our point of view, the Fermi interaction is the fundamental one and _ others 
should be understood as an apparent form. Since we have no knowledge about the 
structure of a boson as a composite state of fermions, we must assume some relation bet- 
ween the fundamental interaction and others, as will be shown in § 2. 

In § 3, it will be shown that our assumption III) gives mutual relation of phases, 
47s, and only those cases in which 8-decay interaction is of V-A type are allowed but 
all other cases in which /-decay interaction is of S—T-—P type are inconsistent with the 


experimental fact: the preferential backward emission of e* from /4* meson. 


§2. The boson-fermion interaction as an effective form 
of the Fermi interaction 


According to our assumption II) mentioned in § 1, weak boson-fermion interactions 
should be considered to be an effective one induced by fundamental Fermi interactions. 
So, we treat the decay processes containing a boson as follows. The effective transition 
matrix element of decay such as c->a+6 induced by the V—A type Fermiinteraction, 


which is restricted by the assumption III), is assumed to be written as 
M= (27) *0"(k—p—q) (Halp) A£F75)7.0(Q)) 
-O,(k) + Ark, (27%) 
where f~1.” This assumption is tentatively applied to processes without leptons, i.e., 
A-decay process. This form of M is formally equivalent to that of gradient weak boson- 


fermion interaction. The factor A, an effective coupling parameter, may be a function 


of masses of participating particles and has the dimention of length in natural unit. 


§3. On the uniqueness of the universal V-A coupling of 


the Fermi interaction containing leptons 


We have acquired some experimental facts concerning the asymmetry and polariza- 


tion of particles emitted in decay processes. 


* It is to be noticed to which ordering of four fermion fields we refer as the V—A type Fermi interac- 
From the assumption III), it is possible to forbid the mixing V-A and S-T-P among the weak in- 


tion. 
The order of four fermion fields is taken so as to allow such a separation. 


teractions under consideration. 
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Especially important facts are the following : 

a) The e~ (e*) emitted from a polarized nucleus is almost completely antipolarized 
(polarized) along the direction of its motion.” 

b) The e* emitted in =*—p*—e* chain predominates backwards relative to the 
direction of motion of #2* meson.” 

b,) The e* emitted in this chain is almost completely polarized along the direction 
of its motion, according to the experimental result obtained by Liverpool group.” 

b,) While, according to the experimental result obtained by Columbia group, the 
direction of polarization of e* is contrary to that obtained by Liverpool group.” 

c) The polarization of 4 emitted from =* is not yet known. 

d) The similar situation as **—z*—e* chain is also realized in KX —#* —e™ 
chain.'” 

e) Most = mesons are emitted preferentially in the direction of (Pi, Pa) in 
A->p+7- decay, and the asymmetry coefficient @ jeer ea (0.44 ~0.52) .” 

Using our assumptions, we will investigate which combination of weak interactions 
is consistent with the experiments mentioned above. It will be seen that the interaction 
of P-decay, v-capture and Ay (or e~) +»+> can be simultaneously of either V—A or 
S-T-P type. If we assume that these interactions are of V—A type and admit the 
experimental facts a), the possible combinations of weak interactions are limited to the 


following cases. (See Appendix 1.) 


i) fGoAt7s)i> $n) Gel—Ts)7e We): 
ii) FP.1—Ts) 77> Pu) (A 7’ Wy), 
iii) FfGpA1s) ip Yn) Prl—7)i> fr)s 
iv) FG x17) 7, Ha) Gop Te! Gy) he 
v) Fr Ati) Ya) (Oe A—7s)7> Hr’) 
vi) FoAETd 7 Yr) Gul—Ts)%> Hr); ’ 
i), ii) and iii) are as same as those in (3-1), and 
iv) eis (Ga(L£75)0 yn) (Pp 0 Gp), 
v) Fro LET) Tp Pu) Gel —7s) 7% hy’), 8 
vi) FGp(L FT) Tp Pa) Gr —Ts)7> Pr), 


i), iv) and v) are as same as those in (3-1), and 


v IG A=1) 0 $40) GOW, 

tee ; ‘ yO a. 3-3 
iii) Fr), Yn) Gull +red7, ry), a 
vi) fY, (1 a6 vs) lp Wa) (by er ae 1s) Tp h,) > 


and 
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i), ii) and iii) are as same as those in (3-3), 
v) is as same as that in (3-2), and | 


 ppe as 3-4 
i) DLE FI) 0 Gu) Gy 0p) | (3-4) 


a, FG, 7s) lp Ws) (f.(1 +75) Tp py), 
where the double signs are taken with the same order, d/ =, 00 by (GC by (Ct 
=i? = be of. ‘ 
Cera Gals at), i> =7> or 7,7, and ”, is the annihilation operator of 
in. 


Here we assume that the lepton number conservation law is valid. Then, if 
= ¢,, we should assign the same lepton number to ~ and e~, say, +1, and if ¢,/=¢,', 
we should assign the same lepton number to ~~ and e* say, +1. The expressions 
(3-1) to (3-4) are irrespective of the lepton number assignment. (See Appendix 1.) 

Apparently, Cases (3-1) and (3-2) (or Cases (3-3) and (3-4)) give the same 
prediction concerning *(or K) —f—e chain, which is consistent with b), b,) and d). 
(About these arguments, see refs. 11), 1) and 3).) The behavior of * (or K*) 
—f~* —e* chain predicted by (3-1) (or (3-2)) and (3-3) (or (3-4)) where ¢,/=¥¢, 
in the latter case, are schematically illustrated in Fig. 2, Case II,b (see also Fig. 1 
in ref. 12)), and Fig. 4, Case JJI,5 in I respectively. 

In order to explain the largeness of the asymmetry coefficient @ in A>p+77 decay, 
CA p) (pn) is desired to be of V—A type.* Thus, iv) in (3-1) and (3-2) should be 
respectively of the form 

fGxAL+Ts) 1, Yn) GoAtirs)7, Pp) 
and 
Fa FIs) Pn) (Pp (1 £75) Hp) ** 
But, as we can see later, the forms of iv) have nothing to do with our conclusion. 

If we take S-T—P as the coupling type of /3-decay of nucleons and admit the fact 
a), four cases that give the consistent phase relation are obtained. (See Appendix 2.) 
Even if we do not know the stumbling block inherent in the S~T—P coupling (which 
has the difficulty to account for the suppression of /3-decay of bosons), we can discover 
that this coupling should be discarded. Because, the S~T—P coupling gives the forward 
emission of e* in =* (or K*) —p*—e* chain, which contradicts definitely with the well 
established fact b). (See Figs. 2 and 4 in I.) Therefore, we can assert that the p- 
decay interaction of nucleons is of V—A type. This is consistent with the experimental 
results obtained by Herrmannsfeldt, et al.” and Goldhaber, et al.””, but inconsistent with 


the result of He® recoil experiment by Rustad and Ruby." 


* Using the gradient boson-fermion coupling, we get the asymmetry coeficient aeq7z(A°>p+n27) ~ 
0.9, and on the other hand, in the case of the direct coupling, Qcat~0.19 ) 


** Assuming x fn (1475) 01d) (Pp O,/ ¢,) to have the same form as iv), we can get the required 
t=1 
forms of iv) in (3-1) and (3-2). 
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One of the differences among (3-1) to (3-4) is found in value, i.e., for 
(3-1) and (3-2) p=3/4, and, for (3-3) and (3-4) where p=, p= (9/16) 
(1—a@), where @ is the asymmetry coefficient for 44” decay.” The small discrepancy 
of the experimental value of ’” from 3/4 in naturally explained by the latter.’ However, 
it may be said that (3-1) and (3-2) are more preferable than (3-3) and (3-4) 
because of the following reasons. First, the universality of the Fermi interaction contain- 
ing leptons is realized in (3-1) and (3-2). Secondly, if we assign the same lepton 
number to /~ and e* (which means ¢,/=¥¢),°) in order to forbid many unwanted pro- 
cesses such as pt—ettet+te, ptoet+y and p*+H7(—e*+H(, etc., then (3-3) 
and (3-4) are definitely excluded, because ~=0. This assignment, combined with the 
assumption IV), implies that all leptonic processes must be accompanied with neutrino.* 
The experiment to measure the polarization of 4#~ emitted by 7* or K™ and the inves- 
tigation of #--capture interaction” are important to check the validity of the universal 
V-A. coupling. 


§4. Concluding remarks 


From our discussion it was concluded phenomenologically that the type of §-decay 
interaction must be of V—A and that the universal V—A coupling of the Fermi interac- 
tion containing leptons is a natural consequence. Our result is consistent with the assump- 
tion made by Feynman and Gell-Mann,” Sakurai? and Marshak and Sudershan*’, and the 
analysis by H. Umezawa, et al.’” According to our standpoint, the gradient weak boson- 
fermion interaction should be an apparent form of the V—A type Fermi interaction. 

Assigning the same lepton number to “ and e* and admitting our assumptions, we 
are led to conclude that leptonic processes are always accompanied with neutrino. There- 
fore we can get the same selection rule as pointed out by Feynman and Gell-Mann* and 
Oneda :*” 4Q (baryonic charge) = +1 when JS=+1. 

We are interested in the process dependency of effective coupling constants when 
Fermi interactions are reduced to effective boson-fermion interactions. The A-factor in 
(2-1) for various decay processes which really occur should be of the same order (107° 
~10~'). Observing available decay processes, we can find some tendency that the larger 
the Q-value is, the smaller the decay probability becomes. This Q-value dependence 
seems to us to be the character of A-factor in (2—1). 


One of the authors, K. Fujii, wishes to thank the Yukawa Yomiuri Fund for a 
research grant. 
Appendix 1 


In (3-1) to (3-4), if Qontr.g, then o,°>—7* 7, hr! 7, o% The 
expressions of (3-1) are obtained as follows. Accepting i) in (3-1), 


* While, if we assign the same lepton number to “~ and e~, only for (3—1) and (3—2) we are allowed 
to assign the lepton number to neutrinos in the sense of the two component theory.1®) 
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7 (a2 paev! )==7 5 yeast ts 1 for a=p and n, 


(A-1) 
n (a: pnev')=+1 for a=e, v, “and A. 

As //-capture and /-decay really occurs, 

(4: pnev) =-+1 for a=e and A, (A-2) 
and 

n(a: epvy’)=-+1 for a=n, p and A. (A-3) 
From (A-1) and (A-3), it follows 

(a: pnpy) = F1 for a=p and n. (A-4) 


Therefore, irrespective of the type of //-decay interaction, the /-capture interaction is of 
V-—A and the double signs must be of the same order as i) in (3-1). From (A-1) 
and (A-2), it follows 


(a: epyy')y=+1 for a=e and A. (A-5) 


Assuming that the type of /-decay interaction is of V—A, we get the expression ii) in 
(3-1). Then, 


n(a: epyy’) =+1 for all a’s. (A-6) 
From (A-1) and (A-6), the following is deduced beside (A-2) and (A-4), 
4 (a: pnpev)=+1 for a=v and py. 


Then iii) in (3-1) must be chosen. In the same way, we can get iv) to vi) in (3-1). 


Appendix 2 


If S-T-P coupling is taken as the [-decay or -capture interaction of nucleons, the 


following cases ate possible : 
i) > f(Pp(A—7s)0 Yn) (Pe(1—7s)0 fy’), 
ii) SG. (1—75)0 dy) (f, 0 fy’), 

iti) ae (1—7s) Pn) (Pu (1 +75) Py), 


(A-7) 
iv) ee f(Px—7s)0 Pn) (Pp 9 Pr) 
v) fGoA+Ts) a) (Ge(1—7s) fr’), | 
vi) fA +75) 0 $a) (Pu(1 +75) fy), 
i), ii) and iii) are as same as those in (A~7), and 
iv) fn ca +75) Ton) CA i b,) ? (A-8) 


v) SG 1-75) 0H) God —75)0 $!), | 
vi) fGp AT) a) GuA+70) $y), 


132 K. Fujii and K. Iwata 


i), iv) and v) are as same as those in (A-7), and 


ii) fGeA—-7s) 7p Pu) Pr Ty! Pr)» (o0) 
iii) > fGp(1—7s)0 Yn) Gu —7s)0 $x) 
vi) fGrA+7s) $0) GrA—-Ts) $y) 
and 
i), ii) and iii) are as same as those in (A-9) and \ 
iv) and v) as same as those in (A-8), and (A-10) 
vi) SfG,A—7)0 ¢s) Gu —7s)0 Hy). 


The reason for that (1+7,;) disappears from the nucleon part of /-decay interaction is 
that the tensor coupling is necessary as far as the allowed Gamow-Teller transition in /- 
decay occurs. 

Any choice from (A-7) to (A-10) gives inconsistent results with experiment. 
Because, if the spiralities of e~ and y* emitted respectively from #* and = (or K~) 
are determined according to the corresponding interaction Hamiltonians, we see that e~ 


from {{* mesons should be emitted preferentially forward with respect to the direction of 


e 9 
motion of “4* mesons.”’” 
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Charge-transfer energies or electron exchange energies ate calculated between conjugated chain 
molecules and between benzene molecules which are formed into molecular compounds from the same 
species and called self-complexes. The essential nature of these energies is that these are largest in 
the displaced position that one molecule is laterally translated by some bond lengths to the other, 
and values of which cannot be negligible relative to the dispersion energies at the normal separation 
distance at which the molecular compounds are formed. The results obtained explain the observed 


phenomena qualitatively. 


$1. Introduction 


It is well known that the dispersion forces are most important and dominant among 
the inter-molecular forces. The dispersion energies between large chain molecules have 
been calculated by Coulson and Davies” with the method of LCAO MO (molecular 
orbital constructed by a linear combination of atomic orbitals) or by Haugh and Hirschfelder” 
with the Free Electron method, which are both oriented from London’s treatment’? of 
dipole-dipole interaction. They showed not only the long range character of the disper- 
sion forces but also certain interesting directional properties of these inter-molecular fields 
of forces. However, at the distance where large chain molecules form the crystal, the 
charge-transfer forces may become effective. Mulliken” has already suggested the importance 
of this force in electron donor-accepter system and explained the charge transfer spectra 
of many molecular compounds. Also he pointed the laterally displaced systems observed 
in benzene crystal or graphite are energetically stable from the view-points of the symmetry 
natures of electronic wave functions. Recently Suzuki,” considering the charge-transfer 
forces, treated Quinhydrone which is the molecular compound between Quinone and 
Hydroquinone, and obtained the reasonable stabilization energy of this complex, with the 
use of the estimated resonance integral from the charge-transfer spectra. 

I shall show that the charge-transfer forces between conjugated chain molecules and 
those between benzene molecules give some part of inter-molecular forces at rather short 


distance, and there occurs the lateral translation of one molecule relative to the other, 


when molecular compounds are formed. But in the following treatment the limitations 


are that the same kind of two molecules approach each other in the manner that the 


two molecular planes are parallel and the overlaps between the atomic orbitals of distinct 


molecules are as large as possible. 
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§2. General formulation 


The ground state N of any molecular compound AB is expressed as 
P y=a¥ +5,b,F,. (1) 
In eq. (1), Yo is a no-bond wave function and has the form 
P= (A) 7, (B), (2) 
where ¥,(A) and ¥,(B) are the ground state wave functions of A and B molecule 


respectively. ,; is a dative wave function corresponding to the transfer of one or more 


electrons from B to A, or A to B. 
Pv =F (A*—B-)4+27,(A —B*) (3) 
The coefficients a and 5, are related neglecting overlap integrals, 
a+2,b7=1. (4) 


Considering energy relations, the second order perturbation theory will give an adequate 


approximation for the loose complex, 


ee \ Fy HE yloxsE— 2 (5) 


where 
E,= | HY de, E,= | YAY ,dv, 


p 


H,,= \ Uv HY dv. 


H is the Hamiltonian operator for the entire set of nuclei and electrons. , is equal 
to the sum of the separated energies of A and B in the ground state, and E, has a 
similar meaning in the charge transfer state. The stabilization energy owing to the 
charge-transfer is 

Hy 


AE=~ — 
E,—E, 


(6) 

An essential requirement in eq. (1) is that V, shall be of the same group species 
as ,, otherwise H,, is zero and there is no resonance. This means usually that Y, 
must be of the same spin type as V,, and of the same orbital species to the overall 
symmetry of the complex as a whole, under the group-theoretical classification correspond- 
ing to the complex. Now suppose that V,(A) and %,(B) in eq. (2) are both totally 
symmetrical singlet state and /, and ¥, must be of this same type. 

In the following LCAO MO approximation is used. 


§ 3. Between conjugated chain molecules 


For the conjugated chain molecules which have 2m carbon atoms, Coulson® has 
derived the general formula: the molecular orbitals 
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2m 


o, = a Cha raN (A) 


where 


= I) ; NAT 
ON Pye om 


and corresponding energies, relative to the energy of an isolated 2pz carbon atomic orbital 
X, ate 


<<, =2.0,.. ccs ae (8) 
2m+1 
where §,,, which is negative, is the resonance integral between adjacent atoms. In these 
formula, the overlap integrals are neglected, between distinct atomic orbitals. Such neglect 
is known to make very little difference to the energy in the ground state. 

The ground state of each molecule is obtained by filling each of the lower m orbitals 
with a pair of electrons according to the Pauli principle. The charge-transfer states are 
given by configurations in which one or more electrons are translated from the lower 
bonding orbitals of one molecule to the upper anti-bonding orbitals of the other molecule. 
From these configurations only the most effective ones are selected, under the symmetry 


consideration, that are shown in Fig. 1. 


Case a Case b Case ¢ Case d 
m—2 Oo —— a ects a Fa a et 
a ea re ag ee 
be -e-e- ——-e- -—-e-e- eo 
nm-1 —@-e- --oe eo =e —~-e-  -¢-o- ooe- -3—< 
1 <~<o -3-o- -e-o- oo eo -°©o -@ee- -<o- 
A B A B A B A B 


Fig. 1. Charge-transfer states between conjugated chain molecules. 


Only Case a and Case b shall be considered, because Case c and Case d give only 
negligibly small contributions, proportional to the square of the overlap integral. 


Then the stabilization energy of the complex, relative to the separated molecules is 


= H' (a)? ek aan, 
Cm $2 (B) =C, (A) Coil (B) aa Pe (A) 
HAGIh eelitetiiayys 2 (Sc, 


“: (9) 
Em? (A) —&,,(B) Emel (A) = ees (B) 


“| 


where H’(a) and H’(b) are interaction energies of Case a and Case b respectively. The 


denominators of eq. (9) 
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Em+2— Em = Em Soot 


=— 28, {cos —™ — 2+ cos" ——_ a} (10) 
2m+1 2m+1 

are easily verified. Next H’(a) and H’(b) 

must be expressed in terms of atomic orbitals (7) 
and their coefficients in each molecular orbitals. é 
Now suppose, as in Fig. 2, that the two 

molecules A and B approach in such a manner (i) 
that their molecular planes are pararell with ¥ 
each other, and then the interactions between 
electrons of 7Z-orbitals are o-type. It is 
obvious that such an approach obtains the 
greatest interaction energies. Here, by example, the explicit expressions of H ‘(a) and 


3 


Co Co 
Cu Cor us 
eS Car 
Se 


Fig. 2. An approach of two conjugated chain 


molecules. 


H’'(b) corresponding to Fig. 2 are given 
2m 2m 
H! (a) = | S¥ ena (A)Z (AH! 3 ns2u (B)Z, (B) de 


: Amo (a) [0]+ Asn 1 (a){1]+ verse Mee ot (a) | 2m— 1], (11) 


where H’ is the perturbing potential. 


[«]= \ 1 (A) H'%,.,(B) dv (12) 
and. 
Asmnirt (a) = 2 {mu (A) Coch sae “I> Coats x (A) Cm+2.y (B) 
= Si ein Pe) eS (m+2)(u+x) kee (4#+4) wale (m+2)m _) 
wal 2n+1 2m+1 2m+1 2m+-1 
a (=1)*"Jcos 7+ sin at he [sin a 
2m+1 2m+1 2m-+-1 
SPOR a m(«+1) 4 [sin ade = | iy 
2m+1 2m+1 2m+1 ) (13) 


The similar equation for H’(b) is obtained 
1 
Tf (b) = 2 Ay, 9 (b) [0]+ Aon» (b)[1]+ faee c2 + Asn, 2m—1 (b) [2m —1] (14) 
but 


| mK ; 
Asm,» (b) = — jos 7 + sin- ones 7 |sin i 7 


2m+1 2m+1 2m+1 
= tos sin ETD foin ot “t. (15) 
2m-+1 2m+1 2m+1 
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The charge-transfer energies between the chain. 1 2 3 Din 

molecules of the same kind will depend on 

the size of molecules and on their mutual ip4 

coordinates. The size of molecules is simply 1 2 3 5 on 
iM le 

determined by 2m, the numbers of carbon amo! 

atoms. For the purpose of numerical evalua- normal separation of chains. 

tion, the interatomic distance within each «: lateral translation of one chain relative to 

molecule is taken as 1.4 A, which is the value the other measured by the unit of one 


bond length. 


appropriate to a conjugated carbon chain. 
Fig. 3. Mutual coordinates between two 


The mutual coordinates used are shown in 


molecules. 
Fist 3: 
The resonance integral |«] is separated into o-type and 7-type ones, 
[«]=cos*4|«],,+sin?4[«],. (16)- 


[«],, and [x], are estimated in unit of 9, and §,, which are respectively resonance 
integrals of 2pz and 2po atomic orbitals at the distance 1.4 A in the field of molecular 
chain. Mulliken estimated 9, as —40 kcal/mol from the spectroscopic data, while 3, 
is here roughly supposed to be —100 kcal/mol, considering this value to be proportional 
to the dissociation energy of the single bond between carbon atoms. Moreover it is 
assumed that these values depend upon the distance in the same manner as overlap 
integrals between corresponding atomic orbitals do. 


As an example, butadiene shall be considered. 
(€) 
H’=4A,,[0]+ A, o{1]+Ay2[2]+ A,3[3]. 


The same formula stands for a and b. 


(2) «=1 

H’=4A,,[0]+4(Ayo+ 412) [1] +3 (Aaa + 44s) [2] +3 Asol3) +3 41314] 
(3) eo 

H’=44,.[0] +4 (Ayr + Ags) [1] +4 4s0[2] +4 41113]4+-4$ 41214] +3 Acc[5] 
ve ak=3 


H’=1(A,3[0]+ Ayo[1] + 4i1[2] + 41,03] + 41114] + 412[5]+ i316) 


In practical calculations, the smaller terms than [4] are neglected. For larger molecules, 
similar equations are easily obtained. Results of calculation are tabulated in the following. 

From the foregoing calculation, we can know (1) the charge-transfer energy exceed- 
ingly decreases as normal separation of molecular planes increases. That is, the dispersion 
energy is the function of the reciprocal power of distance, but the charge-transfer energy 
is the function of the exponent of distance. Therefore the charge-transfer force is 


notable at rather short distance, where the valence force is about to play its part. (2) 
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Table 1. Am.~ (2): upper, 42m,x (b) : lower 


0 0 0 0 
: 0 0 0 0 
1.1186 1.5267 1.7060 1.7993 
Z 1.1186 1.5267 1.7060 1.7993 
= 11186 —1.2240 —1.1140 —0.9851 
1.1186 1.2240 1.1140 0.9851 
— 1.2240 —2.0913 —2.6354 
. —1.2240 2.0913 2.6354 
1.5267 2.0913 2.187 
> 1 467 —2.0913 —2.1879 
1.1140 2.1879 
° 1.1140 2.1879 
y ei, —2.6354 
1.7060 2.6354 
—0.9851 
8 —0.9851 
1.7993 
4 — 1.7993 
Table 2. Charge-transfer energies between conjugated chain molecules (kcal/mol) 
2m k oe 4A 8A 15 A 
0 1.25 0.1307 0.59379 0.149-4 
1 29.7 0.2667! 0.138-S 0.181734 
4 2 | 0.308 0.30778 0.252-4 
3 0.308 0.30478 0.2357 
* 0.935 0.34171 0.10672 
0 3.06 0.27471 0.13878 0.292724 
oI 72.2 0.181 0.28978 0.296724 
2 0.594 0.44578 0.234-%4 
6 3 0.288 0.16478 0.911735 
4 0.287 0.15478 0.467725 
5 0,593 0.43978 0.20874 
* ) 2.62 0.140 0.55872 
8 4 | 1.06 | 
10 a | AR] 
: | 7.95 


0.130-1=0,130 X 1071 


*: Dispersion energies obtained by Coulson and Davies. 
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The charge-transfer force is most effective at the definite lateral translation of one 
molecule to the other, which is characteristic of the sort of the molecule. Such circum- 
stance is known from 4A,,,, values of Table 1, that is to say, the charge-transfer force 
is strongest at the relative position where the absolute value of A.» is highest. How- 
ever, other factors, e.g. the steric hindrance and the repulsions of attached hydrogens, 
etc., are considered together. 


$4, Between benzene molecules 


It is known that benzene molecules are constructed into a crystalline state as shown 
in Fig. 4. In Fig. 4 (a) the full line and dotted line both show benzene molecular 
planes which are perpendicular to paper, and the former is translated by about a half 
width of molecular plane relatively to the latter. The similar structure is observed in 
graphite. A pair of A and B molecules in the frame of Fig. 4 (a) is taken out as a 
unit and will be discussed, which is looked 


down as Fig. 4 (b) generally and at the 


(b) 
case of graphite as Fig. 4 (c). Coulson” 
has pointed out: If the dispersion force is 
dominant among the inter-molecular forces i = 


between two benzene molecules, the less- 


Nes er 4 


directional nature of this force makes it 


possible for benzene molecules to pack together wt 


¢ 
Le 


as closely as possible. And any molecule , o dea) 
: : : Fig. 4. Crystalline state of benzene molecules. 
(say the full lined one of Fig. 4 (a)) is 


surrounded by twelve neighbours, four full lined ones at the same horizontal level, and 


four dotted lined ones above and below. This illustration is satisfactory, but cannot deny 


dy 6, ba 


1 =i — 
i t Z axis is perpendi- 
fer 31 Tecate 1 ay 1 cular to paper 
t —1 
2 —— 20 () a! 
5 Vile 
tif oe —1 1 <I 1 
db, ds d, 


Fig. 5. Molecular orbitals of benzene. 
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other sources for the structure of benzene crystal. The charge-transfer forces between 
benzene molecules support this structure energetically, which is not expected from the 
view-point of the dispersion forces. 

The molecular orbitals of benzene obtained by the linear combination of six 2pz 
carbon atomic orbitals are shown in Fig. 5. The numbers on vertexes of the hexagon 
are the coefficients of linear combination respectively to each 2p% carbon atomic orbital, 
and the numbers inside are normalization factors, where, needless to say, overlap integrals 
are neglected. The symmetry natures of these orbitals are tabulated in Table 3 referring 
to the coordinate system in Fig. 5. (x) is reflection in the xz plane and o(y) is 


reflection in the yz plane, and i is inversion in the center of symmetry, but the signs 


Table 3. Symmetry natures of benzene molecular orbitals 


| o 1 | $-) G2 $-2 Os 
ee oll ED Sen ie : 

a(x) | - ae = § 2 rt 

o(y) | + = 4 s J 7 


| + | = = + + = 


in i are taken reversely, considering A—B complex shall be formed. A and B benzene 
molecules in the non-bonded state of eq. (2) fill 0, and 04, orbitals with two electrons 
respectively, according to the Pauli principle. In the charge-transfer state of eq. (3), 
one or more electrons in the occupied orbitals of one molecule translate into the vacant 
orbitals of the other molecules. As the case of the chain molecule in which only the 
one electron translation is considered, the relative position of A and B molecules is 
restricted under the symmetry natures of molecular orbitals. That is, only Case (d) in 
Fig. 6 corresponding to Fig. 4 (b) is allowed, and others are forbidden. For example, 
considering the charge-transfer from A to B in Case (b) in Fig. 6, the electron donating 
orbital of A and the electron accepting orbital of B share the same nature in o(x), 
but are different in i, therefore the both orbitals cannot be combined and the charge- 
transfer is forbidden. Only in Case (d) in Fig. 6 in which the complex has the sym- 
metry of i only, the above mentioned two orbitals keep the same nature of symmetry , 


Case (a) Case (b) Case (c) Case (d) 
ann 
\ \ =, 
4 > \ 
/ \ > ? 
=f \ / \ / 
ee Never walt 


Fig. 6. Relative positions between A and B benzene molecules. 


and the charge-transfer is allowed. The charge-transfer states corresponding to / (A*— 
B~) are shown in Fig. 7, under symmetry consideration. 


Then the stabilization energies of complex, relative to the separated molecules are 


Charge-Transfer Forces in Molecular Compounds 141 


Case a Case b Case c Case d 
3 Sees Sse eee ees 
-20= == Sac: ond —S _ =—- =o 
ee a 2 en no 
0 —e—@= = -8— —e—-  _-e-0- -ee- -2-e- —~-6— 
A B A B A B A B 


Fig. 7. Charge transfer states corresponding to Y(A*+—B-) between benzene molecules. 


Table 4. Charge transfer energies between benzene molecules (kcal/mol) 


oe CS 2) Come 


\ 

») 

/ \ 4 
ae! s/ Nw TSt/ 
0.632~° 0.9457! 
co") GN oe 

\ \ \ 

) : ) 72 
w id “ Ke \ y, 
Noe Sy \ ee a, ---~—/ 
0.2767 0.9467 0.172 
1 1 1 
a—, Lar 2s, a ) (Point 
eh CA, 
ny lees ‘ ( > 
\ \ / as 
¢ y ¢ ) / \ 
0.6127" 0.156 a 


0.9617! 


0.632-2=0.632 X 10-2 (kcal/mol) 
(i, j) i: lateral translation along x axis in unit of 1/2, bond length 
j: lateral translation along y axis in unit of 3/2 bond length 


142 S. Aono 


H! (a)? H’ (b)? H'()?__, «Hd? 
a &,(B)—&,(A) | &(B)—€1(A) | &(B)—& (A) — €-2(B) —& (A) 
A ay H’‘(b)? H'(c)? 4 H’(d)? z 
*E(A)—€(B) | &(A)—En(B) | -& (A)—&(B) 2 (A) — 6, (B) ! 
(16) 


The denominators in eq. (16) are all equal to —2.. Numerical calculation is performed 
quite similar to the case of the chain molecule, but the normal separation between two 
molecules is 3.5 A (3.41 A in graphite). The results obtained are shown in Table 4. 

This calculation shows that the charge-transfer energy is largest at the relative 
position where two benzene molecules are looked down as naphthalene molecule. Though 
the absolute values of the charge-transfer energies will change corresponding to the values 
of 8, and f,, it is hardly considered that their relative values will vary. Therefore the 
above mentioned Mulliken’s speculation corresponding to the displaced structure of benzene 
crystal or observed in graphite, is not fully satisfied by this calculation. While the 
dispersion energy between two benzene molecules was estimated as 1.50 kcal/mol by 
Hirschfelder with the use of Free Electron method, which is the largest value at the 
position where two benzene molecules are situated quite face to face as Case (a) in 
Fig. 6. This value is about ten times the largest value of the charge-transfer energy, 
so the charge-transfer energies give small contribution to the inter-molecular forces between 
benzene molecules. The dispersion forces are essentially caused by the intra-molecular 
phenomena that the dipoles associated with the electronic excitations in one molecule 
interact with the ones of another molecule. On the other hand, the charge-transfer 
forces are oriented from the inter-molecular phenomena, charge transfer or electron exchange 
between two molecules. It is naturally assumed that these two kinds of forces can 
contribute additively to the formation of molecular compounds, because they are characterized 
by perfectly different electronic behaviours. From this assumption one can explain qualita- 
tively the observed displacement in benzene crystal, or in graphite as the energetically 
stable structure, that is the intermediate structure between the perfectly overlapped one 
and the naphthalene type one, when two molecules are looked down. 

$5. Conclusion 

The above treatment is not necessarily satisfactory on the following point: the use 
of semi-empirical orbitals neglecting overlap integrals and the wave functions with single 
product of molecular orbitals, rough estimation of 3 and §,, or quite neglecting the 
effects from o orbitals. Moreover, the zigzag structures of chain molecules are simplified 
into the linear ones, and the repulsive forces by attached hydrogens are not considered. 
Nevertheless the interesting results are obtained, the lateral translation between two 
molecules is quite essential in the charge-transfer state, which is in good agreement with 
the observed phenomena, and the charge-transfer energies are certainly effective at the 


distance of molecular compounds to be formed. It has been usual that if the shape of 
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molecule is large and very far from being spherical as seen in benzene, the packing is 
considered simply geometrical as close-packing, quite similar to the case of the metallic 


crystals or the case of the interaction between the rare gases. Now one may believe that 


the laterally displaced structures observed in benzene crystal, or in graphite, are energetically 


stable ones, considering the charge-transfer forces. 


1) 
2) 
3) 
4) 
5) 
6) 
7) 
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Use is made of the branching ratio, lifetimes, and asymmetries of the angular distributions in 


S-decays to show that the selection rule |4I|=+ is consistent in the decays of S-hyperons as well as 


in the A-decays. The invariance under time reversal is assumed, and the asymmetry parameters, the 
average degrees of polarization of Y and the ratios of the coupling constants are obtained. 


S 1. Introduction 


According to the generally accepted assignment of the isotopic spins to the hyperons 
and #-mesons, the decays of these particles always accompany the change JJ, =} in the 
isotopic spin component. The change |4JI) in the total isotopic spin has been discussed 
by many authors” on the basis of parity conservation or non-conservation with the as- 
sumption of time reversal invariance (T-invariance).') The recent experiments on the 
angular asymmetries in J-”) and +'-decays” seem to indicate violations of parity conservation 
and charge conjugation invariance” in the hyperon decays, whereas no evidence of a 
violation of T-invariance has so far been presented. 


In the A- and @-decays the observed ratios of the decay probabilities, 
w(A->n+7°) /w(A) =0.32 +0.05,” 
w (G*—7* + 2°) /w(0,°) = (1.80 £0.23) X 1073, © 
and 
w (4, 27°) /w(0,°) =0.14 + 0.06", 


suggest that the selection rule |JI|==} is dominant in these decays, while the contri- 
butions from |4I|=3/2 and 5/2 are about 10% of the former in the case of (-decays," 

In this paper we emphasize that the analysis of the observed asymmetries together with 
the branching ratio and the lifetimes of the various modes of +'-decays, can serve to 
determine the kinematics of the decays and the polarization of +'-hyperons before the 
decay, We are thus in a position to examine the |4JI| selection rules to some extent. 


§ 2. Derivation of asymmetry parameters 


Let the s- and p-wave amplitudes of the final nucleon pion systems be (a,, 6.), 


(4, by), and (a_,b_) corresponding to the modes of decays (the spin of +'-hyperons is 
assumed to be 3): 


Isotopic Spin Selection Rule in 2'-Decays 145 


S*—n+n*, Stspt+n and S-3n4n7-, 
The asymmetry parameters for these decays are given by” 


a@;=2Re (4;*b;) / (\a; ); (1) 


4 |b; 


where 
ey = 
respectively, The simple isotopic spin consideration shows that 
a_=(a,4+2'%a,)/A and b.=(6,+2'76,) /A, 
with either 4=1 for |4I|=1/2, or A=—2 for |4I|=3/2.* 


Actually the measurement does not give a’s directly but their products with the 
average polarization degrees P* and BP” of 3'* and 37 :° 


Pp*a,=— (0.36+0.21), 
P*a,=— (0.37 £0.19), (2) 
Pp-a_=— (0.13+0.26). 


However, the ratio of the first two is independent of the unknown quantity P*, and 


from eqs. (1) this ratio is expressed as: 
P* &,/P* a, =s5= {Re (a,*by) / (|49|?+|bo|*) } / {Re (aE 4) / Clas P+ 1641) }- (3) 


The branching ratio of +'*-decay and the ratio of the lifetime of to that of 
<'— are given by 


w(S*—>p+n°) /w(L* n+") = R= (aol? + [by))/ Cae? + [O41 (4) 
and 
t(3*) /c (27) r= {Ia, |? + |b, ?+2 Cao]? + [b01”) 
+2 2 Re(a* a,+5% by) } /# (las |? + [64 1? + | aol? +1 50/?)- (5) 
The observed values of s, R and r are: 
s=1.03 £0.80 from (2) and (3), 
R=0.88 + 0.12” (6) 
r=0.470+ 0.094.” 


The assumption of T-invariance and the neglect of the final state interaction** 
imply that a’s and 6’s are all real. In this case eqs. (3), (4) and (5) determine the 


relative values of these amplitudes and consequently a’s by eqs. (iy. im fetins. Of 5, ik 


* The case of |4I|=5/2 gives branching ratio R which is inconsistent with the observed value. 

** More exactly, it is sufficient for the following discussions to make approximations cos (6,—08)==1, 
but not always cos d,2=1. Here 04 is the scattering phase shift of the pion-nucleon system in the final 
state a= (I; J=4;1). These approximations are correct within the errors of about 10%.) 
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and r. The results of some manipulation show that 
a,=+2€[(1—&)/{(1+5)?—4557} ], 
pias (7) 
a_=a,{(1+2sR)—(1+5)V 2R /s} /#r(1+R), 
where 
€= {(14+2R) —#r(1+R)}/2v 2R. (8) 


The degrees of polarization P+ and P™ can be obtained from (2) and (7). The 
numerical results obtained with the aid of the empirical values (6) are given in Table L. 
The table also lists the results with the simple set : 


sa1, Re) jand ek, (9) 
which are in conformity with the data (6) within the experimental errors. 


Table I. Asymmetry parameters and average polarization degrees 


| AI|=1/2 | 4I|=3/2 
i IS Ae empirical set (6) simple set (9) empirical set (6) simple set (9) 
ay + (0.68-+0.29) +1/% 2 + (0.32+0.35) +1/2v 2 
ay + (0.70+0.28) +I1)/ 2 + (0.33-40.36) +1/2v 2 
a = (0.84-+0.23) 1/42 = (0.49+1.3) +11/8v 2 
P+ =F (0,530.22) = (0,500.30) | +=(1.0+0.6) 
P- + (0.15-+0.31) + (0,180.37) | == (0.13+0.27) 


§ 3. Discussion 


Our test of consistency is to see if |P*|<1, as is required from the definition of 
the polarization degree. From the table we see that the selection rule |JI|=3 is quite 
consistent. The results for the case of |JI|=3/2 contain errors of more than 100%, 
and nothing definite can be said. However, if we could disregard the errors, |@.| and 
|@,| are smaller than the corresponding valus of |Pa'|’s, and |4I|=3/2 could be ruled 
out, 

The analysis of the decays of (-mesons and -hyperons has already shown the 
necessity of the dominant |4J|=4. Our results show that the same selection rule 
|4I|=3 is acceptable for the 3-decays, and that the contribution from |4I| >3 could 
be zero in this case as well as in /-decays, in contrast with the case of @-decays where 
such contribution is not negligible. 

Finally it is interesting to obtain the ratios of the coupling constants from the 
above analysis. If we assume the interaction Hamiltonian for the decay to be either 
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Hy =2* (9.49475) nt* +2* (g,+99'r5) pr® 


+3° (9_+9175)n2- +h. c., (10) 


or 


Fn =2* 7. (94 +9475) 09,0 +347, (Gg + Ge! 7s) pI, 
+3-7,(9_+9"7,)n9,2-+h. cy (11) 


the s- and p-wave amplitudes are directly connected to the coupling constants g’s and 


g’’s as* 


b,/ac=[{ (M—m)*—p"} /{(M+m)?— 2°} (gl */9e*) 


b,/a,= {(M+m) /(M—m)} 
me Mom) ie G70.) 
respectively, where M, m and y are the masses of 3’, nucleon and pion, respectively. 


Table II. Ratios of coupling constants 


| Direct coupling (10) | Derivative coupling (11) 
9+! 194. | + 3.9 +25 | +0.46 3.0 
| 
90/90 | +24 + 4.1 +2.9 “£0.48 
9-’l9- | +18 Sen5.4 | 2222 =F0.64 


Table II gives the results of the calculation with the empirical set (6) for the case of 
|4I|=4. It is to be noted that the assumption of |4I|=—4 definitely excludes the 


“ one-to-one law,’ which seems valid in the weak interactions involving neutrinos. 
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The current and the magnetic moment of the nucleon are calculated, starting from a relativistic 
one nucleon formalism. It is shown that the relativistic corrections are very important for the current 
and the magnetic moment of the nucleon, especially for the vector part of the former and the scalar 
part of the latter, where the effects of the recoils are of the same order as those given by the static 
model. The numerical result, obtained by means of Chew’s method, shows that the neutron current 
and the scalar part of the magnetic moment may be consistent with the experiment by choosing suitably 
the coupling constant (g?/4=0.08 or 0.10) and the cut-off momentum (kmax=5y, 64 or 7). How- 
ever the effect of the nucleon recoil to the vector part of the moment is destructive and gives a 
definitely smaller value than the empirical one. It is also shown that the renormalized charge is 
equal to the unrenormalized charge in our approximation which contains only the terms up to the 


order of 1/M. 


§ 1. Introduction 


There are many meson theoretical approaches about the magnetic moment of the 
nucleon, among which the static model by Miyazawa” gives a fairly reliable result on 
the vector part of the moment. In the usual meson theory, however, it is very difficult 
to explain the scalar part of the moment in the static framework,” unless we introduce 
a strong correlation between the mesons in the cloud” or an effect of the virtual nucleon 
pair” phenomenologically. Then, there are many efforts to explain the magnetic moment 
by considering the effects of nucleon recoil,” isobar states,” nucleon pairs” and strange 


particles.” But it is found that these effects are too small to overcome the above 


discrepancy or rather increase shee 

Since the magnetic moment, especially its scalar part, has its origin in the relativistic 
structure of the Dirac particle, and the effect of the nucleon recoil® has been studied 
only pure kinematically, it is significant to reinvestigate the problem in the relativistic 
framework. Then, we assume that the core nucleon behaves like the Dirac particle 
interacting with the meson field relativistically, and thus take into account the effect of 
the virtual nucleon pairs through the zitterbewegung of the core nucleon. The importance 
of the effect of the virtual pairs has been pointed out by C. Iso in his phenomenological 
argument.” The effects of the strange particles are not considered in this paper, for it 
may be expected to be very small.” Furthermore, we remove the S-wave term which 
the Foldy transform of the 7;;-interaction. It is because the effect of this term 


appears in 
ssed considerably in the pion-nucleon scattering and that suppression cannot 


must be suppre 
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be expected in the usual treatment,” so the careless handling of this term may produce 
an unreliable result. 

In order to define the current and the magnetic moment unambiguously, we consider 
the real nucleon moving in the external electromagnetic field and separate the interaction 
into two parts”: one is the convection current interaction and the other is the magnetic 
one. Such a separation can be easily done by means of the technique developed by 


Suura.'” 


It is possible to renormalize the electric charge in our formalism. However, it is 
shown that the renormalized charge is equal to the unrenormalized charge up to the 
order of 1/M, if the real neutron has no charge. This vanishing of the charge of the 
real neutron is the consequence of the requirement of the gauge invariance and it is 
satisfied only approximately in our case. 

The numerical estimate can be done analogously to Miyazawa,” but we consider 
only the renormalized Born terms for the relativistic correction terms. It is shown that 
the relativistic effect is very important. As is seen in the final section, the nucleon 
current and the scalar part of the moment become consistent with experiment. This 
result was never obtained in the static theory. As to the vector part of the moment, 
the correction is destructive and reduces the value obtained by Miyazawa by about thirty 


per cent as a whole. 


§2. Derivation of current and magnetic moment 


Since we are interested only in the effects of the relativistic recoils, we start from 
the Hamiltonian 


4 


H=—ial + 8M+H,+ ford (x) — Efe _62(x) — ROM. (a) 
2 


H,, is the free Hamiltonian of the meson field, x is the coordinate of the core nucleon 
and $, the vector in the isospin space, represents the meson field. M is the renormalized 
nucleon mass and f, is the unrenormalized coupling constant. The counter part for the 
mass renormalization — 0M is introduced explicitly. The interaction f,/.76(x) involves 
the S-wave pair interaction besides the of coupling term in the static limit. As is well 
known, this pair term produces too strong interaction to explain the S-wave scattering at 
low energy. Then, it would be better to remove this term completely rather than to 
handle it carelessly. We shall see later that the term — (,9f°/2M)¢®(x) should be added 
so that our formalism agrees with Chew-Low’s model in the static limit. The real nucleon 
with momentum p is described by the wave function ¥,,(x), which satisfies 


H®,,(%) = M?+p'¥,,(x), 
PY, (x) =(—iV +P,) ¥, (x) =p? , (x), (2) 


where P and P,, are the total and the meson field momentum, respectively. 
In order to define the current and the magnetic moment of the real nucleon, we 
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introduce the interaction with the external electromagnetic vector potential A, which is 


given by 
H!=—e( 5% )ad(x) +6,| A(s) OXI 9)ads, (3) 


where e, is the unrenormalized electric charge. Let us consider the matrix element 
; ; 
M(p’, p) of (3) between the real nucleon states with momentum p and p’; 


a 


\ Por) |! () | Py (x) pede. (4) 


M (p’, P= oa 


The requirement of the gauge invariance enables us to write M (p’, p) in the form 


M(p’, p)= = oi Yu | i4@) (p+ Pp’) +is-qxA(q)}(*+*) 


+idu-qxA| |u)+-higher order terms in p and p’, (5) 


where |u)> and |u’) represent the bare spinors corresponding to the initial and final 
nucleon state, respectively. A(q) is the Fourier transform of the vector potential ; 


— i [ .-ae 
A(q) =e |e AC te, 


where q=p’—p is the momentum transfer. In (5), the constant e is just the renormalized 
charge (not necessarily equal to e,) and Of is the anomalous magnetic moment of the 
nucleon” in unit of the nuclear magneton. 
The remainder of this section is devoted to find the concrete expressions of e and 
Oj. For this purpose, first we separate the motion of the center of mass in the real 
nucleon state with momentum p. According to Suura,'”’* we can write V,,(x) as 
P', (x) =exp {ipX} ¥,(x), (6) 


where 


X=3{H™, (HH, x}+ | H,(@) @—s)ds]} 7) 


is the coordinate of the center of mass of the system and Y,(x) represents the real 


nucleon state at rest. The matrix element (4) can now be written as 


M(p', p)= ork \ CD (x) |e "* H! (x) ef *| B(x) ) dx. (8) 


Next, we perform a unitary transformation 


P(x) =exp(—iP,x)%,, VO,=0. (9) 


* Our Hamiltonian (1) differs from Suura’s by the term — (fo?/2M)P¢°(*), but the relations (6), 
(7), (8) and (36) in his paper hold without any modification. 
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to eliminate the core nucleon coordinate. Correspondingly, any operator O is transformed 
according to 


OS o= pre x O g tte x 


Especially, an operator F(z) of the meson field is transformed like 


F(z) =F (s—x). (10) 
For example, the Hamiltonian and the center of mass are transformed to 
A= —ior +8M—aP, +H, +fap.79(0) — 6s (0) — 36M, (11) 
2 
X=} (A=, (4H, x} + | H, (=) =ds]}. (12) 


The Schrodinger equation is 
H?,=M9,. (13) 


Applying the transformation (9) to (8) and noticing that exp (—ipX ) -exp(ipx) is 
independent of x, we perform the integration with x in (8) explicitly and obtain 


M(p’, P) =| {—eo( 7 )ad+ eA | (9XxI¢),ds! | %) 


+ @,| a qx A\[9x (sxF9) kds|) 


—ip! (@,| (X—x) —a(— + ad +e,A \ (@XI9)sdst | M) 
+ip@,| {-a(42% )adtea | (6xP¢),ds| (X—x)|@). 


The further reduction can be done by means of (36) in reference 11), i. e., 


= es F i 
(X—x) |) = (R-—2) 9), (15) 
with 
R= | H, (s) sds. (16) 


Thus, we have 
MOP) p)= (Ol —e( AT) jad + (4: (pt p) +ia-qxA)} |) 
+ (G,| eA (PXFo) sds |D.> 


+(@,| “qx. | [9x (ox16) hes) 
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* 


+ (9,| (— ) (Ra) (aA) | 0, 


sp Pol | (FX APS), @ Pt p')ds|%) 


"(| {(Rq), | @XAP9).ds} 1%) 


+ ap Pol (@- P+ PA| GX!) .ds|%). (17) 
Now we shall eliminate the odd Dirac matrices in eq. (13). For this purpose, we 
apply the transformation” : 


| @,) =e | Qo» —=exp an (i0,0P,, +f, 017%) | | Bo)» (18) 


where M,=M—0M is the unrenormalized nucleon mass. The Hamiltonian (11) is 
transformed according to 


PSH He 2H GG asda iG, (3G, A) lees. (19) 


The expansion is characterized by the powers of 1/M,. This expansion is not always 
justified for (f,/M,) is not so small in the ps—ps theory, though one of the authors 
(A. K.) has shown that the higher order correction of 1/M, is not so large to invalidate 
the expansion in the case of the pion nucleon scattering.” Then, we obtain 


Poo oe ) N 
H’=0,M+H_-+ p,( >= 0 \oPrd—0,0M 
Pee padi Bid om) t 2M, Ppa 


fo . § zr { CRA ey ee 20 
auaray, (oP, FA} + hey b (20) 


The leading term of the dotted part is the odd term of order (1/M,)*. To eliminate 
the odd term f,p,(t7/2M,) from (20), we further apply the transformation : 


|b.) =e") exp | 


— fo ora 
re Tr ha @1) 


which does not change the even terms of (20) up to the order of (1/M,)*. Noticing 
that /, is equal to unity for the large components, we have 


H|)=0, 


' ies fo fo ) 
fete ie oS ota oP, teh OM, (22) 
5 2M, 2M, 8M, 


neglecting the terms higher than (1/M,)*. If we neglect the (1/M,)* term in (22), 
it is reduced to the Hamiltonian of Chew-Low’s model with the recoil kinetic energy 


P,2/2M, where the constant 0M has no significance at all, since 0M causes only the 
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shift of the energy and the result of Chew’s model is not affected by this energy shift. 
Applying (18) and (21) to (17) and noticing that 0M is the order of zeroth power 


of M,, we obtain 


mip’, p)=(\| (8 )(22®) {4 PtP) +i9-4x} 


e 


+. qx A| [9x G@XP9) hds— | (@X APO), = PtP) ds 


— (2 ){ (Rg), | Ox AP 9) sds} 


=k alt 2 )(He-o-qx At (LAPS, Rq}. ) 


ui A) —(_%_)(P,.p+p)A\| (@xF¢ ds| 
+( it.) (2x), (Rq) (64) —(S4--) (Pa p+P)A | (@X08)sde ||) 
+0(—.) (23) 
Mie 
where we used the fact that for the true vectors such as P, and (oXVo)dsz, 
(|P2|)=0, «| | (9x ¢)ds|)=0. 
The replacement of M, by M does not alter the result up to the order of (1/M)* in 
C23). 
The final expression of M(p’, p) can be obtained, using the relations : 


(|| @xAP9), (= p+p)ds|)=4(A-ptp) (|| @Xa08)ads|), (24a) 


(|| (Rq), | @xAPé)adsl |)=3qX4¢|(RX | (@XP9).ds+h.e)), 


(24b) 
(|{4P,, Rq}.|)=kqxA(|(RXP,—P,XR)|), (24c) 
(|(tX9)5 Rg) (GA) |)=$qg X 4¢| (7X9); (RX) |], (24d) 


which are derived in the Appendix. Hence, we have 


M(p', p) = -() (A-p+p’) ¢| (22)+2{ (7X 20¢),ds 


+(4) Pa: | @XP#).ds} eat t )iqx AC |—M| [9x @xv¢) hds 


+(4**)o4 (Rx | ox P8)dsth.c.)+ 2 (1% . 


eee ee nn nn nn nnn ne ee ee ee 
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+( al * T=) (RX P. +h.) +(4) (tX¢), (3x R) Dez) 


where the first term is the convection current interaction and the second the magnetic 
moment interaction. 


§ 3. Evaluation of nucleon current and magnetic moment 


Before the evaluation of the expectation values in (25), we note some gross features 
on each term of (25). In the following discussions it is convenient to divide both the 
current and the moment terms in (25) into the scalar part and the vector part in the 
isospin space. That is, 


M(p, p)=—(—"-) (A-p+P') (A +Q)|) 


2M 
—i( 2 qx AC (Hat He) |) (26) 
where 
Q,=1/2, : (27a) 
Q=+( 2.) | Gexarp)sds+(t)Pa | @XP8)ods, (276) 
He. babe hy : 
‘ Z +( y ) 2 +( =) + (RXP,—P, xR), (27¢) 


“es —M [6x (sx 09) dst (Rx | ox Po ideee he ch) 


+7,( < )+( =) < at (7) = (RX P,,+h.c.) 


fo 7d 
+(4-) (tx), (6XR). (274) 


The first terms in (27) and the third term of (27d) are the ones given in the static 
model and the others are the relativistic corrections. The second term in (27d) combined 
with the first term has a clear-cut classical image. If we notice that j,(%) = —e,¢XVo 
is defined in the coordinate system with the origin at the core nucleon, we can write 


the magnetic moment /4,, due to the meson cloud as 
H=hlsXjn@—reds, 


where r, is the coordinate of the core nucleon in the center of mass system and is given 


is 11 
by, according to Suura, ) 


r,=—R/M. 
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Then, we have 
nf aeneG aks pa clad (; zB ast 
p.=(—*-){M| 2x jn(2) dz ul x Tah ) [> 


which is exactly the sum of the first and second terms of (27d). The last term of 
(27c) and the fifth term of (27d) are the corrections corresponding to the core recoil 
discussed by Holladay” in a pure kinematical manner. The second term in (27b) can 


be written in a form 
[ees | j.(=) ds], 


except a certain constant factor, then it can be interpreted as the charge induced by the 
distortion of the current of the field by the nucleon recoil. The other terms are originated 
with the zitterbewegung of the core nucleon and it is difficult to sketch simple physical 
pictures. 

The requirement of the gauge invariance enforces that the real neutron has no 
charge, then the neutron expectation value of the total charge operator Q=Q.+Q, must 


vanish. 
(n|Q|n) = (n|Q.|n> + (n|Q,|n)> =0. (28) 


In the covariant 7,-theory this condition (28) is satisfied exactly in each powers of 
coupling constant fj. In our case, however, it is expected to satisfy the condition only 
approximately. It is because the higher order effects are included partially in our calcu- 
lation. Moreover the neglect of S-wave pair term and the introduction of the cut-off 
procedure may violate the condition (28). If we define the quantity 


‘nlO|n> 
A= ANQin) (29) 
(PIQ\ Pp) 
this gives a measure for the violation of the gauge invariance. It is shown in Table I 
that 4 is satisfactorily small for the reasonable choice of the coupling constant and cut- 
off momentum. In the static model the condition (28) is not satisfied at all. 


If the condition (28) holds, the renormalized proton charge is given by 


e= ey PIQ| p) =2e ( p|Qe| p>- (30) 
As is seen in (27a) Q, is exactly one half up to the order of 1/M, then the renormalized 


charge is equal to the unrenormalized charge in our approximation. 

To calculate the expectation values of each term of (27), we neglect the (1/M,)* 
term in (22). This neglect gives a certain effect of order 1/M to the expectation values 
of f,, but it may not be important numerically. As for the recoil energy P2/2M we 


keep only the diagonal parts in meson occupation numbers specified by angular momentum. 
If we expand the meson fields as usual : 


ox (x) = as (2,) me [ der =e agi m | DPe,t.m (x) > 
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me(x) = 3} (—) (24) (at, — at, Joh 3 
a icp 2 kim ki—m | ¢ nm) s ( 1) 


the final Hamiltonian we use is 


Sn | ax a - k 
A= me o,* Dein Axim + fe ) : oe uh ) Om Co (kim ain dijep,) ) (32) 


l,m ; V 62R k V 20, 
where 
Pram (%) = (2/R) MR fi (kr) Yim (9, ~), (33) 
and 
O;* = O,+ as 5 (34) 
2M 


R is the radius of the normalization sphere and U(k) is the cut-off function in momentum 
space. 

We can reduce the expectation value of each term in (27), which contains the 
meson quantities, to the expectation values of bilinear products of the P-wave meson 
operators, since the real nucleon state |)}) contains only P-wave mesons in the cloud. The 


result of such a reduction is: 


an ME \ [ox (2X09) hds=M 1 (at* af, —agk din tag ag*,+diedzm), (35a) 


km Op 


1 pt Cs 1. — 9 + — + _— 
oy, (R x \ (o x Vo) 3da+ h. Ce) > pa m (dg dna ink Am, i dint cies aie Din diay) ? 
sm 


(35b) 
2 ; i a (yy NE : 
i=l estaP3) d =—+(- \)Sumue (+) wl Gbie 
Se ) ROU eS ) lamar) 
* , d y / » / 
x (hey + in) eX (E anglice ) Othe )), (35c) 
1 ie Ge Nuke oy, Lata ae O(k—#) O(k—K! 35d 
sy Pa | OXPP)ade=i(Z) St (lat XaenhIE-K)IE-K), B54) 
ie 1+7, > = ( ace ST, Ra 35 
( M )o( 2 ) (6) 2 M fo Oe Cin Tim) ( e) 
z apn . 
()(42*) xP. -P.xR) =(S*) Fp Simevaindim, (358) 
2M Zz 2 kim 


where O0(k—k’) means (R/7) 04. Only exception is the last term of (27d), which is 
reduced to a form linear in the meson operator. 

The expectation values of the products can be obtained from Miyazawa’s result"? with 
slight modification coming from the replacement of «, by w,* in our Hamiltonian. These 
values consist of the perturbation term with the renormalized coupling constant and the _ 


integral terms containing higher order effects through the scattering cross section. In 
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evaluating the relativistic correction terms we neglect the latter integral terms, which are 
always 20~30% of the perturbation term. The expectation value of the last term of 
(27d) can also be divided into two terms, the perturbation term and the integral term 
and we neglect the latter term in this case, too. With this approximation, the expecta- 


tion values of each term of (27) are given by 


== 2 )| T y= yg ) 1 ( 4 | 4 6 1 6 
2 hae ey x V ad: 4K, o— K, o——— Kes 
Quad (2) | (ax 209) .ds1=(—Z-) (—) | Ais Ka 


K,°s = Lis] > 


2M 
Oa < (sy) P= \ (8x1) .ds|)= —( Z) a ( amy, a K.*., (36) 
rem tl(gh)teedo=( LE) 2 dead 2) dR 


Hae (Ge) RXP.-P XD (To) 


n= (\(L) (Rx | xP Pade th.c.).l)=—-(-—) +S _ ks, 


4a ge 5s 
a= (|(s) Hata o.|)(3-) + Kes, 
Ha = (|(-) (RX Pa othr a aS , 


om (I( Sep) xa (xR). =—(E)-* (2 Jae, 


Pn =<|(—M) \[ex (5X9) eda] =(_ £ ) 3 (3) MK, 
+($*) G+h), 


for the proton spin-up state, where 


Ky ke (ee Ly, =(£Y® dk) 


= 
J wo w*t J w@*t dk 


I,* and [,* are the double integrals, containing the cross section, i. e., 


i? = | (035 (Kk) = O31 (k) — O43 (&) +04, (A) o%* oo del. 
ie wf (a* +a,*) 2k 


L*es 1 | (Gx (k) emit (k) — O38 (k) — on (k) ))o* Ox Udkdl 


3 
wn. of o7* (aX +07) PR 
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v%* is the relative velocity in the center of mass system. g is the equivalent ps— pv 
coupling constant. The expectation value of 7, or o, is given in two ways,’ correspond- 
ing to eq. (23) or (28) in Miyazawa’s paper : 


A: Qui or Pan =| Sabie (| ts y= = ee ) : K,,4+1,*, (37a) 
2: 2 2 47 T 


Bi Ong or ge Chey ce) ee Meenas OE 
F Quiz 1B (| 2 », «| 2 » 2 Jeo an R ? ( ) 


where the intergal [,* is given by 


L*= d | (G5 (k) — 5 O51 (k) —3 0 (R)) %* [4 dkedl 
1 
ir: 


0, (0,* +0,*)7 RP 
and 
3 ee oe . 
I= F+— | + On, dk (38) 


is the unrenormalized coupling constant. (|o373;|) is also given by Miyazawa, i. e., 


Shera = (ta. (39) 


$4. Numerical results aad discussions 


We calculate the quantities in (36), (37) and (38) with the straight cut of the 
meson momentum, which may produce the ambiguity of a few per cent in the numerical 
values. The results are given in the following tables for the several choice of the coupling 


constant and the cut-off momentum. 


Table Ia. The nucleon charge and the values of 4 for the case of g°/4z=0.08. 


or Qi Total Q A 
PR Qy2 Qos | : | 

kmax et 3 | B ev v ys | B ix B 
Sp | 0.303 0.097 peek) | — 0.055 0.607 0.401 — 0.097 | +0.110 
6u 0.269 0.087 0.507 —0.091 0.685 0.503 —0.156 | —0.003 
7 ps 0.253 0.080 0.668 —0.128 0.783 0.611 —0.221 | —0.100 


Table Ib. The same quantities as the above for the case of g’/47=0.10. 


0.504  —0,009  —0,004 
0.613 | —0.147 | —0.102 


0.768 =0:208 Wi O12 11 
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i) Nucleon current (charge) 

The vector parts of the nuclear charge are given in Tables 1. 

A or B refers to (36a) or (36b), corresponding to the two different ways to 
calculate (|t;|) or ({o3/). Q,, is overestimated slightly for the case A, sic we have 
neglected all o,,; except ,, in evaluating I*’s. This neglect makes I,* definitely large. 


“3 


These values of 4 seem to be satisfactory in the qualitative sense. 
ii) Scalar part of the magnetic moment 


The scalar parts of the magnetic moment are given in Tables II. 


Table Ila. The scalar parts of the magnetic moment for the case of g°/4x=0.08. 


\ Ls) Bs? Total yw; 
— = =n = Ar ~ Bs3 = 
kmax A B | A B . ” ; 
Ty 0.303 0.097 / 0.082 0.026 0.090 0.475 0.213 
6u 0.269 0.087 | 0.118 0.038 0.147 0.534 0.272 
7 ps 0.253 | 0.080 | 0.166 0.052 0.219 0.538 0.351 
Table IIb. The same quantities as the above for the case of g°/4z=0.10. 
S Ls} Hs2 Total ps 
kmax No A B | A B A B 
Dy 0.219 | 0.124 0.074 0.042 0.113 0.406 0.279 
6u 0.154 | 0.113 0.085 0.062 0.183 0.422 0.358 
7 ph 0.099 | 0.105 0.082 0.086 0.273 0.454 0.464 


The empirical value is 4,=0.439. The theoretical values are also in agreement with 
the empirical one qualitatively. 

iii) Vector parts of the magnetic moment 

The vector parts of the moment are given in Tables III. 

The empirical value is /4,=2.35. The theoretical values are certainly less than the 
empirical one, though it is not a discouraging discrepancy. The disagreement is due to 
two facts. First the modification w,—>«,* reduces the value of ty, from Miyazawa’s 
value (by about 30% in the case of 9° /47=0.08 and &k,..=6). Second the main 
correction /4,) is destructive and not so small. As is mentioned in the previous section 


Table Ila. The vector part of the magnetic moment for the case of 9 /4z=0.08. 


- | | res 
| fos / Total pa, 
Moy My? i) vg ] 
kmax A } B | A B 
Sp 1.145 | —0.447| 0.272 | 0.073 | 0.055 | 0.018 | 0.156 «| ) 1.25 | 1.22 
6m 1.462 | —0.627| 0.262 | 0.115 | 0.079 | 0.025 | 0.209 | 1.50 1.45 
1.770 | 0.271 @he uL75 1.67 
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Table IIIb. The same quantities as the above for the case of g?/4z=0.10. 


; = Hei Hr2 Lvs Hv Lev6 
kmax A B A B 
= SS ee EE a in ——- 
Dp 1.376 | —0.559 0.294 0.099 0.050 0.028 0.194 1.45 1.43 
6u 1.749 — 0.784 0.283 0.156 0.056 0.041 0.261 LP Hey/Au 
7 ps 2.105 — 1,026 0.274 0.225 | 0.054 0.057 0.339 1.97 eye 


a 
there is a bit of ambiguity of order 1/M in calculating ,, but it may not alter this 
Situation essentially. 

In conclusion we summarize the results as follows: first, the relativistic corrections 
are very important for the nucleon current and the magnetic moment, especially for the 
vector part of the current and the scalar part of the moment. The results are roughly 
in agreement with experiment, except for the vector part of the moment, with suitable 
choice of the coupling constant (g°/47=0.08~0.10) and the cut-off momentum (Anax 
=S5u~7p"). The vector part is about a few ten per cent smaller than the empirical 
value, and this comes from the reduction of the main term and the destructive effect of 
the main correction. Second, the renormalization of the electric charge and the nucleon 
mass are both considered in this paper. The renormalization of the nucleon mass has 
no effects within the approximation made here. As to the charge renormalization it was 
shown that the renormalized charge is equal to the unrenormalized charge up to the 
order of 1/M, if the neutron charge vanishes. The vanishing of the neutron charge 
was realized only approximately. 

The authors wish to express their thanks to Prof. H. Tanaka and many other 


members of the Institute for valuable discussions. 
Appendix 
Here we shall derive the relations (25). To justify (25a), it is enough to see 
(I4yl)=C)| XP P2345) =A, (A-1) 


where @ is a pure numerical constant. We consider the transformation of the space 
coordinate from (i, j, k) to (j, i, &), and refer the operator corresponding to this 


transformation as S. Then, the manipulation below can easily follow : 


(yl = | EXP #)sedel =C(S7'S | (X19). 435-5] 


=< \ (x (2!) XP 18 (2) ax de! |= | (-(2) XV ,6(2))exds|) 


=(|Aj|)- (A:-2) 
F ' sigs eat pine N . 
Since the only Hermitian symmetric tensor in the charge space is <;7,+7,7,; or 0;;, this 


ensures the equation (A-1). 
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To justify (25b), it is enough to show the relation, 
C1Bly = CLARes | XP i8).4} [= Ala (oeay—7s) > (A) 
where |u) is the bare spinor. In general, it can be written as 


(\Re| XP s6)ade|)= Cul Bore, + 795) [u), (A-4) 


where $f and 7 are the pure constants. The simple calculation shows 
[Ref XP sP)adel)* = (ul (BF 2, 0¢4+7* 94) |) 
= (ul (R,| (XP ,9)sds+2i| XP .9)sz5ds) |u) 


=(u| (80,0;+70,,+ 2ia@d,,) |u). (A-5) 


(A-5) implies 
pipe (A-6) 


Next, we apply Wigner time reversal T,, to (\B,|), i.e., 
(5|Bys|s’) = (u(s) | (7 +7*) |e’) 
=(G|Tz'T,BuTe' Te |¥)=—(—$ [Bal —3) 
= —<a(—s')bG-+r*) ja(—9))- 


Putting s’=s in (A-7), we have 


(A-7) 


r=-r* (A-8) 


(A:6) and (A-8) ensure the relation (A-3). The other relations can also be verified 
in the similar manner. 
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The exchange magnetic moment operators for the two nucleon system are obtained in terms of 
photopion production matrix elements. This photopion production part is closely related with the 
pion-nucleon scattering part, thus enabling us to use x—WN scattering experimental values in the ex- 
pressions of the exchange magnetic moment operators. These formulae are applied to cases of triton 
magnetic moment and those of heavy nuclei (the Fermi gas model). In the former case, rather a 
good result is obtained, while in the latter it is found that exchange magnetic moments can explain 
only small parts of the deviations from the Schmidt lines. 


§1. Introduction 


Many attempts have so far been made both phenomenologically” and theoretically?” 
in an effort to derive possible forms of exchange magnetic moment operators which come 
from virtual pion currents between nucleons. Every theoretical approach employed for 
this purpose, however, has been based on the perturbation expansion method in the 
coupling constant. 

A different approach is here made to this problem. It involves an expansion only 
in the number of pions exchanged, and the pion self field of each nucleon is treated, in 
principle, exactly. In this paper we compute the exchange magnetic moment involving 
the interchange of one pion. 

The process of our interest as a whole can be regarded as photopion production at 
one of the nucleons followed by the exchange of that pion with another. The matrix 
element of the photopion production process has been studied by Chew and Low” based 
on the static pion theory. Their analysis indicates that the matrix element is composed 
of two parts, one arising from pion current and the other from nucleon current.* The 
exchange moment coming from the former part, the pion current term, is essentially equal 
to the result of the lowest order perturbation theory, which has already been calculated. 

In § 2 we obtain the contribution from the nucleon current using our method, and 
the correction to the triton magnetic moment is evaluated in § 3. § 4 is devoted to the 
calculation of the exchange moment in the case of heavy nuclei by the Fermi gas model. 


§2. Exchange magnetic moment due to nucleon current 


In this section we shall derive the magnetic moment operator arising from the 


* This division of current is different from the conventional one. See ref. 4). 
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nucleon current by means of the S-matrix calculation. Before that, however, it is neces- 
sary to summarize the result of pion current contribution to the exchange magnetic 
moment. Chew and Low have shown that the photopion matrix element that comes 
from the pion current is essentially equal to the usual renormalized Born approximation 
if small terms are neglected. Under this approximation the exchange magnetic moment 
from this part agrees completely with that of previous calculations using the perturbation 
method in the lowest order. 

Using the static approximation throughout and restricting the number of pions ex- 
changed to one, the exchange magnetic moment operators due to pion currents for the 


two nucleon system can be derived in the 


following way. The processes we consider 

here are the interaction of the photon with 

the charged pion currents existing between a ae «ea. See 
pair of nucleons and the interaction with the : 

so-called interaction currents. This term, t 


nucleon 


which is evidently connected with the former = ------ pion 


é = F 3 —~~~~— photon 
on the gauge invariance requirement, 1S 


written as :* ae 


—jins A= —ef (7,¢.—72,) (6A). 


The diagrams for these processes are shown in Fig. 1. Note that in the latter process all 
of the photon, pion and nucleon interact precisely at one point of space. S-matrix for 
the first of these diagrams is this, 


Smeal ae oT [2 x r®],3(0) 


x | dPkdq[ (6k) (oq) { (ki +9, P20 (kg) 


= (ky +90) F 0° (k—q)} ad ef (Gx) — kx) ] 


oe wo, 


which yields after a little calculation 


Te wee 4 
= rh 7 },0(0) {ox (GF) Y—a!x, (G7) Y 


+[6x 6], Y—[x® x6], (6) ¥ 
—[x® xo], (67) ¥—[x® xx (12)}. 


In the same manner, we get the contribution from the second diagram 
? 


Sa Spl Xt },4(0) 


X {[x® xa], (67) Y— [x xo], (o@P)Y}. 


We take S=c=y=1; y is the pion rest mass. 


| 
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In the above expressions x“ is the position vector of the i-th nucleon and x denotes 
the inter-nucleon distance with the magnitude x. Also f is the renormalized coupling 
constant and Y means the Yukawa potential. All the other notations will be understood 
without explanation. 

Combining these results and noting that the quantity, which multiplies the factor 
2710 (0)H where H is the magnetic field, makes the magnetic moment operator, we get 
after some more calculations the following expressions for the exchange magnetic moment 
arising from the pion current. 


MQ, =—-F [2x 2, 
87 


a 2 @) ao) 5 
x| | x3 (X-O St a) (a+ 1 J=[oxo%,| en? 
oe 


+ | @? +2) xx} 72) | (1) 


where 


é 


1 9 1 1 ji re 
(o® ; x) (o® . x) 


=i(a®..6). 
oe 


T2=3 


No consideration has been made in the above evaluation of the pion-nucleon interaction 
expected to occur at the different places along the nucleon line. 

Next, and this being the main subject in this section, we consider the effect result- 
ing from the nucleon current. The photon interaction with this current is composed of 
both scalar and vector parts with respect to the isotopic spin space, and it has been 
shown by Chew and Low that the former is negligibly smaller than the latter. Therefore, 
we do not consider the scalar part here. The vector part of the interaction is of the 
form i(/,—/4,)/2-7;(0H), and the comparison with the ‘pseudovector coupling term 
ft«(Gk) of the pion-nucleon system indicates the possibility of replacing photon by a 
neutral pion with charge index 3, momentum H and coupling constant i(/#,—/.,) /2. 

With this interpretation, we can treat the photopion production part in the seme 
way as the pion-nucleon scattering process. The scattering matrix for the P-wave pion 


from state (ik) to state (jq) by a fixed nucleon is given by” 
(jase) = 270 (hy —q,)| A(by) #155 R) (09) 
+ By) {7:25(6q) (oe) +2,7,(ok) (oq)} 
+ C(k,) =)7:(6q) (ok )| OS 


(ik) are the charge index and momentum of the incident pion and (jq) are those of 
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the scattered pion. Values of 4, B and C for zero-energy (k,=0), which are necessary 


here, are given by 


A(0) =C(0) = 


© dp 
187 \ Oo, Ox ( Pp) 


w 
B(o) = 2 | Pon(p) 
67 J} Wy 

where 63, is the total cross section of the 
pure 3/2, 3/2 state. In this expression, ©, 
and o,, are neglected. The numerical value 
of the above integral is 

es 

| 3 -dp=3.7. 

ees - 
The most striking feature of this scattering - 
matrix element is that it can be computed 
using the experimental value of scattering 
cross section. By virtue of this fact, every 
effect occurring at the scattering part can be ag 
included in our result. 3 


The S-matrix for our process is given by 


Sai Pe Pe |g Cialsl3, HD 


xX 270 (qo) fr (o® q) eigx® ‘ a? a ; 1 t 
(27) Jo — 


+terms with nucleon numbers interchanged, 


which describes the processes that the external photon (considered as the equivalent pion), 
interacting with one of the nucleons, produces the virtual pion, which is subsequently 
absorbed by another nucleon. The Feynman diagram of this process is shown in Fig. 2. 
After a little calculation of this S-matrix, we obtain the final result of the exchange 
magnetic moment due to the nucleon current. 

They are: 


i=1 


where 


a 2 
M,= (<9 +2) {320 40") 


a (6 46%) |) 


ee 
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M,= (739 +7) (6% +6) f(x). 


pal ier (2) x x-G%)—G® $e 
M,= (=P — 2) {3 ( 5 ) race oD) 


x 


M,= (2 —2) (6 —o®) f,(x) 


9 


Mize i 79 x 7], 3 x(x. 5X a) —[{o% x off, (x) 


M, -- [7 4 sO. [o™ x o®]f, (x) 


and 
ye Ae ee 
fal) =fal®) =fal3) = ——— Pf 40) 
fo@) =~ He 4(0)( 4 4 Je. 


Thus, the total exchange moment is given by the sum 


es ) N 
M as MS ° Mee : 


§ 3. Application to the triton magnetic moment 


As the first application of the results obtained in the last section, we consider the 
triton magnetic moment. As is well known, its experimental value is larger than that 
in the Schmidt model, and the admixture of states in the triton ground state cannot be 
responsible for this deviation since it has the wrong sign for this effect. So we want 
to see to what extent the pion effect could explain this disagreement. Assuming that 
only one pion is exchanged between a pair of nucleons and that none between the other 


pairs at the same time, the magnetic moment operators for the three nucleon system can 
be given by 
M(123) =M(12) +M(23) +M (13). 
Now the expectation value of M“™ for the three nucleon case with respect to its 


ground state gives the numerical value ~+0.3 n.m. and this makes a part of corrections 


due to the exchange pion current.” 
To obtain a further correction, we must evaluate the expectation values of the 
operators listed in eq. (2) (nucleon current contribution). For the triton wave function 


we take 

a(1)n(1)  @(2)n(2)  a@(3)n(3) | 
a(1)p(1) a@(2)p(2) a (3) p(3) ¢(123) 
B(1)n(1) B(2)n(2) 8(3)n(3) 


= 
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where @, 9 and p, n are spin and isotopic spin functions respectively. ¢(123) is the 
part of space variables, totally symmetric with respect to inter-nucleon distances X,9, Xs 


and x,,, and assumed to be of the following type :” 
(123) =exp[ —A (Xp +%q+%1) | 


where / is a parameter with the numerical value of 0.5033. 

Calculations are easy and we obtain the result that the contribution of M.{Q, terms 
to the triton magnetic moment is about —0.062 n.m. This value, when added to the 
above mentioned one which is the same as Villar’s result, forms evidently the total 


modification to the triton magnetic moment caused by the exchange currents. 


§4. Exchange magnetic moments of heavy auclei | 


In this section, we calculate exchange magnetic moments of heavy nuclei. For the 
nuclear wave function, we choose the Fermi-gas model wave function %, since it is the 
most convenient one for calculations and since our main purpose of this calculation is to 


see the degree of the deviation of nuclear magnetic moments from the Schmidt lines. 
The function is 


1¢,(1), aaa at > @x (1) 
pop 
1¢,(N), ee » ey(N) 
where 
tay \nachen Wee tthonans p() Feit (1) 
g;(1) Re ing @(1) 151) or (eee 3(1) {FO 
etc. for 1S i<N—1 
and 
R(r) Y/ (0d) (1) (POD saya 
() Yi) a(a) Rey fla 
gy (1) = pears <= | @) 
'-1 (Hd) @ x 2 nay {pq Bis 
ROY) |) apy tt 8) ¢) — 9/7 ¥08) 8a) } {20% jig. 


In this model, N—1 is assumed to be even and both spin and isotopic spin are saturated. 
R,(r) in the above expression satisfies the following relation : | RiP) jkr) rdr=cd (k—k,), 
where ky is a Fermi momentum, V corresponds to y= 1 a unit “=1, where r, is 
(nuclear radius) X (nucleon number) ~"/°, 

Since the wave function used here is merely an approximate one, we cannot hope to 
obtain the exact value, but at least it will give a correct order of magnitude since the 


Fermi gas model wave function proves to be a fairly good approximation, when nuclei 
are rather heavy and a reasonable cut-off length could be chosen. 


In table we list the final results for proton (in nuclear magnetons) thus obtained. 
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where k, is the cut-off momentum (in unit “=1). Signs must be changed in these 


values in the case of neutrons. 


§ 4. Concluding remarks 


Our new approach to the problem of the exchange current effect on the nuclear 
magnetic moments has been shown to be very powerful in that the experimentally 
observed results have been adopted in a framework of our theoretical description. No 
detailed knowledge of the interaction mechanism is needed, only if the initial and final 
states could be exactly specified by observation. 

The main reason for that we have chosen H® as one of the examples is that it was 
the first nucleus to indicate the possible effects of the exchange currents in it. The 
deuteron magnetic moment anomaly, on the other hand, can be explained without taking 
into account any of the exchange effects, that is by means of the admixture of the D 
state. But, this is not the case with triton moment. No reasonable admixture is enough 
to account for the anomaly of triton magnetic moments. Because of this fact, we at- 
tempted to explain this situation in terms of the exchange effects. 

It turned out, however, that the whole exchange effects gave us the correction of 
about ~+0.24 n.m. to be added to the proton moment, and this amount of correction 
evidently exceeds the observed value by about 0.054 n.m. This deviation might possibly 
be brought to zero by the admixture of reasonable states. 

As to the second example taken up in our calculations, a brief remark must be 
added. Although the shell model has succeeded in explaining various properties of atomic 
nuclei, the quantitative discussions of them using this model cannot yet be made with 
satisfaction. For instance, the deviation of the magnetic moment from the Schmidt lines 
still remains large. For this deviation many sorts of explanations have been given, some 


of them being the modifications of the shell model itself. However, even in such an 
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attempt it is important to notice how large the exchange magnetic moments are. These 
exchange moments of heavy nuclei were calculated by Miyazawa,” and indicated the 
possibility of obtaining the large value enough to account for the deviation from the 
Schmidt lines. But, Miyazawa’s calculations were performed using a very large coupling 
constant (>0.25) and neglecting the finite extension of the nucleons, and this fact 
would be the cause of the large values of his results. Our results, on the contrary, 
turned out to be small, so we must seek for some other way of the explanation thereof. 

To conclude our discussions, the authors are very grateful for the guidance by and 


many useful discussions with Dr. Miyazawa. 
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In order to explore reliable consequences of the pion theory of nuclear forces at high energies, 
the p—p scattering at 90 Mev is discussed in detail. No evidence has been found to show that the 
static pion-theoretical potential is appreciably modified up to 100 Mev. This conclusion is in conflict 
with the prediction of very strong spin-orbit coupling potentials recently made by Signell and Marshak 
and by Gammel and Thaler. Especially, it is discussed that the former spinorbit term was introduced 
in reliance on Gartenhaus’ potential at small inter-nucleon distances. These spin-orbit coupling 
potentials are unreasonably strong from the pion-theoretical view-point. Importance of the °P,—*F, 
coupling parameter in determining angular distributions is pointed out. 


$1. Tntroduction 


In our previous works on the verification of the pion theory of nuclear forces,” it 
was shown that all experimental data at low energies (below 20 Mev) are quantitatively 
explained without any exception and that the explanation is made mainly in terms of the 
outer part (outside the pion Compton wave length) of the static pion-theoretical potential. 
The effective —N coupling constant was thereby determined as 9,’/47=0.080 +0.010. 

One of the essential points of our method of approach used there is dividing the 
inter-nucleon distance into three regions and distinguishing the reliability of the potential 
or the nuclear interaction in each region. It was also discussed that the validity of this 
method is necessarily limited to regions of not so high energy (below about 100 Mev), 
since: (i) Interactions at small distances, which cannot be defined by the present-day 
quantum theory, become important at the high energy. Especially angular distributions 
are seriously dependent on the inner interactions mainly through the P-wave phase shifts. 
(ii) Dynamical corrections to the outer part of the potential are expected to be appreciable 
and modify it. Properties of the corrections are directly related to properties of the pion- 
cloud and possibly to the nature of the nucleon core as the source of the cloud. 

It may be said that we have had no serious attempts to explore reliable consequences 
of the pion theory of nuclear forces in the high energy region.** Works on nuclear 
forces around 100 Mev can be roughly divided into two groups. The first group” 
consists of works of phase shift analysis.” Unfortunately, these works seem to be performed 


under purely phenomenological and mathematical programs, utterly ignoring the expectation 


* Some of preliminary results of this paper have already been published.) 


** For references, see ref. 1, p. 71. 
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that some of characteristic features of the low energy pion theory should still remain 
unchanged to some extent at high energies. Moreover, we cannot expect the phase shift 
anlysis to give a unique solution easily. The second group includes studies of 
mechanical application of the static pion-theoretical potential, even to scattering at 300 
Mev. Since the reliability of the potential as discussed in detail by our group” is 
scarcely taken into account in these works, we can get rather few informations on the 
pion theory of nuclear forces from them. 

In the present paper, consequences of the static pion-theoretical potential established 
at the low energies are compared in detail with the p—p scattering experiments at 90 Mev. 
As will be discussed in § 2, D- and higher-angular-momentum waves (with the impact 
parameter larger than the pion Compton wave length) are given their phase shifts by 
the outer part of the pion-theoretical potential. Low-angular-momentum waves are, on 
the contrary, treated more or less phenomenologically as is usually done in the phase 
shift analysis. In the same way, the n—p scattering is analysed by Watari in the 
following paper.” In our analyses in § 3, we shall find no evidence to show that the 
static pion-theoretical potential is modified appreciably. As will be discussed in § 4, any 
velocity-dependent interactions would be, if any, still weak around 100 Mev, and no 
drastic changes may be expected in the structure of the pion-cloud at such energies. 

Analyses at higher energies (~150 Mev) have been done by Tamagaki.” It is to 
be noted that although his analyses are necessarily far more phenomenological than ours, 
he obtained the same conclusion on the velocity-dependent interactions. In §5, remarks 
will be made on spin-orbit coupling potentials recently introduced by Signell and Marshak” 
-and by Gammel and Thaler.’” 


§ 2. Phase shifts and the applicability of the pion- 
theoretical potential 


According to the division of the inter-nucleon distances into three regions,’ we 
classify phase shifts also into three categories : 

I) Phase shifts with L2 are treated pion-theoretically, where Lb is the angular 
momentum. Since the D-wave impact parameter 6,=1.7*, they are almost completely 
determined by the outer part of the one-pion-exchange potential in the region I (x > 1.5). 
They are shown in Table la. 

II) P-wave phase shifts with 6,=1.0 are under the influence of both the one- and 
the two-pion-exchange potentials. They are reliable as far as their qualitative features are 
concerned. Not only *P-wave phase shifts °d}, ‘0! and °0%, but also €, (the *P,—‘F, 
coupling) and *dT (the eigen °F,-wave phase shift) fall into this category. 

The phase shifts of this category would be most sensitive to various dynamical 
corrections to the potentials, since, if such corrections are important, they would be in 


the region II (1.5 >x> 0.7) where the pion cloud has not the classical nature as in 


Unit of length of the inter-nucleon distance is the pion Compton wave length (1.415 x 10-1 cm) 
Throughout this paper notations and terminologies are the same as those in reference 1. 
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Table 1. Phase shifts and coupling parameter with L=1 at 90 Mev. They are 
calculated by assuming the Taketani-Machida-Onuma potential with g,?/4z=0.08.5a).6») 
This TMO potential is considered here as the standard since in the region I it has 
the correct asymptotic form, and in the region II it has all features that are common 
with various potential having the plausible pion-theoretical basis. 


a) L=2. Pion-theoretically quite reliable. 


state | phase shift value 
1D, 19 | 0.047 
| 
3F, 398 | — 0.026 
SF, | 398 | 0.003; 


b) L=1. Only qualitatively reliable. The inner interaction is cut of by a hard 
core of the radius x,=0.33. 


state | parameter value 
3Py tye 0.377 
3p, ag —0.187 
3P,+3F, 35% 0.065 
sgt —0.035 
€, (coupling parameter) —0.704* 


* The sign for €: was missed in the calculation of ref. 6b. 


the region I and is relatively complicated. However, their features are strongly dependent 
on the cut-off momentum of the virtual pions and we cannot derive the corections on 
a sound basis. 

At 100 Mev the phase shifts of this category are also appreciably dependent upon 
the inner interactions in the region III (x 0.7) if these interactions were to be strongly 
attractive. However it is possible and reasonable, as the first approximation, to make 
the inner interactions of secondary importance for the phase shifts, up to about the excitation 
energy of the nucleon isobar, by means of some cut-off procedures in the coordinate 
space. As will be discussed later, it seems to us that undue reliance on the region III 
interactions of Gartenhaus’ potential” made Signell and Marshak”) introduce an unreasonably 
large spin-orbit coupling potential. 

Aside from non-static potentials which are sufficiently small to be neglected in the 
region I, the static pion-theoretical potential has the following characteristic features in 
the triplet odd state as can be seen from Fig. 1: 

The one-pion-exchange potential 
+ (negative tensor and attractive central) two-pion-exchange potential. 
It may deserve special attention that this. potential predicts the *P-wave phase shifts such 


as 


tj Ono. 20, 0 0, and , *0y~2|"0.). (1) 


These relations come from the fact that the one-pion-exchange tensor potential is much 
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stronger than the central 
potential as well as the two- 
pion-exchange potential, for 
in the Born approximation 


we can expect 0.05 


sgjoc — V,—4V 2), 
aio (7, e)y  \2) 


and 
39200 — (V.—0.4V 7) 
+ (contribution coming 
from the *F,-coupling). 
We shall call *P-wave phase 
shifts satisfying eq. (1) 


those of “ pion type ie 
One example of the 


—1.0 1.of —0.1 > 


pion type phase shifts and 
the mixing ratio are given 
in Table 1b. The potential 


is cut off by a hard core Fig. 1. Pion-theoretical potential in the triplet odd state. 


OPEP : The one-pion-exchange potential with g,*/4z =0.08. 


with a radius x~0.3. TMO,*) FST KMO"™®) : one-plus-two-pion exchange 
III) Singlet S-wave phase potentials. 

ale 5 ae es For derivation methods refer reference 7, and its Fig. 

shift ‘0, is dependent upon TV anna: 


the unknown interactions in 
the region III(x\0.7). 
Therefore we should treat it phenomenologically and determine in such a way that the 


total cross section (or do/d{2) at 90 Mev may be reproduced. 


§ 3. Analysis of p—p scattering at 90 Mev 


When we assume the “P-wave phase shifts of the pion type, the angular distribution 
do/dQ and the polarization cross section P(())do/d®Q are effectively dependent only upon 
one parameter, the *P,—‘F, coupling €,. We would discuss the €,-dependence of do /dQ 
and of P(#)do/dQ and determine €, rather phenomenologically. Reasonings for this pre- 
scription are: (i) The strong positive tensor force predicts eq. (1). (ii) Of the L=1 
scattering parameters, only €, is proportional to &* at low energies, while the others at 
least to k*°, which seems to mean that the energy dependence of €, is complicated and 
may be sensitive to the inner interactions. 

First, we will discuss the polarization cross section, which is shown in Fig. 2, since 
P(@) is due to the tensor force and is very sensitive to €,. Main terms of P(@)do/dQ 


a 
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can be expressed as follows :* 


P(0) do /dO~ (6/k«*) f(€,) {[01|20]+[27|01]+ (3/2) [11 2a] + (3,/2) [27|11]} 


where 
(However, for convenience sake, 
we use °0) and °0! instead of 307 


and *0f as done in ref. 1.) 


and 


() =cos 2€ + - A sin 2€. 
2V 6 


It is very convenient that 
P(@)do/d2 is separated into two 
parts: the €,-dependent part f(€,) 
as shown in Fig. 3 and the 0- 
dependent part in the parenthesis. 
When the *P wave phase shifts are 
of the pion-type as eq. (1), all 
terms of the 0-part in eq. (3) are 
positive. This coincidence of the 


sign, which is due to the strong 


X sin @ cos /, (3) 


[Je J’ 0 |= sin 0% sin°0e, sin (908—80"), o=a, f, 7, 


(mb/ster) 


0.5 


P(d)do/dQ 


0 (degree) 


Fig. 2. Polarization cross section at 90 Mev with the 
phase shifts in Table 1 and eq. (9). & 
the experimental data, see ref. 2. 


is varied. For 


positive tensor potential, is quite favourable to the large experimental P(#/)do/dQ2 and 


gives 


P(@) do /d2~f (€,) sin cos# X 1.03 mb. (4) 


Including the higher wave contributions as given in Table 1b, we have 


P(A) do /d2=f(€,) (0.98 +0.15 cos’#) sin cos@ mb. (4’) 


By comparing eq. (4’) with experimental data at 95 Mev” 


P(@)do/d2= {(1.00 + 0.11) + (0.16+0.13) cos’ 4} sin# cos@ mb, (5) 


we have f(€,) ~1. 
we can estimate by Fig. 3 


When we take the negative sign of €, corresponding to Table 1b, 


€, (pol) = —0.15. (6) 


It may be worth while to point out that P(#))do/d is very sensitive to €,. 
Actually, for €,=€,(TMO) =—0.704 in Table 1b, P((/)do/d& is vanishingly small 
owing to the accidental coinsidence with one of the roots of f(€,). 

We shall then discuss the angular distribution do/d2 : 


do /dQ= (1/k x?) S1B,P, (cos 4), (7) 
L 


* For complete formula, see reference 5). 
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where B,= 2 JV) (sind*)* is ‘ 
independent ee €,, while the other i 
B,’s are not. Main terms of B, 

are shown in Table 2.* 

As shown in Fig. 4, the 
angular distribution becomes un- 
isotropic as €, increases. biSeise es ite) 
due to the following situations. 

As can be seen in Table 2, all 0.5 
main terms of B,, the coefficient | 
of WEs(cos @)=4.G cos’ —T).,,. ate 
non-positive at €,=€,(TMO), 


e,(1MO) 
while they become non-negative 


without any exception as €, in- 
creases. Therefore, as long as the 
*P-wave phase shifts of the pion 0 
type are assumed, the angular 
distribution becomes more and more 
unisotropic as €, increases. We 
can estimate therefore / 


€, (a. d.) S—0.35. (8) Fig. 3. f(€) defined by eg. (3). 


The small inconsistency between €,(pol) = —0.15 and €,(a. d.) <—0.35 is, however, 
not serious at all. The 0-dependent part of P(#)do/dQ is sensitive to the magnitudes 
of the *P-wave phase shifts (proportional to 0° for small 0). Therefore if 4 increase a 
little, do/d2 and P(@)do/d@ can be reproduced with one value of €, consistently. An 


example is as follows. Let J denote a variation of the contribution from the tensor 


Table 2. Main triplet terms of By and their €:-dependence. Notice that as €» increases, none of them decreases. 
By= > , Prputpr (Es (Jello, 


Jp. Mp 
where 


(Jo|J’0’) = sin °9?, sin 88?/, cos (85°, —89°)). 
SS 


Jo Jo" (Jo\J’ 0’) | brpztor (€2) “ae 4 Same | oa 
11 ll 0.035 3/2 a positive constant 
01 2a 0.023 | 4 [cos €)+V 8% sin €9] 2 ~OC + + 
01 2r —0.012 — 4[—sin €:+V 8 cos €5] 2 +>+>~0 
11 2a —0.012 9 [cos €) —V 34 sin €9] 2 | +>+>~0 
11 2% 0.007 9 [sin €.+ V 2% cos €o] 2 ~0<+>~+4+ 


z ante We ¢ 
The complete expression is given in reference >) 
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potential to the *P-wave phase shifts. _ (mb/ster) 
In eq. (2) let us assume that d= 
0.01 which corresponds to about 10% 
increase of the tensor contribution of 
the TMO potential. Such uncertainty 
of 10% may come from the uncertainty 
of the coupling constant or from the 
two-pion-exchange potential whose pre- 
cise shape is not known as shown in 
Fig. 1. According to eq. (2), 4 
adds 44 to *0}, —24 to °0} and 0.44 


to *0%, and we have 


e 60 90 
a1 __ 551 ae 
Oo= 0.42, “03=— 0.21 eee) 
and “je ==0.07- (9) Fig. 4. Scattering cross section at 90 Mev with the 
phase shifts in Table 1. € is varied. 146) is assumed 
This example of the *P-wave phase as 0.50. For the experimental data, see ref. 2. 


shifts gives precise fits to do/d@ and 
P(@)do/d2 as shown in Figs. 2 and (mb/ster) 
5, and for €,~ —0.15. 


§ 4. Concluding remarks 


i) Conclusions 

In § 3, it was shown that p—p 
scatterings at 90 Mev, and consequent- 
ly those below ‘100 Mev,* can be 
reproduced by the outer part of the 


dg /dQ 


static pion-theoretical potential, if the 


inner interactions in the triplet odd 


state are made non-essential by a hard- 0 o be 8 
core cut-off. Watari also obtained the Mees 
same conclusion in his analysis of ae Fig. 5. Scattering cross section at 90 Mev with the 
scatterings at 90 Mev,” that the shifts in eq. (9). €: is varied. 109 is assumed as 0.40. 
higher wave phase shifts, which are For the experimental data, see ref. 2). 


pion-theoretically reliable and classified 

under the categories I and II, are not contradictory to the experimental data. Namely, 
we find no evidences to show that the static pion-theoretical potential is appreciably 
modified below 100 Mev. No drastic change may be expected in the structure of the 


pion cloud. Our conclusion is conflict with predictions of strong spin-orbit coupling 


* A complete fit to p—p angular distribution at 18.2 Mev has already been obtained.") 
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potentials recently made by Gammel and Thaler, and by Signell and Marshak.” This 

point will be discussed below’. 

ii) Remarks on Gammel and Thaler’s spin-orbit coupling potentia 
Gammel and Thaler discussed about an evidence upon the spin-orbit coupling force 

making use of Stapp’s results of the phase shift analysis of 310 Mev p—p scattering. 

They proposed a potential which has the following qualitative features : 


pPD* 


(a) a very weak central potential. 

(b) a positive tensor potential with a long range (range~0.9). 

(c) a negative spin-orbit potential with a short range (range~ 0.2). 
In the features (a) and (b) Gammel and Thaler’s potential is similar to the pion- 
theoretical potential. Hwever, above as low as a few tens Mev, Gammel and Thaler’s 
phase shifts cease to be of the pion type. For example, as shown in Table 3, they do 
not satisfy the inequality (1) but 


39¢ ~ 391 > 0 > 401. 
This deviation from the pion type is due to their strong spin-orbit potential : 
Vo Vie ae, 
V,=—7122.5 MeV, 
a =0119: 


Although the range is small, its depth is exceedingly large. Actually, the spin-orbit 
potential is much larger than the pion-theoretical tensor potential inside x~0.9. There- 
fore, Gammel and Thaler’s potential gives fits to the experimental data in quite a different 
way from ours. In view of the success of the low energy pion theory, it would not be 
reasonable to suppose that the pion-theoretical aspects of nuclear forces may be lost even 
at as low as a few tens Mev. 

iii) Remarks on Signell and Marshak’s spin-orbit coupling potential” 

Recently Signell and Marshak also concluded the necessity of a spin-orbit coupling 
potential, which is added to Gartenhaus’ potential,”) by analysing the p—p scattering 
data at the energy region around 150 Mev. We think that the disagreement between 
their results and ours comes from different ways of thinking about the inner interactions 
and, consequently, different methods of the cut-off in the triplet odd state as will be 
explained below. 

As discussed in detail by Japanese physicists, the inner interactions in the region III 
is beyond the applicability of the present-day pion theory and should be treated phenomeno- 
logically." Availing ourselves of this point of view Gartenhaus’ potential was commented 
already in 1956.” However, Signell and Marshak state that “there is no free 
parameters in Gartenhaus’ potential,” and do not give proper phenomenological prescription 
to it. It is interesting to compare phase shifts by Gartenhaus’ potential and those by 


ft See also note added in proof. 
* The author wishes to thank Dr. Gammel and Dr. Thaler for sending him their preprint. 
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Table 3. Phase shifts and coupling parameter with L=1 by various potentials. 


39h spt 392 3a € 
Gammel-Thaler,!2) 80 Mev 0.188 —0.223 0.168 0.002; 0.326 ; 
100 Mev 0.153 0.251, 0.196 0.008, 0.291 
Gartenhaus®):®) 90 Mev 0.93 —0.12 0.21 —0.06 0.23 
Signell-Marshak®) 100 Mev 0.276 —0.244 0.171 —0.026 —0.258 
Ours (Table 1b) 90 Mev 0.377 —0.187 0.065 —0.035 ~| ae 
(eq. (9)) 90 Mev 0.42 aOi21 0.07 —0.04 0.15 


Signell and Marshak’s one tabulated in Table 3. Large positive value °0) and small 
magnitude of *d; of the former are due to the very deep inner central potential (—1.5 
Bey at x=0), and cannot explain the p—p scattering data. From Table 3 we can im- 
mediately see that the effects of this inner potential, which is absolutely meaningless 
from the pion-theoretical point of view, are largely cancelled by the spin-orbit term 
and that the resultant phase shifts are not similar to Gammel and Thaler’s but rather 
to ours. It should be pointed out that, contrary to Signell and Marshak’s discussions, 
Gammel and Thaler’s spin-orbit potential reproduces the experimental data in quite a different 
way from Signell and Marshak’s, as shown in Table 3, especially, in different signs of €,. 

The cut-off that Signell and Marshak adopted in their calculation is the zero cut- 
off at the nucleon Compton wave length. If a phenomenological cut-off by a hard core 
were assumed throughout, the results would be much changed, since the outer part of 
Gartenhaus’ potential is quite the same with ours. Therefore, we could conclude that 
it is the undue reliance on the inner interactions that necessitates Signell and Marshak 
to introduce such a large spin-orbit potential. 

It is interesting that the same conclusion is obtained by Tamagaki in his analyses 
of the 150 Mev data.” 

iv) Remarks on the phase shift analysis 

In most of the works of the phase shift analysis,” the “P,—°*F, coupling is completely 
neglected for the sake of simplicity, namely it is assumed that ¢,=0. 

Our analysis in § 3, however, shows that both the scattering cross section and the 
polarization are very sensitive to €,. One could hardly claim that results of the phase 
shift analysis are correct, if €,=0 is assumed identically. Moreover, the presence of the 
strong tensor force of the pion-theoretical potential does not justify neglecting the 
3P,—F, coupling at all. Simplicity of analysis should not be preferred to the reality of 
the pion theory. Rather, correlations between 39%, 803 and €, should be investigated under 
the assumption of such a strong positive tensor potential as the pion-theoretical one. 


The author would like to express his sincere thanks to Dr. R. Tamagaki and Dr. W. 


Watari, whose stimulating discussions have helped him develop this work. 


eb S. Otsuki 
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Noce added in proof. Recently Sato calculated the spin-orbit coupling potential making use of an ex- 
tended method of ref. 10b, ie., taking into account the features of the pion clouds somewhat phenomeno- 
logically*. His result is almost the same with the previous one**, which was calculated by the Fukuda- 
Sawada-Taketani formalism!“ treating the normalization condition correctly for the Tamm-Dancoff ampli- 
tude. The latter potential is: 


V is=— (9°/4z)? (2p°/M) [3+2(¢-e)] 
X (1+2/x+2/x2+ 1/x5) exp (—2x) /23- 


These pion-theoretical spin-orbit coupling potentials are by far weaker than the phenomenological ones, 
especially in the outer region, as shown below : 
Magnitude of V7,s(aside the negative sign, in unit of 0.0lze?, » the pion mass) 


x=1 x=2 
pion-theoretical one 0.77 0.0058 
Gammel-Thaler 5.4 0.019 
Signell-Marshak 57 0.30 


Thes fact supports our conclusion that any velocity-dependent interactions are not strong at least 
around 100 Mey. 


* A private communication from S. Sato. 
** See ref. 7, p. 126, 
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The n-p scattering at 90 MeV is discussed in detail from the standpoint of reference 1). Around 
this energy phase shifts with L>2 are almost completely determined by the one-pion-exchange poten- 
tial. Phase shifts with L=1 should be qualitatively determined by the pion theory. Phase shifts 
with L=O must be determined phenomenologically. In this paper, using P-wave phase shifts obtained 
by Otsuki from the p-p scattering data, we determine scattering parameters in the (J=1)-states from 
the n-p scattering analysis. In section 2, n-p polarizations are analyzed and some sets of phase shifts 
with J=1 which are consistent with the experiments are determined. Using thus determined phase 
shifts, angular distributions are calculated and compared with experimental data in section 3. Thus 
we find that the static pion-theoretical potential can reproduce n-p experimental data at 90 MeV 


without any serious modifications. 


§ 1. Tnatreduction and summary 


In 1956, we discussed low energy properties of nuclear forces in detail from the 
stand-point of the pion theory and found that all of them are quantitatively explained 
in terms of the static pion-theoretical potential with the coupling constant 9,°/47=0.080 
+0.010.2** This success is a motivation of the present work, where it is discussed that 
the pion-theoretical potential established at low energies can reproduce n-p scattering at 
90 MeV. At 90 MeV, we have two kinds of useful data, namely, single scattering 
cross section do/d over all angles” and polarization cross section P(())do/d2”. We 
already calculated the former in 1953” and the latter in 1956”, assuming a static pion 
theoretical potential derived by the perturbation approximation up to the fourth order 
(the TMO potential). The agreement between the experiment and the theory was fairly 
satisfactory. However, the success of the low energy pion theory, seems to require 
more detailed investigations for high energy phenomena; especially we must find out 
which features of the potential bring about the agreement and what modifications of 
them make the agreement better. 

Around 90 MeV, corrections to the static pion-theoretical potential, e.g. the rela- 
tivistic effects, the recoil effects, etc., become appreciable. Unfortunately we have few 


knowledge about these corrections. But they are expected to be still small in the region 


* The preliminary report of this work has been published in this journal®). 
** Notations and terminologies used in this paper are the same with those in ref. 1). 
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I. Even in the region II they may not be so large as to mask main features of the 
static potentials in this region. Therefore it is desirable to investigate properties of these 
corrections phenomenologically by analyzing experimental data with the static pion-theoreti- 
cal potentials established by low energy analyses. 

On the other hand, at the energy of 90 MeV, the impact parameter of P-wave is 
b,~1.0 and so the angular distributions are not seriously affected by the two-nucleon 
interactions at small distances”. The analysis from the standpoint of the pion theory 
developed in reference 1) is still meaningful in the problem of nuclear forces. 

The p-p scattering has been analyzed by Otsuki”, and his results are referred to in 
our analysis. 

We shall divide phase shifts into three classes as was done in the p-p analysis. 

(1) Phase shifts with L>2. Impact parameters of the wave of L>2 are larger 
than 1.7. They are almost completely determined by the one-pion-exchange potential in 
the region I. Their values are tabulated in Table I. These values are reliable within 
the uncertainty of several per cent due to the ambiguities of the coupling coustant and 
the inner region nuclear forces. 

(II) Phase shifts with L=1. P-wave impact parameter is 6,~1.0, and their phase 
shifts are mainly determined by the potential in the regions II and I. Therefore their 
qualitative features can be regarded as pion-theoretical. 

(III) Phase shifts with L=0. S-wave is much influenced by the interactions in the 
region III. So we can say nothing about S-wave phase shifts from the present-day pion 
theory. They must be determined phenomenologically as will be done below. Their 
magnitude and energy dependence may reflect some features of the interactions in region 
III. 

We analyze the n-p polarization in section 2, and determine phase shifts with L=1 
and 0 which are consistent with the experiment. Using thus determined phase shifts, 
n-p angular distribution is calculated and compared with the experimental data in section 3. 

Our conclusions are as follows : 

(i) According to the method of approach proposed in reference 1), the static 
pion-theoretical potential established at low energies can reproduce n-p experimental data 
at 90 MeV without any serious modification. Considering Otsuki’s result of the p-p 
analysis, we could say that the pion theory of nuclear forces explains all experimental 
data up to 100 MeV correctly. 

(ii) In the triplet odd state, however, the tensor potential seems to be more or 


less stronger than that of the TMO potential. This result is consistent with a disposi- 
tion of the FST potential. ; 


Table I. L=2 phase shifts in radians determined by the pion theoretical potential 


Triplet Singlet 
) nn °F, Ds Po ae 
0.29 0.00 — 0.03 0.00 0.05 — 0.04 0.00 
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(iii) We find some evidences to show that in both the triplet and singlet even 


states, hard-core-like repulsive interactions are present at small distances. 


§ 2. n-p polarization 


The polarization cross section P(#)do/dQ is analyzed firstly by the following reasons. 
(i) Phase shifts of the singlet states do not contribute to this phenomena. 
(ii) P(@)do/dQ2 includes phase shifts in a form |04|0,|==sin 04-sin 0,-sin(0,—0,). 
Therefore higher wave phase shifts are safely neglected since [04|0,|~0° for |04|~|0,| 
~O<1. 


The n-p polarization cross section is expressed in terms of P,(cos @) as follows : 


PO) do fd k™ sin 7 ag P; (cos). 
LE 


Experimental results by Stafford et al.” are as follows and are shown in Fig. 1. 


P(@) do /d2=sin 6 (1.17 0.11) P+ (2.34 + 0.23) P,+ (0.21 40.26) P,| mb/sr. 
( for 95MeV) 
avsin O{ (1.12 40.09) P,+ (2.28 + 0.21) P,| mb/sr. 
When the phase shifts by the TMO potential shown in Table II-a are assumed”* we 


have the following result : 
P(@)do/d2=sin 6/0.61P,+2.14P, + 0.03P,—0.02P;| mb/sr. 
Signs of the main terms, a, and a,, coincide with the experimental ones. However, 


since their magnitudes are a little smaller, we would discuss this point more in detail 


below. 
Main parts of a, and a, which we denote as (a,), and (a,),, are as follows :** 


Main parts of a, (due to interferences between different parity states) 


*P,— (5,+°D,): 9 (€,) {L07|1¢]—Lo7|17]}, 
eee Mapas): (3/4)9(€) {[18|1@]—[18|17]}, | 
i (Can 
CGP, +°F,) — (5,4+°D,): Fé» €)[1a\2¢]+F(«, +%/2)[1¢|27] 
+F(€,+7/2, €,)[17|2a]+F (€,+7/2, et 


where 
[ Je|J/o’ |=sin 39° sin °0,,°' sin (70, —°0,,°’), 


9 (€,) = (3/2)/2) sin(2€, +), tan i ivy 


* The first term of the last equation in page 1204 of reference 5) is not correct. The revised equa- 
tion is given here. The equation in page 1206 should also read as P(do/d2) Sein Xi-5F (cosh) mbies, 
Consequently, curves of Fig. 1 of reference 5) are not correct. The first two sentences (from “However” to 
“+easons.”) of the second paragraph in page 1205 must be deleted. 2 

** Complete expressions of coefficients a7, are given in a forthcoming paper by R. Tamagaki”. 
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dg 


oe (ab (mb/ster) 


90° 


@ (c. m. system) 


Fig. 2 


Fig. 1 n-p polarization cross section. Experimental 
data are those of ref. 3). The curve is calculated 
using phase shifts given in page 9. 


and 
F(€,, €:) =cos® €,+4 cos? €, cos? €,+ (5/8,2) sin 2€, (cos* €,— (2/3) sin*€,) — (5/6) sin*€, 
— (7/12) sin? €, sin® €,+ (9/8)/6)sin® €, sin 2 €,. 


V2 9(€) is plotted in Fig. 2. 


Main parts of a, (due to interfences between same parity states) 
CS, +°D,) — CS, +*D,): — (9/4) 9(€,) [18 \17), 


| (a) 
CS,-+-*D,) =<), : (5/4)97(€,) {[23|1a@]—| 28 \17]}. 


As was discussed in detail in the p-p analysis” the triplet odd state phase shifts 
should rather be restricted in our analysis. Since the pion-theoretical potential has a 
strong positive tensor part and a weak central part, the phase shifts should be as follows : 


3A 1s 8120 8) es 3) 1 18) 1} 
0, > 0, “0, <0, “d,*=0, and *d,'~2]*0,'|. 


Their values, which are calculated with the TMO potential (set a) and which are 
modified in the p-p analysis (set 6), are tabulated in Table II. Table II contains also 
some sets of phase shifts which are taken in the present analysis. 

We should like to determine phenomenologically the °S,+-°D, state scattering para- 
meters, °0,", °0,7 and €,. For the moment, we take values of the set 6 for the triplet 
odd state parameters. By virtue of the continuity of the phase shift from its low energy 
value, °0,*>0 and *0,7<0. Then in order to reproduce (a,),>0 we have g(€,) >0. 


There is another possibility in which 7(€,) <0 and °0,%, 90,7 >0 or °0,°, °0,7 <0. However, 
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these phase shifts do not continue smoothly to the low energy values. 
The phase shifts must satisfy the following conditions : 
(i) In order to reproduce the magnitude of a,, we have 


(a,) ,=0.28 —0.21. (for 95 MeV) 


Higher wave contribution which is mainly due to the *P-waves is estimated from the 
p-p polarization data as 0.02—0.03. 
(ii) . Moreover, because of 1/2 y(€,)<1.5, we have 


eae 35. (for 95 MeV) 
7 (€1) 

This quantity depends only on °0,* and °0,". 

(iii) To reproduce the magnitude of a, we have 


(a), 0.146 0.116. 


Higher wave contributions are estimated using values of phase shifts listed in Table I 
and Table II, set 6, and are found to be about 0.01. 


(iv) In order to reproduce the n-p total cross section, we have 
(sin 70,*)? + (sin 70,7)? =0.5 +0.1. 


The higher wave contributions, which are estimated from the p-p cross sections and Table 


I, are subtracted from the experimental value. 


Table II. YT=1 phase shifts (in radians) 


SP, SP, SP, 3Fy € aS, 


a(TMO) 0.377 —0.187 0.065 —0.035 ~0.704 0.825 
b 0.42 074 0.07 —0.04 = 0,26 0.40 
c 0.45 —0.22 0.07 | —0.04 ~0.35 0.36 
d 0.52 -026 | 008 | —0.04 —0.18 0.00 


d’ 0.52 —(:26 0.08 | —0.04 —0.35 0.27 
| 


First, we determine regions of °0,* and °0," which satisfy conditions (ii) and (iv). 
Then, by the condition (i), values of 7(€,), and therefore values of €, are determined* for 
each curve of (2/97(€,)) (a,))=const. These procedures are graphically shown in Fig. 3-a. 
Next, for each set of (°0;*, °0,7; €,), (a)o is calculated, and compared with the condi- 
tion (iii). Thus the regions of (°0,%, °0,7; €,) which are compatible with the n-p 
polarization experiment are determined. The shaded portions of the graph in Fig. 3-b 
are such regions of (°0,%, °0;"; €,). 

In the previous calculation” we took GA—= 0152, *0,.=—0.28- and €,=0.1. This 


point is outside the allowed regions. 


* Corresponding to one value of g(€,), there exist, in general, two values of €, among which the larger 


value is denoted as €s and the smaller €<. 
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0; 
= ss ‘= 
Z . (4,), 
&\é 
- 0.40 ¢(e,)~—1.5 — ¢,=0.17 40.09 
“9.60 g(é)=—1L0 ~ 1.35 ¢,= —0.16+0.09 or 0.51+40.09 
- 0.80—> g (e,) = —0.75~ —1.00- ¢,= —0.30 £0.05 or 0.6540.05 
-- 1,00 g (e,) = —0.60~ —0.79 > e, = — 0.37 £0.03 or 0.73+0.03 
1.20 ¢ (<,) = —0.50~ —0.66- ye = —0.40+0.03 or 0.77+0.03 
1.40 g (¢,) = —0.42~—0.56 > ¢, = —0.4440.02 or 0.80 + 0.02 
30F 
——> 
0 0.2 0.4 0.6 0.8 1.0 


=().2 


(4))9=0.12 3 (2) 
(4,)9= 90-11 5 (¢.) 
—0.4 ‘ (q&q)p=0.12 3 (e,) 


(49) =0.11 3 (s,) 


=—{\6 


(a,))= 0-11 3 (s.) 


(b) 
Fig. 3 


§ 3. n-p angular distribution 


In the previous section, we determined the regions of the phase shifts with J=1, 
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which can reproduce the n-p polarization data. In the present section, we shall investi- 
gate whether these phase shifts can reproduce the n-p angular distribution. 
The angular distribution do/d®2 is expressed as 
do i 
Saree B,, P;(cos @). 
A remarkable feature of the experimental data is that do/dQ is V-shape and is 


symmetric about 80° in the centre of mass system, and so 
Bp om nd. 8. 0. 


We have two allowed regions of the phase shifts determined in Fig. 3-b. In one 
of them, where |*0,*| < |°0,"|, B, is always negative in disagreement with experiment. In 
the other, |°0,*| > |°0,"|, we take several sets of phase shifts as typical examples which 


are shown in Fig. 3-b and Table III as Q—(®. 


Table II. n-p angular eistributions 


33,0 33,7 | é, 1g) | By | B Se ee 

fe oes | 025 ‘ oo ae Met | yee * re 1.76 

B 074 | 0.23 | A osiero20 5,86 | bey eect g oe | 108 

| 

@® | 083 =023 ae 017) 6.41 | aa 24 aie 1.58 

© eee oh onze MOOS onlwees0.17 5.98 | ares So de sPeonmuce 

® 0.67 ~035 | ee —020 || 5.87 | ea | we ag 1.64 

| | | | 

me oy4en|t odor) FOF) yer0.17 | 6500 Tn elon Ip actin) See Gaeiiee 

ec rn eer ee eae 


Assuming T=1 phase shifts in Table II-b and L>2 phase shifts in Table I, the 
angular distribution is calculated and the results are given in Table III. Some curves are 
plotted in Fig. 4. | 
Since the 'P,-wave phase shift is fairly dependent upon two-pion-exchange potential 
in the region II, its magnitude has some ambiguities within —0.25--—0.15, and 
determined by the condition to reproduce the total experimental cross section. ae 
For our phase shifts B, is rather large, and gives almost flat angular distribution 
between 60° and 120°. This aspect also appears in the experimental data. B, is 
always small and negative. These general trends of B, and B, are almost deterred 
by the higher wave phase shifts and should be regarded as one of the characteristic 


2 a es Pes 
consequences of the one-pion-exchange potential in the region oe 
B, and B, are sensitive to the mixing parameter €,. If €, becomes smaller B, 

2 
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becomes larger and the hump of the angular distribution curve at about 90° is 


suppressed. At the same time, —B, becomes large and gives rise to a large backward 


scattering cross section. 


10) T T Lj T T ai aaa 3 i a i ae a Tr rT T a a 3 : 
0° 40° 80° 120° 160° 0 40° 80° 120 160 
@ (c. m. system) @ (c. m. system) 


Fig. 4 n-p angular distributions. Some curves of Table III are plotted. 


To restrict the coefficient B, within 
a narrow region around zero, the value 
of €, must be limited to be about 0.1 + 
0.2, though in this case, theoretical 
curves have smaller forward scattering 
cross sections than the experimental one 
(see Fig. 4c). 

do /d® at large angles can be dimi- 
nished by the use of somewhat larger 
value for —*0,* and 'd,'. These circums- 
tances are exemplified in Fig. 5a. 

It seems impossible to increase 
theoretical value of do /d at small angles 


0° 40° 80° 120° 160° to that of experimental values by some 

@ (c. m. system) (deg) change of T=O state scattering para- 

meters. Therefore we are obliged to 

change values of T=1 parameters to some extent. Such slight changes are done under 
the restriction that they are consistent with the p-p scattering. B, and B, are sensitive 


to the mixing parameter €,, and larger values of €, give smaller backward and larger 
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forward n-p scattering. Changes of the *P-phase shifts are also effective to shift the 
angular distribution curve as shown in Fig. 5=b: 


mb mb 
ster f ster 
(a) (b) 
20 20 
154 iC 


0° 40° 80° 120° 1607 02 40° 80° 120° 160° 
@ (c. m. system) @ (c. m. system) 
Fig. 5 n-p angular distributions. 
a) Effects of changes in T=O parameters. Curve () is that of Table III with €,=0.09. Curves c¢ and 
c’ are calculated using T=1 phase shifts of ‘Table I-c, and 


191,= —0.25 for c 
35,%=0.60 38,7 =—0.28 €,=0.09 and 
1§1,;=—0.20 for c’ 


b) Effects of changes in T=1 parameters on n-p angular distributions. Curves are plotted for 


3§,%=0.62 36,7 = —0.25 and €,=0.18 
and T=1 parameters of Table II. 
c) The curve is calculated using phase shifts 


given in page 9. 


From these results, we may estimate 


the final values of the phase shifts, which 


reproduce all nucleon-nucleon scattering 
data at 90 MeV and consistent with 


0° 40° 80° 120° 160° the pion-theoretical potential. They are 


@ (c. m. system) as follows : 
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pe AE, me Died eee es sp *P, . *8,4°D, 
state 
0.00 —0.03 0.00 0.29 *4,.=0.08 —0.26 0.52 *0,7=0.74 
phase shift : 
39% = —0.04 0,7 = — 0.28 
€,= — 0.18 €,=—0.07 
state ‘Gs oP et? oe oe 
phase shift 0.00 —0.04 0.05 —0.17 0.20. 


Theoretical curves calculated with these values are plotted in Fig. 1 (polarization) 
and in Fig. 5-c (angular distribution) and are compared with the experimental data. It 
is to be noted that the quantitative conditions imposed by the 95 MeV polarization data 
in section 2 are too strict for scattering parameters at 90 MeV, since P(#) increases 
rapidly with the energy. If the conditions for the polarization are not taken seriously 
the fit of the angular distribution may be improved. 

Of course, these values of phase shifts have some uncertainty. The phase shifts of 
high angular momentum waves have some uncertainty mainly due to the uncertainty of 
the coupling constant. These uncertainties may have some effects on the results through 
the interference terms with L=1 and 0 waves. The lower wave phase shifts also have 
the larger uncertainty, which comes from the experimental error and the unfixed part of 
the nuclear forces in the regions II and III. Last values of the triplet odd phase shifts 
are a little larger than that estimated from the p-p analysis”. At any rate, we see from 
these results that the values of phase shifts, which are reasonable from the pion-theore- 


tical point of view, can reproduce all the nucleon-nucleon data at and, consequently, up 


to 90 MeV. 


$4. Coneluding Remarks 


We have analysed the n-p scattering data at about 90 MeV, and obtained phase 
shifts which are not contradictory to the outer part (i.e. regions I and II) of the pion- 
theoretical potential. 

In the outer part, especially in the region II (0.7< x< 1.5), the pion-theoretical 
potential has not precisely been obtained, although its general features are well known. 
Comparing our results with the TMO potential’, our phase shifts seem to suggest a 
little weaker fourth order potential (two-pion-exchange potential) than that of the TMO 
potential. This indication is compatible with the propensity of the FST potential and 
KMO potential.* 


Our phase shifts with L>1 are calculated or estimated under the assumption that 


* See, for example, Fig. 1 of ref. 6) ot Fig. 6 of ref. 7). 
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the potential inside one half of the impact. parameter (x<6,/2) hardly contributes to 
them. Our results seem to show that this assumption is actually realized. This fact 
means that the inner potential is not so drastic as to cover the effects of the outer 
potential. The inner interactions may be tentatively represented by “ hard core + square 
well potentials ”’. 

Our small values of °0,“ and '0,° indicate the existence of hard-core-like repulsive 
potentials of radius x,=0.2—0.3 in the even state potentials.”*) Such repulsive interac- 
tions should be regarded as intermediate energy properties of the nucleon core. On the 
other hand, we have some information on the low energy properites of the inner interac- 
tions by analyzing low energy data of the two-nucleon system, especially those of the 
deuteron. Detailed investigations on the energy dependence of the inner interactions are 


desirable since it might be connected with the structure of the nucleon core. 


This paper was written basing on the discussions with S. Otsuki and R. Tamagaki. 


The author wishes to acknowledge their co-operation. 
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It is the purpose of the present paper to discuss the intuitive meaning and the interrelations of 
various coefficients of the Ising model. Thus, the intuitive meaning of 6;,, and 67,,, the coefficients 
of cluster sums, is explained. Wan der Waerder-Oguchi’s expansion is generalized to the case with 
magnetic field and the coefficients @,,,, is defined and its meaning is explained. The equivalence of 
the calculation of the susceptibility with zero field and the non-crossing chain problem is proved. 
And a method of evaluation of these coefficients is proposed which requires least labour of counting 
the graphs by using the decomposition of Fuchs’ coefficients 7,. 

The data of our calculations are listed in the Appendix. From these data an outlook of this 
problem to the future and the theory of the imperfect gases are discussed. 


§ 1. Introduction and summary 


Ising model was first introduced as a model of ferromagnetics.'’ The problems of 
alloys and lattice gases were shown to be equivalent to this problem*” in their mathema- 
tical structure. For the problem of square lattice without magnetic field, the partition 
function was obtained by Onsager," Kaufman,” and Kac-Ward," and the spontaneous magne- 


tization by Onsager” and Yang.” 


The problem with a finite magnetic field is, however, 
not yet solved in a closed form.* 

To obtain the partition function in the case where the magnetic field exists we must 
either solve the eigenvalue problem” by the perturbation method" or count the number 
of graphs, and express the partition function in double series of the temperature and of 
the magnetic field. Domb et al.'"~" deviced a method with use of symmetry by which 
the labour of counting is decreased and they published many elaborate results concerning 
the quantitative aspects of thermodynamic quantities by calculating the double series- 
expressing the partition function. The method of expanding the logarithm of the partition 
function (per spin) in series, however, is better than the method of expanding the par- 
tition function itself. The calculation of lim (1/N) log S)+--can be performed straight- 


N+ oo 


forwardly by picking up the only terms proportional to N (total number of lattice points) 
after enumerating the graphs on the lattice plane which are defined in respective cases, 


and some of its coefficients have an intuitive meaning as explained in the text. Further 


Hence we do not know where the magnetization curve has the singularity as the function of the 


magnetic field, and what kind of the singularity this is. This is important to the problem of imperfect 
gases,°5—27) 


OO 
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it is easier to calculate thermodynamic quantities from this series. 

We may cite various methods of expressing the logarithm of the partition function : 
the expansion by cluster sum 6, (Fujita and Katsura™ to 6,), that by irreducible cluster 
sum /, (Miinster and Sagel’® to 84), that by Fuchs’ coefficient 7, (Miinster and Sagel'” 
to 7,), that by semi-invariant 4; (Kirkwood') to /,, Chang™ to /,), etc. 

In the present paper the meaning of some of these coefficients is clarified and the 
interrelations between the coefficients are discussed. In particular, we give the direct and 
intuitive meanings of 6,,,, 6,,, and 2,,,,, which are the coefficients of x” and f’” of cluster 
sums and the coefficients of generalized Van der Waerden”)-Oguchi®”’s expansion in the 
case with magnetic field (detailed definitions are given in § 3 and § 4) respectively. 
Then the equivalence of the calculation of zero-field susceptibility to the problem of non- 
crossing chains is proved. Moreover, Fuchs’ coefficients 7, are decomposed in two parts 
and one of these is evaluated exactly using the number of nearest neighbors. With the 
aid of this decomposition of 7,, a method to obtain the coefficients with least effort of 
counting the graphs is proposed. For example, the semi-invariants up to /, can be obtained 
without counting the number of graphs by this method. This method can be applied 
also to other lattices, but the application will be postponed to another opportunity. 
Various coefficients obtained in this way are summarized in the Appendix. The forecast 
of the possibility of expressing the partition function in a closed form, the theory of 
condensation, the ring approximation of the irreducible cluster integrals, etc., are discussed 


in the last two sections. 


§ 2. Preliminaries 


We consider a two dimensional square lattice whose total number of lattice points 
is N. The interaction energies between nearest parallel spins and between antiparallel 
spins are —J/2 and J/2 respectively. J>0O and <0 correspond to ferro- and antiferro- 
magnetic cases respectively. We denote the external magnetic field by H, and magnetic 
moment of spin by m(=g/9). 

Let 

x= eX — eo LIkT and yaace se (2 2 1) 
The partition function per spin is given by 
xt "A (xe, y peeling (Tr eK dsisg + CUise HN (22) 
Noo 


where s;(=+1) is the spin variable, >) 5; s; means the summation over all (2N) nearest 
bonds, Ss, means the summation over all (N) lattice points, Tr means the summation 


over all (2) configurations. 
The free energy per spin as a function of temperature and magnetic field is denoted 


by F and is given by 
F/kT = —log A(x, y) +1/2:log y+log x. (2-3) 


The magnetization per spin, M, is given by 
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(1+M) /2=y 2 log A/ Oy. (2-4) 
(1+M)/2 is denoted also by / (specific density in the case of ong pen gas) in this 
paper. The free energy as a function of temperature and magnetization is denoted by 


F* and is given by Legendre transformation. 


_M 
it se eles! coals tye shear J) qual 


log y (x, M) + logx. (2-5) 
ae Lo a ome © 2 


We have following symmetry relations. 


A(x, y) =yA(x, ~), (2-6) 
Jae 
gigs Shai iy fi doe (2-7) 
Ce en ee 
on ea Bee oe 2-8 
A(x, 1) (—, 1) (2-8) 


Here we explain the topological terminology used in the present paper (cf. Riddel 
and Uhlenbeck"”) . 

A graph is a set of points, lines and connected graphs which are not connected 
together. A connected graph is a set of points and lines connecting some of them, of 
which any two points can be bridged by appropriate connected lines. An articulation 
point of a graph is a point where the graph could be cut into two or more disjoint 
graphs by cutting all the lines going to this point. A star is a connected graph which 
has no articulation points. A closed graph is a star which has no points to which three 
or more lines end, i.e., a closed graph is a polygon. 


Evidently the next subsumption relation holds : 
graph D connected graph D star D closed graph. 


In “integral”? formulation of cluster sums and of irreducible cluster sums, ( (3-4), 
(5-1), Fig. 4 and Fig. 5) the points in the above definition are referred to particles. 
While in the direct formulation of cluster sums and of generalized van der Waerden- 
Oguchi’s coefficients (in Fig. 1. 2. 3) the points are referred to positions on the lattice 
plane, and hence permitted lines are those which connect nearest lattice positions only. 


$3. Cluster sum 
(2:2) can be rewriten as 


a on no Zz 
log A(x, y) =(1/N) log S IT e~ KOs: 5) ffe—CO—50) (3-1) 


From (3-1) it can be shown that one +— bond contributes x, one ++ or — — 


bond contritutes 1, one — spin contributes y, and one + spin contributes 1 to each 


of 


configuration in the summation \}). The non-vanishing terms after taking the limit 


ol 


lim (1/N) log S}---are the ones which are proportional to N in the summation. We 
V-> oo 
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define the coefficients of expansion 6, and 5, by 
los “A(x, 'y) = S" 5, y’ 
=Db Ay, 


and let 
b, = Sy Bim Xe 
m 
er 
ef b,,=6 
Qy=6 
yee = e-o-9 2 
re 4 Be 4 aay 
b o= —85 
+ by = 25 
ee-0 : Qy=~7 
‘Ne — 22. 
O-8-—® — 
Es mE e—C-® =2 tf 
=8 e—6 


bg =130 


Fig. la Fig. 1b Fig. 1c 


nm, 


The graphs to be counted in a) byiny b) bys and c) 2 


bonds is given by 2/—(m/2), m and m in a), b) and c) respectively. Fig, 1 
In a) / is 


shows the case of 2 bonds and Fig. 2 shows the case of 4 bonds. 


The number of 


the number of points (positions) and in b) / is the number of particles, 


a) the disjoint points are contained. 


exist. 
singly, doubly and tripply occupied points respectively. 


the number of end points and that of junction points with three branches. 


The details concerning signs are explained in the text. 


In b) and c) disjoint points should not 


In b) multiple occupancy is permitted and @, O and ©@ means the 
In c) n is the sum of 


IQS 


196 


ees Adtt 


ut 


bs 1¢= — 1651 
by 99 = 33609 
Fig. 2a 


Thus we see the meaning of b 
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oni 
8 


b,,=83 Q4=50 
Z 2 
8 8 
8 8 


BPPAB EET att f 


Q 317 
8 =o 
—4 Qe, =526 
521 
boy = — 831 Qu= r2 
Qi4=0 
b,,= 4299 ~~ 
Fig. 2b Fig. 2c 


im? ym is the number per lattice point of the graphs each 


of which consists of | points and 21— (m/2) bonds formed on the lattice plane, and in which 
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Fxg the nearest lattice points are all connected (Figs la and 2a). 

Seve Alternatively speaking, 6;,,, is a number per lattice point of closed 
2 : : 

graphs each of which consists of m lines, and whose area 


+33 3 is | in the unit of one lattice spacing” (Fig. 3). Here and 


throughout in the text a number per lattice point means the 


Big 3.91 67.4 coefficient of N in the expression as the function of the total 


! is the area in th it : : 
eo number of lattice points N. 


On the other hand, we apply the Urcell-Mayer theory of 


of one lattice spacing. 


m is the number of Ae 
Roads. the imperfect gas’? to the Ising model,” 6, is expressed 


by 1 


1 
integrand / : d 
2 z 3 


configuration 


il connected 


a ae IT fis (3-4) 


where S means the summation over all man- 


ners of arranging different | particles in N G23) —4f? _4f 
points allowing the multiple occupancy, and (@+G) aap afi 


connected 


>) IT means the summation of the products of =4fe af? 


fizs over all connected graphs (clusters) which 
consist of | points (Fig. 4). This connected 


graph is irrelevant to the structure of the 
lattices. Here 
0 4f? 


fuy=—1, i and j particles are at the af? 
same point in the graphs, . 
=f=x~’—1, i and j particles are at the 


nearest points in the graphs, 


; —1—12f? 1—8f+ 16 
=0, otherwise. f i 
= : RAVAN %3. £. =2—24f + 30s? 
When we express 6, as a polynomial Fig. 4. 


in x, the term of the highest degree of x 

corresponds to the graphs which consist of / disjoint points, and hence the highest 
power of 6, is x”. The term of the lowest degree in x corresponds to the graphs in 
which the number of bonds connecting / points on the lattice plane become maximum, 
ie., 4,/[ or the least even integer exceeding 4,//. It is denoted by —2|—2)/J | 
in Gauss’ notation. [x] means the greatest integer not exceeding x. When we express 


b, as a polynomial in f, the lowest term is f°, and the power of the highest term is 
2+ 2[—2)/1 |+1. | (3-6) 


by, is the cluster sum in the case where J=0 (or T=) and is given by —(—)’/l. 

From the results (cf. Appendix) we have found the direct meaning of the coefficient 
of 6, in f” as follows. We shall denote this coefficient by b,,,. First we enumerate the 
graphs consisting of m bonds on the lattice plane such that in this graph disjoint points 
are excluded and the nearest points need not be connected. The sign + or — is multi- 
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plied to each graph when the number of connected graphs which construct the concerned 
graph is odd or even respectively. Next we arrange the same / particles in each of 
graphs. In this arrangement the multiple occupancy is permitted but the vacant lattice 
point should not exist. Again the sign + or — is attributed to each lattice point accord- 
ing as the number of particles on that lattice point is odd or even respectively. Then 
the first and the last signs are multiplied to each number of concerned arrangements. 
The algebraic sum of the resulting numbers for given | and m (for 1<m<2/+2{—2,/1 }) 
is bim (bo=—(—)'/l). This manner of countings is illustrated in Figs. 1b and 2b. 


§4. The generalization of van der Waerden-Oguchis expansion 


In this section we consider the generalization of van der Waerden*”-Oguchi’s™’ high 
temperature expansion to the case where a finite magnetic field exists. From (2-2), 
log A—log x— (1/2) log y 


giv. oN N 


=> bog II ch K(1 +s, s, th K) II ch C(1 +5, th C). (4-1) 


By performing the product we have 


1 : 
=—— leg | ch’? X ch®e 

N og |c c 

Xd) G435)5 5 thK+ S$) (% 5) Ge 5) th?K+---) 

sg=ti <t,j> <tj> <é/j!> 

x (1 +5) 5 th C+ 335) % Su th?C+ ---)], (4-2) 


where 5} in the first parenthesis means the summation over all bonds, and S* in the 
<ij> 4 
second parenthesis means the summation over all lattice points. When we develop the 


product, the coefficient of th”K th" C is given by 
>it, 55) (Sy Sy) batt (S;¢m) Sj(m)) Sp Spy ~7 *Sz(n) . (4 a 3) 


Each term of (4-2) corresponds to a graph which consists of m bonds and n points and in 
which the nearest lattice points are not neccessarily connected. After performing summation 


> jnok vanishing terms correspond to the graphs which consist of either the junction 
set 


point of even bonds or the double point of one end of the bond (s, s,) and point s,. 

Thus we see: the coefficient of th”K th"C in the expansion of (1/N) log|---], 
which we denote by Qin, is given by the number per lattice point of graphs which consist of 
m bonds and in which the sum of the numbers of the end points and the junction points with 
three (generally odd) branches‘! is n (Fig. 1c and 2c). Here the nearest points need not be 
connected but the disjoint points should not exist. Especially 2.,.,;,=0, and Q,,, is the 
number per lattice point of closed graphs consisting of m bonds. 


} Short range order of Cus Au was treated by Oguchi.2® 
TT Such a branced graph is seen for example in for example Q4,;= —56. 


a a 


2 ae << 
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Hence : 
log A—log x— (1/2) log y=2 log ch K-+log ch C+log 2 
eS So 2 th’ the €: (4-4) 
n=0 m=0 


The magnetization M and the susceptibility y per spin are given by 


M=—th C—3)[22,, th C 
aE >I { seme (n se 2) Qn 49, =, th”*'C]th” K, (4 ; 5) 
kT y=1—th? C439 [22+ (—8.2 on +122y,) th? C 


>in 2) (nes 1) 2,,-2, 2 Dm 
n=4 


a(n 2) (4-7-1) 0,75. Vth” G/ th™ K. (4-6) 
In particular the susceptibility in the case where H=O is given only by 2,,,. 
Rivet ges hoe 20s th” Kk, (4-7) 


™m 


Thus it has been proved that the calculation of the zero-field susceptibility is equivalent 
to the enumeration of non-crossing chains (two end chains) of given length in the lattice plane. 


Since 6,,, is the number of closed graphs as explained in § 3 it is evident that 


a bi,= Dom (4 8) 


Z=1 


§ 5. Irreducible cluster sum 


The introduction of irreducible cluster sum to lattice system was originated by Fuchs.” 
The irreducible cluster sum 9, is defined by 


ts 1 ar 
Pf) a Ne (5-1) 


where >} II means the summation of the products of f;;s over all stars (irreducible clusters) 
which consist of k+1 points. The meanings of f,; and S are as same as those in the case 
of b,. 6, and (3, are connected by'”"”*? 


= Sie ars pe (5-2) 
2 k 


Np ! 


14+M . ' 
The magnetic field as a function of pa with , is an antisymmetric func- 
2, 


tion of M=1—2p. Therefore >}%, 0" should contain the term —log (1—/), and 
hence the coefficient of §, in f° is —1/k. Defining /,/ and /3,,,’ by 


B.= a +B: Be = OP em ider (aps) 


we have 
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ee We < 


configuration 


© ' ie 
se sy 
asparatt <0, oS thine eet 


24f4 — 24f% 0 
© 


Cae 


3193= —2+24f7+ 80f3+24f4 


Fig. 5 
log A (x 0) = —log (1 —p) ee 7 k - B,! pith, (5 -4) 
fe at i im k+1 
log y(x, ¢) =log e—log(1—p) —4 log x— S38," (*. (5-5) 


The “integrand” of (5.1) and the configurations which contribute to the integral 
are shown in Fig. 5. 
The result of calculations of 9’,,, shows that the lowest power of /,,,’ is —{ —k/2] 


in Gauss notation and the highest power is 


m* (k) ='2(k—1) +[—2V k—1] +1. (5-6) 
First several terms are given by 
k a oo Se ee a a a ee 
m* (k) Ls <Qoewdhe Se of qoePCyoidervenngsy 


These properties of 3,’ concerning the highest and the lowest powers in f” will be 
used in the next section to discuss the properties of 7’. 


$6. Fuchs’ coefficient 7», and its decomposition 


In this section we express the free energy F* and other quantities in terms of Fuchs’ 
coefficients 7, and decompose 7, in two parts 7’, and 7,’ and we give the explicit form 


of i 


The free energy F* as a function of magnetization and temperature is 
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F* Bx! 
ae log e+ (1—() log (1—9) +8Kp—S} PEO + Ig x 


ket ky] 


201 


(6:1) 


In general, the term 8K is replaced by 2zKp where z is the number of nearest neighbors 
of the lattice. The sum of the third and the fourth terms on the right-hand side of 


(6-1) should be expressed as a function of o(1—p) =(1—M?) /4. 
k-+-1-th coefficient in the expansion in powers of (1—M”) /4 by 7. 
—8Kp si Bx eo ah Tk el LAE 
( ups me, = ra ( (1—)***. 


Hence /7;’ is expressed in terms of 7;. 


ye ee Oe)! 
kR+1 ofbe@+2 AL (kR+1—2A)! 


aa. 


Let 


We 


Then 


as Vs a 
pii—p) — 5 i ee ee Ye= 2D Vin tine 


2 


Inserting (6.4) into (6.2), we have 


ao a, 
<1 km m 4k+1 
Sy ee 


co Am mem oO m*(k) / m x 
patos Mi LS T. 1) is nes a Maer io oAh gard oes) ou Pics ae, 1 VA heady 


eS m D md pet 2/2] k+1 2 


denote the 


(622) 


(6-3)7 


(6-4) 


(6:5) 


The coefficient of ¢**'f” in the expansion on the right-hand side is 7,,,/k+1. 
We denote the part due to the first term of the right-hand side by /,,,’’, and the 


remaining part by 7;,,’. 


Then 
fant BL } Jia yh 12 Cw) Gee)! 
2ni pee m  (k!)2 
=O (m=0, 1,:--k). 


Yum does not vanish only when m=k+1, k+2,---m*(k). 


(pS EOS 


(6-6) 


This fact is due to the properties of the highest and the lowest powers of /3;'. We 


rewrite 7; as 


es (k!)? RN? m=1 m m=k 


Hence 7, is decomposed into 7,’’ and 7;’. 


! [Nae m—1 fC m mi (Ie) 
8 (2k)! & 5 (2k)! ¥% (= ) 6 re “ "tat 


(6:7) 


+ Fuchs®) proved (6.3) for k=3. (6.3) is, however, valid from #8, when we consider that yz’ exists 


from 7 and that (6.2) holds. 


202 S. Katsura 


The part of F* due to 7;/’ in the sum in (6.2) can be given in a — form. 


etd +1 
’ in the expression J a 


Of —k+1' 
x (1—p)**1, change the order of the summation 2 and integration, and perform the 


We apply the integral representation (6.6) to 7% 


summation Ss] under the integral sign. Then we have by residue theorem 


m1 


3 - Te ket pa _1-V14+-4fp—P) 

Fea ktr © AO =~“ Fat) 

= 1—™" (hi os 4/1 eed (6-8) 
ieee. erie : fs ard 


Noting d/dK=4(1+f)d/df, then 


Sy a —pyt=al (1M (eh Flog +f —M*))} 1K (6-9) 
k=0 vt 


=2 log (/1+f(1—M?) +1) —2 log 2 
—(1+M)log (V/1+f(1— M2) +M) + (14+M)log(1+M) 


— (1—M) log (vy ‘1+f(1—M®) —M) + (1—M)log(1—M). (6-10) 


The 7,’ contributions to the free energy, the magnetic field, and the susceptibility, 
expressed as the functions of magnetization, are summarized in closed forms and the 


whole expressions are given as follows. 


1+M , 1—M 1—M 
Ege = ligt i cence 
kT 2 ce ee = =e 
ae (+fQ—M))MdK-SUT# (ISM 66-11) 
. hres “r= 
aS tag 
IT gy 
=e At fA SM) +M | 2 i+ eS a! —M? o 
V1+fa—M) —M 2 vs 1—-M 2% Sein 
Nae Oe 1 ay. 
kTy =(1—M"))/1+fC9—M*®)  1—M 
1S (1-M\ MS, (eee 
+ 3n (— == ra Ss hee : (6-13) 


When we neglect 7,’ contribution, the Curie point is given by x,=0.5, which is as 
same as the Bethe-Fowler-Takagi approximation. 
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S$ 7. Coefficients. /” jm Gad |" 


jm 


In this section we rearrange the series S)S)- km 


th cet beet 


coefhcient /",,, and relate it to the number of graphs in the lattice plane. 


1 — V2 kh+1 
reef f” and define the new 


Since the highest and the lowest powers of [,,' are —|—k/2| in Gauss notation 
and m*(k) respectively, the highest and the lowest powers of 7,,,’ are k+1 and m*(k) 
tespectively. When we rearrange the double series 


| Ba — M2 k+1 
Se / km (- m = 
hack eI 4 ) i ae 
into 
>(—4 2) = Veber get cis =M") te 72) 
j=0 1—M i m*(k) Roe th é 


then the series in the second summation symbol becomes a function of only one variable 


fon. We expand this series in porwer series of 
u == th [4 log{1+f(1—M?)} ] 
and denote the coefficient by /’’,,. Further we expand this result into the double series 


of and u, and denote its coefficient by /’,,,. Then 


1— 


(7. 1) = (7-2) = +)= Py! (7-3) 


Ss (Le So 7-4) 


g=0N 


When we put M=0 in (7.4) and insert it into (6.11), F*/kT should reduce to 
the Onsager Kac-Ward solution, and hence /',, is identified to the number of closed 
graphs (=), defined in § 4). 

We could not find the intuitive meaning of /”,,, or /’,,, for j = 1 in general. We 
find, however, that the first non-vanishing /',,, for given j is the number of graphs (in 
this case stars) which consist of m bonds and of least number of points. 

Free energy, magnetic field, and susceptibility can be expressed in terms oa BFS 
and the part due to /%,, can be expressed in a closed form. For example the suscepti- 
bility without magnetic field is, 

peers, 60, —thK i 1—3 th’ K 


iT y  1+thK (1-+th K)* 
4 (1+th? K) (2 th* 2K—1) 2 2 K(kK(K)) 
(14th K)? 


—8>' Pig, tat KS, (7-5) 


m=T 
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where K(k) is the complete elliptic integral of the first kind, and k=2 sh 2K/ch* 2K. 


§ 8. Evaluation of the various coefficients 


The evaluations of the expansion coefficients b, -£,°25 Ferl yee, eee 


valent to each other, ie., a knowledge of any coefficient to m-th power of respective 
parameter (x, f, K, thK; C, y, th C, etc.) enables us to calculate all other coefhicients 
to m-th power of respective parameter. 

Now we propose a method of the evaluation with least labour of the enumeration 
of the graphs. 

First we consider the expansion by /’,,,.. From the meaning of Pom (the number 


of closed graphs) or from Onsager’s solution we see 


ee es P.=2=T 4 (8-1) 


and 
f= 0 for ni le 2s Se oe 


And non-vanishing /', begins from m=7. 


Thus from (7.4), 
Tsa =P oq =4, Tas =P’ =—40, Hse’ =F q' =276. (8-2) 


Using 7,’ given in (8.2) and 7,’ in (6.6), all (,,,’ and 6,,, up to m=6 can be obtain- 
ed from (6.3) and (5.2). The semi-invariant /; can be given by the relation 

Ss] gv (1 —p) May (—4K)* 3 

= ki! | ii i! — 
Thus §;,, 6, 4,, can be evaluated up to 34, b,, Ag respectively, only from data (8.1). 
If we regard the Onsager’s solution as known, the above coefficients can be given without 
counting any graphs. 


Further, we know (directly or from Onsager’s solution) 


P=, Py3=9 2, Po=0. (8-4) 
The first 7’,,, which does not vanish is /” 


2910° 
? a 7 
i Thus the data of / i> Uys, and J’ enables one to calculate all coefficients up to the 
ot ower of —" 9 4 
Pp the temperature parameter (f°, K°, th °K). The necessary 5,,,’s are 
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‘geese - 


bs514=134 (the number of graphs which consist of 8 points and 9 
bonds in the lattice plane). (8-5) 


The direct counting of graphs corresponding to these 5,,, enables us to calculate /’,,, 
Pis¢ Py, and all expansion coefficients 7/jm, ns ims nm ip tokn== ern D hus bye ye 
4;, are completely determined up to 6,, 9,, A) respectively. In the Appendix the values 
of b,, Bx’, 7, 2, 0/, Pj and 4,, are listed. The values up to 6; agree with the values 
obtained by the direct counting of graphs by Fujita and the present author-'? Further 
we transform Domb’s coefficients (summarized in Appendix 1 of reference 11) into 6, 
and compare it with ours and ensured the coinsidence up to /=8. , and 7, are ob- 
tained by Miinster and Sagel’ up to k=4, and we find their 9, and 7, should be cor- 
rected. The semi-invariant /, are obtained by Changt™ up to 4,, which agree to ours 
up to 4,. 

In conclusion we can say: In order to evaluate the expansion coefficients up to m'th 
power of temperature variable, (f, K, th K etc), it is sufficient to know /’,,, up to m=m’. 
for m! =m =1+1 
and hence to know the first /—m! (>0) terms of 5,,, by direct counting of graphs. 


When we regard /%,, as known it is sufficient to know 9’,,, or 6, 
This labour of enumeration is much less than the one which Domb’s method™ requires 
in the evalution of the coefficient of A(x, y). 

Here we have used the symmetry of F* with respect to ¢ and 1—p, the number 
of nearest neighbers, and Onsager’s solution in the form of high temperature expansion. 
Besides we have the knowledge of Onsager’s solution in low temperature expansion and 
Yang’s spontaneons magnetization. When we use the coefficients of these expansions, the 
number of graphs to be enumerated can be decreased further. However we did not use 


these expansions since a systematic method could not be found. 


§ 9. Discussions 


Well-known beautiful results of the Ising model of the free energy and the magneti- 
zation have been obtained in the case where the magnetic field does not exist. In fact, 
it might be a fortunate accident that the partition function without magnetic field is 
reduced to a combination of elementary functions, and we might not expect such circu- 
mustances in the case with finite magnetic field. 

We notice, however, all existing theories (Kramers-Wannier, Onsager, Yang) follow 
the method in which first the partition function of finite system is obtained and at the 
last stage the limit of infinite volume is taken. It is desirable that at a very early stage 
of the formulation the limit of infinite volume is taken and the functional equation for 


+ In Chang’s general results’, 0=06’, z=4, y=8, 71=0, 12=24. 
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the partition function per spin is established. This functional equation may be an integral 
equation or a difference equation. At that time, the coefficients of the double series ex- 
pressing the partition funtion may be obtained by a general formula or at least by a recur- 
rence formula, even if the solution of this functional equation is not reduced to a known 
function. Then we may regard the problem as solved for the first time. 

The method developed in this paper would be the first attempt toward this end. 
There is obtained an expansion of which the coefiicients, such as 5,,,, 6,,, and 2,,, have 
intuitive meanings as described in the preceding sections. In Figs. la, 1b, le and 2a, 
2b, 2c we compare these different manners of counting graphs in the case of 2 bonds 
and 4 bonds respectively. 

Unfortunately we could not find the direct meanings of PP pe Tt T5'  e 
except /’,,, and [’;,, which does not vanish for the first time for given j. It is desirable 
that some expansions will be found in which the coefficient is expressed only by the 
number of stars in the lattice plane. or example, the decomposition of 1%,,, may be 
hopeful. A little promising method to obtain the solution in a closed form is to express 
the partition function in a determinant form whose expansion coefficients are either: 1) 
the number of closed graphs 6,,,, as a function of the number of sides 2/—(m/2) and 
area /, or 2) the number of graphs 2,,, which consists of m bonds and whose sum of 
the number of end points and the number of branch points with three branches are 2, 
ot 3) I’m for j = 1 explained in § 7, etc. The determinant which is formed of any 
of these will be the generalization of the Kac-Ward determinant. The attack along this 


line is expected in future.* 


§ 10. Note on the theory of imperfect gases 


As described at the footnote in the introduction, we do not know whether the 
condensation point is a singularity or not in the problem of imperfect gases. It should be 
noticed that the Yang-Lee theory of condensation’ does not answer to this problem.) 
Here we remark a few points from the results of our calculations. 

In the distribution of the signs of the coefficients of the data of b,,, and 6,,, regu- 
larity is hard to be found. In )’,,, and 7’,,, the terms for odd & are negative and those 
for even k are positive up to yy’, 7x’, and up to 9'-my 72m: These properties, however, 
are destroyed at f’,,..=+9, 8’y,;,=—60. The signs of the coefficients of 2,,, and 
Aen are distributed regularly in our results. 

The 6, as a function of temperature is shown in Fig. 5. 6, (up to 6.) is positive 
below the Curie point, and is an oscillating function above the Curie point. In the 
Husimi-Temperley model, this property has been proved.?” And it is surmised that zero 
of 6,(T) has a lower bound also from Kihara-Kaneko’s results” for the cluster integral 
of the gases with square well potential. It seems from the curve of 6,(T) in Fig. 6 


*Temperley has informed the author that he has attacked along such line (private communication). 
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that the higher 6,’s do not contribute 
above the Curie point. 


On the other hand, 2, has no 
such a property, that is, there exists 


such /3,’s which are negative even at 
sufhciently low temperatures. Thus the 


usually accepted proposition, “cluster 


integral is always positive for sufficiently 
low temperature’, is not valid for ,. 
Hence the relation between the con- 
densation point (denoted by y,) and 
the singularity of S$) 8; o*, which has 


the smallest absolute value (denoted 
by (;), has the next six possibilities, 


and we cannot select the one in the 


étye present stage: 1) The singularity 


does not exist for j0|<1 (the case 


of Husimi-Temperley model*”. 2) 


te The singularity is not real and exists 
2 outside the circle of R=p, (the case 
of ideal Bose gas)**? 3) The singu- 


larity exists on the circle of R=, 


100 


200 F but is not real positive. 4) (7, is a 
axb§ ; 

a singularity, but there exists another 
L 

singularity on the circle of R=. 5) 


1 es! . ° . . 
) a singularity and there exists no 
0.2 0.4 0.6 0.8 it), 04 & a! 
singularity on the circle R=, except 
Fig. 6. Temperature dependence of the cluster sum. 8 y ’ i ee P 
Tie shecia denotes cent IAT 0, 6) The singularity is not real 


6°—68 should be read as by~bg. and exists inside the circle of R=/(,. 


Ring approximation for /%,: In Mayer’s theories of imperfect gases, the evaluation 
of the higher terms of , is generally difficult, not only in the integration but also in 
writing down the integrand. Therefore the approximation of (7, to the ring part of the 
integrand due to Montroll-Mayer” is often used. When we compare the exact §, and 


the ring approximation of (3, 


1 


o, = —St4pt > : 


ass Ae 
3 at 


Byte =— + 12f* + 10f%. 
The approximation is not good. Thus the solution of Green’s linearized integral 
equation” for the imperfect gas, which is equivalent to the sum of the contribution of 


8, in the ring approximation,* would give a fairly different isotherm from the real one. 
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Appendix 
Cluster sum 6, in powers of x. b= 3b, 
b,=%, 
bate 8 | 
2 
b, —— 6x®— 16x + 3132 
3 
Bie x° + 18x"°— 85x"? + 118x'4— 209 36 
4 

by = 8x"°4.43x"—400x44 926x%— 872x184 1476.20 
bus 2x4 40x" 4 30x" 165124 17978 8__ 91 443 

3 
a 22x"? +-136x*— 486x"°— 5644x'5+ 33609x” 
be 6x? +134x"4+ 389x"°— 3986x—13323x” 
i = a a fe ee 
b= 30x44 s<eees 
b= x" 4.25. 599 
b= 2b weaves 
b= 68x84 <<... 
b= 22x°° 1. ost aun 
b= Cee eee 
b= ade = Sestuien 
b= B8x08 +. «025. 
b,3= 30x84. ...0:. 
b= ged SP 
Bia ce 5 gat MRE 
Dy >, Bry X= Xt 28 4 << + 12x"? + — + 130x'*+...... 
Ss) 


> U Bj, X™ = 2x Bx! + 34x54 152x471 4x24 3472044 0.00. 


i} 


Cluster sum 6, in powers of f, b= Stbim Ff", f=Hx?—1, (b,=6, x”). 


b,=1 ; 
6,=—1/2 +2f 

b= 1/3 —4f +4 6ff , 

be= AY 4 OF se oirai SF Rah ey 


b= 1/5 —8f + 64f— 148f4  93ft4 sf 
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10739 4, 
EE, 
6 


+ 6520x7 — 


— 756402 + 85954x2%°—49328x° + a 


601905 » 
aie eee AE 


+ 164790x” 


— 1064393,24 | 37670x° —785091x% + 377040x" — 
2 


(A:1) 


(A-2) 
(A: 3) 
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b= — 1/6+10f 
b= APT de. 
b,=— 1/8+14f 


= 1/9— 
b= —1/10+ 
by=o 1/1l— 
b= —1/12+ 


Coefficients 2,=>'2,, th” 


2,= 


Q2,= 2th K+ 6th? K+ 18th®* K+ 


2,= 
2,= 


—7th? K— 56th? K— 317th*K— 
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— 130f?+ 831f'+ 312f+ 54f° 


16763 __ 
es 
+ 230f?—1572f°+ 4299f*— 4492f*+ 


= if? +3682f¢—— ft 30434f°— 23263f° 


1082f° 


K, 


th? K fs 2th® K 


50th*K+ 138th°K+ 370th°K 
1512th® K — 6493th® K 


116 


: “—~th® K+ 526th‘ K+ 4484th® K + 30248th® K 


521 


_ =o K— 5072th® K —58238th* K 
9812.45 K +.49992th° K 
__ 47644 16x 
3 
(apes 
m 


Irreducible cluster sum 3,’ =3,+ - =S5P' em f™- 


8, 
oa 


Bo! 2 
ts — cae 2 2 
; f 


Bs —_ fi 


= 2f 


10 “ 
aay f 
— 4f—11f*— 3f* 


284 


Sf 204 284 *4+ 0 46fP+ 2f? 
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ie ing 
+ -312f' +9 22f* 
Oe era + 194f°+ 6f 
+ 84f% + fe 
os+ 1. 39% 
-+8f° 
2p 
(Cee) 

Os 10 112 1 49 hi4K-+... 
at ee ti +12th”’K es K -H130¢ ai 

2 


+986th’ K +2586th® K+6746th® K-+--- 


ae (A°5) 


(A-6) 


Zant 


PHN 
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Pope —14f*—128f*— 320f°— 252f'— 34f* 
7 
. 12828, 2937 ». 
ir 7 AL 134f 4 B55f+——f > —— 
8 2 7 - 
21831 
oe ff ee f 
9 3 2 
Vlop 685s 26864 os74fi4 37646f° 
10 5 3 
Bro! —360f*— 9660f'— 75012f% 
11 
a 60 f'+ 5822f7+ 93674f° 
12 
Bip! _ — 1944¢°— 74921f° 
13 
Pa 1944674 37458f% 
14 7 
If 
Beso. — 10702f* 
15 
ae mies F 
16 4 
Pie) _. 
17 
it 
18 
os) 4 as 1 \& 2 © p’ 
as ‘om = — ee — > BY vn = a > >} = 
>t B'm=0 (m=2.3.4-:-) 
Semi-invariant 4;. t= (1—p) a 4:/N in Kirkwood-Chang’s paper is 
A,=2t 
A= 2° 
As —2P°+ 8e 
‘= 2°— 24°+ 84¢ 
4= —2°+ 56°— 600+ 1632 
Ae 2°—1208+ 30604— 21120°+ 432004 
A= —2+248—13800¢+ 178080— 856800+ 13996808 
4,= 2° —504t°+ 588844 — 124320084 105033604 — 392313608" 


=) ro) =; é k 
—=loge + S)S° ag (-'(2) = K*+ 2 K%4 32K 4 425 p84 Soman 
t=1 i! 4 6 45 252 


=1 keg 
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ae 354f°+ 6f*° 


— 9216f°— sactiak +f? 


een e Sra 


+. Rive see == 13 
16 ie 


a 256720f°— Bien ae 


3 
2577344 
— ; fom 


+ 851696f°+------ 


pt Ad ORCS, alee 
3 
724568 ro 


3 


ae 
== 60316f°—------ 
60316 pe Sane (A-7) 
o 


==(=)" =, (m=1, 2, 3---) (A-8) 


(A°9) 


denoted by 4, in this paper. 4;,=>)/,, t”. 


(A- 10) 
(Amity 


+53938080¢° 


214 
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Coefficients 7:/=>)7'um f patt of Fuchs’ coefficients. 


Hee —240f' —34f* 


i, 2461 34 354f°+ 6f” 
2 
10 39 eu 
eis — 6384 f° + term of order f ad 


i. term of order f"°------ 


Coefficients f=) 1,’ wu” and [;= >) Tj, wu", where u=th[ (1/4) log{1+f 


v= wo+2ue —8u' +a —80n'------ 
= 9) = Ga Og oes 
= 6u'° 
ome 
r= wh +20 = oa +12u"° 
Lr; — out = GP a esse 
lr, = 6n° 
i= 


4, 201 fe co gemias 
a 
at? 
112 w» 
——— 
=, 
201 
4. are 
2 
ni! 


In conclusion, the author expresses his sincere gratitude to Prof. Y. Nomura, who 


encouraged and helped the author in many respects. 


T. Oguchi for valuable discussions. 
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Effects of virtual nucleon pairs to the electromagnetic structure of the nucleon are calculated 
covariantly in the charge independent pseudoscalar meson theory. A rigorous expression for the vertex 
operator is derived in 22, which is strongly reminiscent of the old-fashioned first- and second-order 
perturbation theory. Our method of approximation is to replace the complete set of eigenfunctions 
of the total Hamiltonian by some suitable functions, which are chosen in this paper so as to take 
into account the main contribution from the virtual nucleon pairs. These effects are proved, regardless 
of the magnitude of the coupling constant, to reveal itself only through the renormalization factors, 
which affect the relative magnitudes of the isotopic scalar and vector parts of the nucleon-current 
contributions and the meson-current contributions. Numerical results are almost the same with those 
obtained by cutting off the momentum integration at the nucleon mass in the second order pertur- 
bation calculation both for the ratio of the anomalous magnetic moments of the proton and the 
neutron and the mean-square radii of the charge distributions. 


S$ 1. Introduction 


Recently, high energy electron-nucleon scattering experiments have raised some im- 
portant problems concerning the structure of the nucleon’. The fact that the mean-square 
radius of the charge distribution of the proton is rather large while that of the neutron 
is almost zero, which seems to have been established experimentally”, indicate decisively 
the inadequacy of the static pion field theory even as the first approximation, since it 
gives equal and opposite radii of charge distribution for proton and neutron respectively. 

There are also several other low energy phenomena, which have not yet been given 
satisfactory explanations and are beyond the scope of the static pion field theory; i. e., 
anomalous magnetic moment of the nucleon, S-wave pion-nucleon scattering and so on.* 


In this paper we investigate the effects of virtual nucleon pairs on the electro- 


; * After the completion of this work, several works on the dispersion theoretic approach to the S-wave 
Pion-nucleon scattering have been published. See, G. F. Chew, M. L. Goldberger, F. E. Low and Y. 
Nambu, Phys. Rev. 106 (1957) 1337, and K. Tanaka, Phys. Rev. 108 (1957) 1629. 
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magnetic structure of the nucleon. We shall not take into consideration the contributions 
of the strange particles in this paper; since it seems desirable to establish how to treat 
the effects of virtual nucleon pairs for the evaluation of the contribution of the strange 
particles also. So we investigate, on the basis of the relativistic pseudoscalar meson theory, 
the effects of virtual nucleon pairs to the anomalous magnetic moments and the mean- 
square radii of the charge and magnetic moment distributions of the nucleons. Although 
there is a possibility that some of the above mentioned phenomena might indicate the 
limit of the applicability of the present local field theory, it will be necessary, even if it 
were the case, to deduce from the present field theory some critical results rigorous enough 
to establish the inconsistency between the present quantum field theory on one hand and 
the experimental results on the other. 

We shall derive in § 2 an expression for the electromagnetic vertex operator of the 
nucleon, which is correct up to the first order in the electronic charge but is otherwise 
rigorous, and is not specialized to the kind of the virtual field and the type of interactions. 
The equation may also be written in a form strongly reminiscent of the old-fashioned 
first and second order perturbation theory. 

In § 3 our method of approximation will be presented, which consists in substituting 
the complete set of eigenstates of total Hamiltonian by a suitable approximate one. 

In § 4 approximate calculations are performed, in §5 the renormalization constants 
are evaluated and in § 6 numerical results are given for the anomalous magnetic moments and 
the mean-square radii of charge and magnetic moment distribution of the nucleon. 


In § 7 our method of approximation and results are critically discussed. 


§ 2. General expression for the vertex operator 


We will start from the following equation, which may be taken as a defining 


equation for the vertex operator, RC ae ae 0 fe 


—(T (x) $(y) 4, ) Yo 
=ieG ,, (SE) G (xx) I, (xy! 5 EG (y'y) , (2-1) 
where 
G(xx’) =i(TY@) PO) Do, 
G COVER TA. OAS) Do, 


) and A, are Heisenberg operators for the nucleon and the photon fields respectively, 
and (T--:)) means the true vacuum expectation value of Wick’s T-product”. Summations 
or integrations should always be taken on the repeated suffices or variables. 
Although all the formulae in this section hold both for renormalized and unrenor- 
malized quantities, we shall use the renormalized quantities throughout this paper. 
Operating on the left side of eq. (2:1) S,=70,.+M and (]]},” to the left and 
2D, =70,+M to the right, we obtain 
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~— DB, TY) $4.) o 2, 
=(T 4 (x)G(y) ju) ote Zot p7p9 (2—£) 8 E—y) 
+ied(x—€) ayy, {G (xy’) >* (yy) —G (xy) 6M} 
tied (F—y) {>0* (xx/)G(x'y) —OMG (xy) } 7a» (2-2) 
where 
n(x) =DH@), 
I(M=$(y) B,, 
ju ©) = —-C) 4, F) (2-3) 
and use is made of relations 
OFo—yo) Lin OG, (WriHie (F —y) GE) aru 
D (%—Fo) [Ta (X), ju (FE) |= ie PM (x—S) A (S), 
and 
19 ZZ "757 CT 4b: (x) 4 (x) f(y) ro=FDd* (xx) G(x’y) . 
Multiplying the free nucleon wave functions a(x) and u(y), which satisfy 
a(x) D,=0 and D,u(y)=0 
from the left and right of eq. (2-2) respectively, and using 
1G (xy’) S* (yy) —G (xy) 0M} u(y) 
=G(xy') Z,[8(y’—y) By] u(y) =Z.4(x—y) u(y) 
and 
(x) {>1* (xx/) G (x’x) —OMG(xy)} 
=i Z,[D,0(x—2/) |G (xy) =Z,a(x)O(x—y), 
one obtains 
— a(x) D,.OE(T $(%) 9 (9) Au &) Yo Dyu(y) 
=teLyit(x) 7, u(y) O(x—F) 0 (E—y) 
+(x) (T7 (2) 7 (y) ju) You (y). (2-4) 
We porform the same Operation as above to the right side of eq. (2-1) to obtain 
teu (x) 1), (xy 3S) u(y) 
up to the first order in e. Then we obtain 
teu(x) I’, (xy 3 €)u(y) 
=teZyii(x) 77, u(y) O(x—F) 0 (E—y) 
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+a (x) <T 9 (x) 9 (y) ju) ouly)- (2015) 

This is our fundamental equation. One may also derive this equation starting from the 
S-matrix in the interaction representation for the scattering of a nucleon by an external 
electromagnetic field, as was done by Low® for the case of boson-fermion scattering. 

One may also derive the following equation from eq. (2:5), which may be of some 
use. Choosing plane wave states with momenta p’ and p as the wave functions a(x) and 
u(y) respectively and expressing the space-time dependence of Heisenberg operators, 
O7(x) =exp(—i Px) Oj,(O)exp(tPx), where P, is the total four momentum operator of 
the meson-nucleon system, explicitly, we integrate over the variables x, y and € after 


multiplying exp(ik&) to the both sides of eq. (2-4), to obtain 


ien(p’) I, (p', p)u(p) —teZ,u(p') tr, u(p) 


TORCH ALAC ATA DICARA DY 
sks ™ (En Ps te) (En Py te) 


Pyn=VP! Pm=P 


Ce; Ve vs) Cis jules) (Ge Vii ay) 
m (E,+po—t €) (E,, + po’ —i€) 


Pr=—P Py=—YP! 


Ce es alas a) 4 Ve he Py) 


1 eR 
fs 2 (En— fy =i ©) (Ears kyrat €) 


n ™m™ 
Pyr=P! Pm=k 


Lee 28 wey CTE: ve) Cis; Vin FAK Lie. ju o) 
p my ek (En + fo=t€) (E, Skate) 
= Bae bs: (Po, jew) Cas Vi Ue) Cg) 
pep top ate EG amgi aE npeticni®) 
—»S} sy (Fo, iP a) (eo ae Ps) (Pe Vy! PO) ; (2-6) 
ay oe (BE, eee) (E+ alate) 


Dn=—k Pm=—P! 
where k,=p/—p, and 
Vi} (x) =a(p') (x) =Z." |" (x), a (p') I, 
V,(y) =] (9) u(p) =Zo" [a* (p), H'(y) J, 


a(p’) and a*(p) are the annihilation and the creation operators for a auclegn with 
momenta p’ and p respectively, H’ is the interaction Hamiltonian density, and # nS are 
the complete set of the eigenfunctions of the total Hamiltonian peouging to CEEAOEr 
E,. In eq. (2:6) ju, V, and V,; mean corresponding operators at pee orem of thé space- 
time respectively. Apparently eq. (2-6) has the form strongly reminiscent 2 the old- 
fashioned first- and second-order perturbation theory although there are following differences : 
(a) Every possible intermediate state Y’,, which is allowed by the conservation laws will 


enter into the sums. (b) The matrix elements must be calculated between exact eigenstates 
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of the Hamiltonian rather than between first approximations to these eigenstates. One 
may also define a suitable analytic continuation to the complex energy plane of = (2-6) 
analogously to the case of static pion-nucleon seattering”. It is most convenient, for 
this purpose, to take a coordinate system where k,=0, which is the center of mass 
system in the case of electron-nucleon scattering. We define L7,(p’, p;z) by the quantity 
obtained by replacing V,,}, V, and py=P,' by V,;, V, and a complex variable z respectively, 
in the right side of eq. (2-6). Then, 1,2 (p’, P32 has trivial g- and q’-dependences 


and satisfies the following relations : 


Ira (p’, P; z) =I (Pp, Pp i) (2-7) 
lim tha (p’, P > Zz) 
2> Byte 
=ien(p) I. (p’, p)u(p) —ie Zoi (p') = 7." (p)- (2-8) 


One may use these relations to evaluate the electromagnetic properties of the nucleon 
approximately together with some approximations such as neglect of higher configurations 
in the summation over complete set in eq. (2-6). We do not, however, attempt such 


an approach in this paper. 


§ 3. Approximation method 


We calculate the effects of virtual nucleon pairs to the electromagnetic structure of 
the nucleon starting from eq. (2-5) or eq. (2:6). Our approximation is to replace in 
these equations the eigenfunctions of the total Hamiltonian, ¥,, by some approximate 
functions, which should be chosen so as to take into account properly the main dynamical 
properties of the system. 


Now we assume the following Lagrangian for the nucleon-pion-photon system : 
&=—4Z{[4, (7 (@—iet, 4) +M) ¢] 
—[@7(0+ier, 4) —M) ¢, $}} 
—4Z,{{(8,+ieA,)¢*, (8,—ie A,)¢} 
+ 2° (¢* 9+ 9 g*) + (A,dbs)°+ 4° d,"} 


—¥R- oZld rs 


+42Z,0M[¢, $]+4Z,0 2 62+42Z2462 67 > (3 as 1) 


where Y= (4,—i¢,) /)/ 2 is the charged pion field operator, and the renormalization 
constant of the photon field (Z,P™°”") is put equal to 1, since our considerations are 
restricted to the first order in e. Then, ju(©), 7 (x) and 7%(y), defined by eqs. (2-3), 


are given by 


ju(§) =—-219, t7n¥l— 7 eb Tis Pu Gs 


te ee 
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2 (x) — ig re ees Vs v7 d; b = 0 Md, 


9(y) = 19 Z,Z.7 $7 57:6, +0MG, (3-2) 
where 
oh ae 
Lem SU OnO. | 
erie OOM! 


We will assume in this paper the following simple forms for the approximate state 
vectors as substitutes for the true state vectors ; 


PS exp (18) @,, (3-3) 
S=[Frarseb Ml pao, (3-4) 


where d= (@;)'”, @,’s are the eigenfunctions of the free Hamiltonian, H,, with the 


same eigenvalues and boundary conditions as Y’,’s, i.e., 
AY ,=E,% 1, 
HO =o. (3-5) 


(we assume that bound states do not exist) and unknown functional »|¢] should be 


determined so as to make exp(iS)Hexp(—iS) as diagonal as possible. 


We will fix »[d] by the requirement that the expectation values of the total Hamiltonian 


with respect to the approximate eigenstates, exp(—iS)@,, should be minimum : 
0(exp(—iS)@,, Hexp(—1iS)@,) =0. (3-6) 


|| is then given by 


w[6G) J=3 tan 8) (3-7) 


where 


Peat ap ees M,=M—oM, 


as was shown by Sawada’. When eq. (3-7) is substituted in eq. (3-4), exp(—iS) 
coincides exactly with the Foldy transformation function”, although eq. (3:3) is, in 
our case, not a unitary transformation but an approximate representation of the “ cloud” 
of the physical nucleon. 

Although approximation (3-3) is indeed a very rough one as will be easily seen by 
eq. (3-8) below, it does include the main contribution of virtual nucleon pairs in the 
charge independent pseudoscalar meson theory. One may see this circumstance by observing 
that the matrix element for the creation or annihilation of a nucleon pair of 

9&1 


Sy - 
exp (iS) H exp(—iS) =H)+ ThE | PY 51 Ou (Tibi) Pd x 
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aa Wipe pp pret ie joe dxa))pdx+-- (3-8) 
+ | ys, (2M)? 4 [ 


is very small compared with that of 


HA, iO Z lg | Fr .cdidede. (3-9) 

Although the second and the following terms in the right side of eq. (3-8) are not 
included in our approximation (3-3), the contributions to the vertex operator of these 
terms ate taken into account at least up to the second order in gy by the structure of 
eq. (2-5) or eq. (2-6). Thus we calculate the vertex operator using the approximation 
(3-3) in spite of the well-known insufficiency of the perturbation theoretic treatment of 
the above-mentioned terms, since our aim in this paper is to see the effects of virtual 
nucleon pairs. Moreover one may easily take into account the contributions from the 
P-wave pion cloud in the static approximation by replacing the static part of our results 
by that obtained using the Chew-Low type or the dispersion theoretic approach. It is 
to be noted that we may not say about the validity or the invalidity of the approximation 
(3-3) by inspecting eq. (3-8), since the degree of approximation is very much dependent 
on the structure of the equation where the approximation is being made. 

Now let us make replacement (3-3) in the right side of eq. (2-6). One may 
rewrite them in a manifestly covariant form. The same result, however, can be far 
more easily derived starting from eq. (2-5) as follows. 

We assume the existence of a unitary operator U, which connects the ¥,,’s and @,’s 
defined by eq. (3-5) in the following way 


eT =U. (3-10) 


where bound states are assumed to be not existent. One may rewrite, then, the second 
term of the right side of eq. (2-5) as follows: 


u(x) (Po, Ty) 7A (y)juS) Fo) u(y) 
=u (x) (PD, T(x) N(y) ju lS) %) u(y) , (3-11) 
where 
N(x) =U" n(x) U, ete. . (3-12) 
As will be seen by 
7 (x) =U- 7 (x) U=U- a (x, O)e *#*oUJ 
= et Ton (x, O)e* Meo=sie* Hotty (be; 0) et Horo, 
the time development of (x) is given by the free Hamiltonian. For actual calculation 
It Is convenient to represent U by the operators at the origin of the time. 
If one uses exp(—iS) in place of the above U, where S is given by eqs. (3-4) 


and (3:7), eq. (3+12): becomes chesconvattent expression of the corresponding terms of 
eq. (2-6) with the approximation (3-3), where 
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v/] (x) = 192, Ve iis Te (x) eo (x) + OM ee P(x) , 
0(y) = 19 Z, Zo P (ye 7, 7.8: (y) FOMG(y) eH *™ , 
EO) HIE) +5724) , 


Jn 6) SteZ. 0 C7, a (sa a ) oes 


ir; “eee cine w b(4) , 
He) = sor Lebel y 0.0; > (3-13) 


and 
s(x) =Ts 0; w\d (x) |. 


We have dropped the terms proportional to the normal of the spacelike surface since it 
does not contribute to the S-matrix. 
We have obtained, after all, the following expression which we will evaluate in the 


next section : 
ieu(x) I’, (xy 3 €)u(y) 
=ieZ,ii(x)T,7,u(y)0(x—€) 0(F —y) 
+i(x) (PD, Ty (x) 9 (y)ju.(F)%) u(y), (3:14) 
where (x), 7(y) and j,(€) are given by Eqs. (3-13). 


§ 4. Calculation of (3-14) 


The second term of the right side of eq. (3-14) may be evaluated straightforwardly 


using the prescription given by Kinoshita and Nambu and by Anderson”’, since the time-depen- 


dences of 9(x), 7(y) and j,(¢) are given by the free Hamiltonian. 


For the meson current contribution the result is given by 
(PD, T(x) (y) ju" (5) Po) 
= — #5 Z,arsSo(x— rode —8) Oey te 7) -I™(x—y) , (4-1) 
where 


I™(x—) =| dudva(u)a(v) exp| —H(u!-+ 0°) +2 (w-2) 4, eee | 


‘heyp agua qe t| 2) 


4,(x—y) 
4,(0) 


x [1-2 et +0) 4-2) 


and 
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a(u) = (27) “| ai Alaje*, 


Alal=9Z,Z.7 cos w{V/ Jp(0) A) +0MG" sin »[V/ Jp(0) 2) 


and bold-face letters 2 and wu represent vectors in the isospin space. 


The nucleon current contributions are given by 
(D,, Ty (x)= (y) ja ©) %) 


== ZilrsSr@— A) SeE—N7s{l + ool} 


+S,(x—§)7,Sr(F—y) {E+ 7, 1,"} 
+7375 Sp(x—¢) Tals Sr (¢ =} Ea 
+7,5%(x—S) VulsSr(S—y) 15 17}, (4- 3) 


where I", ---, I,” are very complicated functions of J7'(0) I-(x—<¢), Jz'(0) 4-(¢—y) and 
4;'(0)4,(x—y), the detailed forms of which are given in Appendix A. 

As is shown in detail in Appendix B, only the zeroth order terms with respect to 
45'(0)4,(x—€), 47'(0)45(€—y) and Jz'(0)4,(x—y) contribute to the observable 
quantities. As a result of this fact, it turns out that the principal effect of virtual 
nucleon pairs to the electromagnetic structure of the nucleon in the charge independent 
ps-ps theory reveals itself as the change of the normalization caused by the virtual 
dissociation. This seems to be a characteristic property of the ps-ps theory, owing to 
which evaluation of I”, I”, ---, I," becomes extremely easier. 

At a first glance one might think that J;'(0)4J-(x—y) in the exponential factor 
in the integrand of I”(x—y) corresponds to the propagation of infinite number of pions 
from y to x, but this is not the case, as will be discussed in Appendix C. 

Observing that w|¢] is contained in I” and I” only through the forms, 


cos’ w[4]= {1+ (1+90° 4 My) -™"} /2, 
and 


sin’ »[4]= {1— (1+9° 4M.) “"*} /2 


’ 


we assume that )°/°M,~° may be regarded to be very large compared with 1. This 
approximation seems to be rather good by the following reasons. Observing that 

—¢4 05; dj.(x—y) oa (¥., T'4, (x) 0; ( y) ¥) 

—>(¢,, e’T oy (x) 0;(ye"*%) 

= (%,, T 0; (x) d;(y) P=40,, 4dy(x—y), 
one may infer Z,=1 and y,=y Z,Z,~' is exactly the bare coupling constant in our ap- 
proximation, so %/" would be greater than g®. Further the bare mass M, would be smaller 


than the observed mass because of the contribution from the vacuum polarization. Ac- 
cordingly we may assume 


ae 
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Io? M.*> 1, (4-4) 
since the left side would be larger than one for most important values of A even if we 
replace gy,” by g?~47 X15 and M,~? by M7’. 

One may expect, then, that the correct results would be between the results obtained 


below and the ordinary second order perturbation results, since, as will be seen by shia, 


cos?» |x| 


==sin- vA] 


05+— —-—_—— 2 ee a ee ee = cos yA | 
Doro 9o-> 00 


Fig, 1. 


both cos* [4] and sin’ »|/] are monotonic functions of g,?/?M,~° and their upper and 
lower bounds are given by their values obtained letting y,>0 and g,?/?M,°—00 
respectively. 

Adopting the above approximation (4-3), we may obtain I” and I” numerically 


using the values of the renormalization constants obtained in the next section. 


§5. Renormalization 


It is necessary to evaluate the renormalization constants to obtain the numerical 


results. This is performed in the following way. 
If one has evaluated S’(x—y) =i {)(x), f(y)} 9, Z. and OM may be given by 


7 a 2 4 | Tr lara 5’ (x) Jrood ® 
AG. 
and 


dM=M—42. (t-| g (| dx 
4 u Ox, %=0 


respectively. These relations are easily derived if one expresses S’(x) by the spectral 
functions” : 
5? (2) =[" {85 m) ps (m') + 4s m) palm) dm 
0 


and observes that 


es = i a (m°) dm°* 
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and 
sM=Z,| i (M—m) p,+ Po] dm’. 
0 


: : ' 
It is, however, more appropriate for the contents of our approximation to evaluate S} 


than to evaluate S’ for the determination of Z,. In our approximation S; is given by 
—4S (xy) = (Lo, TE) F(y) Fo) > Pos TP (0) $(9) Mo) 


=-3| dudw exp| — 1} wtv?+2(u-v) a. | 


x | Sey) 6(w) bw) 


—47,Se(—y)rse(u)e(v) | —3(u-0)(1+ =#C—9)_) 


47(0) 
Eu pegeeeae 34 
+43 9 (Ure @*) ) 4,(0) |, sabi: 
where 
b(u) = (an)~* | dae cos »[v/J-(0)4], 
and 


c(u) as (27) —3 | ddev sin w[v_4r(0)4] 
: v 4,(0)4 
As was mentioned before, only the terms of zeroth order in Jz'(0)4J;(x—y) give non-zero 


contributions, so the second term in the bracket of eq. (5-1) disappears and eq. (5-1) 
becomes 


Sp(x—y) =Z,S-(x—y), (5-2) 
where Z, is given by 


(5-3) 


im ho) 


La=5 | du exp(—uw 4)6(u) | 


This fact shows that in the pseudoscalar meson theory the main effect of the virtual 


nucleon pairs to the nucleon propagator is merely to alter the normalization. If, further, 
the approximation (4-4) is used, (5-3) gives 


£o= 1/2. (5-4) 


Z, may be obtained by evaluating (V,,, (L)°—/")0:(€) ¥,), which is essentially 
the S-matrix element for the scattering of a nucleon with initial and final momenta p 
and p’ respectively by an external meson field. Putting 
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0.5) =((P-) bE) =49 2.2, PE), 70) ] 
— 01° $,(F) +Z,16.(€) 474), 
we obtain 
(Fy, O,(0) ¥,) > (Py, &*O,(0) e-@,) 
=19 ZZ." (Vy, Ge* 7,7." bO,) 
=19 ZZ," (i(p')ist:u(p)) {1-2 (sin? w),} , (5-5) 
where u(p) is the usual Dirac spinor wave function. Using 


(FOOL = CS CM ee CS) P) |e=o 


= —| (py —p,) 4-12 | 472 —p)igl.(p, p) Ze" » (5-6) 
and remembering that in our approximation 
—[ (pp)? +H] 4e(p'—p)=1 0 and = Z,= 1, (5-7) 


(5-5) should be equal to ig(a(p’)y;t;u(p)) when p’ approaches to p. Thus Z, is 


given, in our approximation, by 
Se SO 
eit 


. cS iy 2 A“ (0))? 1 EE 
=—1—2(27)12 {1 ( Ge (40) ) = °) 92/29 
eo SR ale meee 


(5-8) 


ea|bo 


Approximation (4-4) is made in the last step of eq. (5-8). 

0M is obtained by calculating the S matrix element for a nucleon propagating freely 
or by calculating directly the expectation value of the total Hamiltonian with respect to 
the physical one nucleon state in our approximation, 


Both methods yield the same result 


2 f2 \1/2 
)M=M,{( 1+ Goo )—1} 
i) ( ar irae é 
2 (9, 4560) 42)" ) (5 o) 


where the approximation (4-4) is used in the last step. 


§ 6. Numerical results 


Now we may evaluate I” and I” numerically, using the values of the renormalization 


constants obtained in the preceding section and making the approximation (4-4), to obtain 
[= g"? F 
t= 49" 4e(x—y) 
TL =t9" 4y(x—9) ; 
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Lp=L'%=1,"=1,"=0, (6-1) 
where 


2 9 ES, 8 a os 8 ¥ te e 
(ee =I de WE (1+) g=K3+—)9 . (6 2) 


All the divergences appearing in I” and I” are cancelled out by 0M, and the vanishing 
of I”, -«-, I” is also the result of cancelling of 0M-term with the other terms. 
Vertex operator, eq. (3-14), may be written in the momentum space in the following 


form : 
Ge p) = £5 teh ut) = yeh) > (6-3) 


where k,=p,—p, and it is assumed that [’,(p’, p) is operating on a(p’) and u(p) on 
both sides. We have written Z, in the first term of the right side of eq. (6-3) as 
Z',, which should cancel the over-all divergence, to distinguish from Z, discussed in the 
preceding section, which should cancel the internal divergence and should be put equal 
to 4 in our approximation. 

Inserting eq. (6-1) into eqs. (4-1) and (4-3), H,(k) and H,(k) are given by 

F(R) = Z, filo” (R) + ret AP (R) +752, fi (R) (6-4) 

where f;°,?(k’) are the ordinary second order perturbation results with the coupling 
constant y’, and superscripts (ms), (mv) and (m) mean the isotopic scalar part of the 
nucleon current contribution, isotopic vector part of that and the meson current contribution, 
respectively. It is seen from eq. (6-4) that the effects of virtual nucleon pairs reduce 
the isotopic scalar and vector parts of the nucleon current contribution to 1/2 and 1/6, 
respectively, relative to the meson current contribution compared with the second order 
perturbation results, since Z,=1/2 and Z,=1 in eq. (6-4). Momentum dependent 
effects do not appear because of the vanishing of the terms including dz‘ (0) dp(x—y), 
etc, in (4-1) and (4-3). To determine Z,’, we divide /’,(p’, p) into the proton 
vertex, /’?, and the neutron vertex, /’;"; 


Pips py =e Pps P) teal CP’, p) » 


where 
PP”, p=) =Z/r7.+7.8P (&) + a Oy k, Hy? (Re), 


ARR) =ZfPE) ALAMO +ALL, 


* Which of the two coupling constants, g°/4z and 9°/4z, should be identified with the Yukawa 


constant, ~—0.08, is not clear in our calculation. One may decide it, for instance, by investigating pion- 


nucleon scattering in a similar way. Numerical results given below for the ratios of the anomalous magnetic 
moments of the proton and the neutron and of the mean square radii of the charge distributions of them are, 
however, not dependent on the value of the coupling constant. 
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nr n (£2 I. n (D2 
Lee, p=7,8f © the Hse Ce); 


FH y'9 (RF) =Zo fi” (R) —Zofid” (RP) —Za fp (F). 

By the requirement [°?(p, p)=7,, Z.! is determined to be 
Zo =1—H,? (0). (6-5) 

Further to satisfy the requirement that the total neutron charge should be zero, i.e. 
I's’ (p, p) =0, we modify I?(p’, p) to 

nr ~- n ( £2 n 1 n ( £2 

ie (p', Pp) =7, 1H, (Rk) cant g (0) } orn; Cavky lds (R°). 

We may then write the renormalized vertices in the following form: 


‘Pp, 
Ae si 6 ag 3 eas a) ar Ouukutier ge)’; (6-6) 


where the anomalous magnetic moments /”” are given by 
P= HP (0) , (6:7) 


the form factor of the vector interaction by 
Fev) =|} +H) —He(0), (6-8) 


and that of the tensor interaction by 
FP") =Hp"(R) /Hp (0), (6-9) 


The numerical value of the ratio of the isotopic scalar part to the vector part of 


the anomalous magnetic momont becomes 


[mela Ga 9. (6-10) 

eRe | 
This result is very large compared with the experimental value, 0.03. However, we 
should compare our results rather with the second order perturbation results than with 
the experimental values, since we have essentially calculated the corrections due to the 
virtual nucleon pairs to the second order perturbation calculation, as will be easily seen 
by the fact that letting S=O in eq. (3-3) one obtains the second order perturbation 
results. The value given by (6-10) is better than the second order perturbation result, 
0.77, and almost the same with the value, 0.51, obtained by cutting off the momentum 
of the virtual meson in the second order perturbation result at the nucleon rest mass.’ 


Mean-square radius of the “charge” distribution, (r°),?", may be obtained by 
Fee) =| i oe ee (6-11) 


the numerical values of which are 
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“Gi =. 4x M? x 


de 


1.37---proton, 
lew | pro a 


— 0.84---neutron. 


The ratio of the above value for the proton to that for the neutron is worse than the 


tatio of the second order perturbation result : 


: 0.487, 
2\, @)p.n— ene : (6- 13) 
cr) _* 
: a 4a M* —0.128, 


and is almost the same with the results obtained by cutting off the momentum integration 
. - - 11) 

in the second order perturbation calculation at the momentum slightly smaller than M. 
The mean-square radius of the “ magnetic moment ” distribution, (r°),””, may also 


be obtained in a similar manner as 


9.5M~?, 
(rp = (6-14) 
3.9M-, 
or 
1.01 
eee X10>* cn*. (6-15) 
1.08 


where in eq. (6:14) we have used the results obtained in this paper for the anomalous 
magnetic moments and in eq. (6-15) the experimental magnetic moments are used. 


$7. Conclusion and discussions 


We have evaluated the effects of virtual nucleon pairs to the electromagnetic structure 
of the nucleon in the charge independent pseudoscalar meson theory, starting from general 
formulae (2:5) and (2-6) for the vertex operator. As a substitute for the complete set 
of eigenfunctions of the total Hamiltonian, ¥,, we have used exp(—iS)@, given by 
(3-3), which is considered to include the main contributions from the virtual nucleon 
pairs in the pseudoscalar meson theory. This is our essential approximation, Vacuum 
expectation values of the T-product have then been evaluated, and the rather complicated 
results obtained in eqs. (4-1) and (4-3) could be extremely simplified by the fact that 
only the contributions from the zeroth order terms with respect to J;'(0)4,(x—y), etc., 
do not vanish. This fact seems to be due to the special character of the pseudoscalar 
meson theory. Therefore the effects of the virtual nucleon pairs to the vertex operators 
reveal itself essentially through the normalization factors which affect the relative ratios 
between the isotopic scalar and vector parts of the nucleon current and the meson current 
contributions, respectively, 

In deriving the numerical results we have further used the approximation (4-4), 
which seems to be rather good as was discussed under eq. (4-4). 


Final results are given by multiplying, by various constant factors, the isotopic scalar 


Oe 


Effects of Virtual Nucleon Pairs to the Electromagnetic Structure of the Nucleon 231 


and vector parts of the nucleon current contributions and the meson current contributions 
calculated by the second order perturbation, respectively. These circumstances will be 
similar in cases of other processes such as meson-meson scattering, 7-decay of the neutral 
pion, etc. 

Numerical results are almost the same as that obtained by cutting off the momentum 
of the virtual meson in the second order perturbation calculation at the nucleonrest mass 
both for the anomalous magnetic moments and the mean-square radii of the charge and 
moment distributions. 

Merits of our method would be that the calculation is Lorentz-covariant throughout 
and that the approximation is made only once in eq. (3:3) except in eq. (4:4) which 
may be justified, so that one may rather easily examine the validity or the insufficiency 
of our approximation. 

Our approximate state vectors, (3-3), are of course very crude, and the contributions 
to the vertex operator from the P-wave pions have been taken into account only perturbation 
theoretically. Although one may easily take into account all the effects of the P-wave pions in 
the static approximation by replacing the static part of our calculation by that obtained 
using either Chew-Low approach or the dispersion theoretic one, it would be necessary 
for the inclusion of the correlation between the P-wave and the pair effects, which might be 
very important, to improve the approximation (3-3) so as to take the P-wave cloud into 
consideration dynamically. Such a problem has not been treated in the present paper. 

The authors are grateful to Professor D. R. Yennie for sending them a preprint of 


reference 1) prior to publication. 


Appendix A Explicit forms of I,", ---, I, in eq. (4-3). 


Functions [,”, ---, 1,” in Eg. (4-3) are given by the following equations : 


I" =| dudvdw exp| —}(u!+0° +10! 


where 

J,8= — 34, (u) A,(v) 8 (w) {4p (x—y) 
—} (u?+0°) 4p (x—y) —guv(1t+4e (x—y))}, 

J? =3 A,(u) A, (v) 9 (w) {4p (xy) —3 (u? +0) dp (xy) 
—juv(1+4;(x—y))} 
— A, (u) A, (v) B, (w) {4p (x—y) —§ (uw? +0") de(x—y) 
—}w'd,(x—£)4,(E—y) —guv(1+4e(x—y)) 
—Avw(dp(x—€) +4p(x—y) 4p(y—§)) 
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—Awu(4e(y—§) +4,(x—€) dr (x—y))} 
—4A,(u) A, (v) B, (w) {4,(x—€) 4,(6—y)(6— (u?+ 0? +200?) 


—Suv4,(x—y) —2vw4,(€—y) —2wud4,(x—é) + ay 


—4,(x—y) (1+ $0? 4s (x—6) +$ 0° dP (y—6) 

~4 (w+ 0°+w")) —un (Fe (x4) + de (y—4) —P) 
—Low(5 4y(x—§) +3 4¢(x—y) 4p(y—)) —h (5 4p(y—F) 
+34,(§—x) 4,(x—y)) 


+ [21 (w+ 0) de@—2) +09) Ie(y—2) 


+uv(4e(y—&) +4-(x—F) 4p(x—y)) +ow(4-(x—y) 
+ 4,(x—€) 4,(€—y)) +wu(1+4de(x—))} 
X {u?dp(x—y) 4p (x—€) + (w+ 00°) J-(y—$) 
+uv4,(x—&) +4,(x—y) Ip(y—F)) +00 (14+ 42 (y—4)) 
+wu(4p(x—y) +4(y—§) 4-(F—x))} 
— {u? 4,2 (x—?) +042 (y—§) +0 
+2uv4,(x—¢)4,(F—y) +20wdJ,(y—S) 
+2wud,(x—&)} X {(u?+e%) Ip(x—y) 
+w?d,(x—&)4d,($—y) +ur(1 +4; (x—y)) 
+ vw (4p(x—€&) +4,(x—y) Ip(y—4)) 
+wu(dp(§—y) +4p(x—F) d-(x—y))} 

Js"=A,(u) 4,(v) 0 (ww), 

J: = A, (u) A, (v) {0 (w) —2 B, (we) 
+B, (w)[1—§ (w dP (x—F) +0 42 (y—8) 
+2uvd,(x—€)4,($—y)) }}, 

J;"=2 A,(u) A, (v) B,(w) [4 -(x—) 
—§ (u® + w") dp (x—8) —h 0? dp(x—y) de(E—y) 
—dwu(1+d(x—6)) —Avw(d,(x—y) 
+4(x—£) dp(€—y)) —huv(d,(y—) 
+ 4p(y—x) dp(x—€))], 


eniaainesiaien a 
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J°= —2 A, (u) A, (v) B, (w)[4-(F —y) 
—$ (v? + 0°) d.(E—y) — fu? dp (x—y) 4p (G —x) 
—fvw(1+4,2(y—£)) —dwu(d,(x—y) 
+45(y—€) 4,(¢—x)) —fuv(4,(x—€) 
+4,(x—y) de(y—€))], 
where 4p(x—y) =4;1(0)4p(x—y) etc. and 4,(u) and B,(u) are given by 


A, (4p(0))16 (x) = (27) -* | du exp(—i6 (x) u) 4,(u), 


uv 


A,|6 (x) |=—9Z,Z,""4 sin w[6]+0M cos of 4], 
A,[6(x) |=9Z,Z,"! cos w[g]+0M6" sin »[ 6], 
B,|0(x) ]=sin’ »[ 4], 

B,| (x) |==(24)~ sin 2»[4], 

B,| 6 (x) |=" sin’ »[ 4]. 


and »{¢] is given by Eq. (3-7). 


Appendix B 


It will be shown in this Appendix B that only the zeroth order terms with respect 
of 4, (0) 4y(x—y), 45° (0) 4e(x—S) and 4;°(0)4,(¢—y) in I%) 1",.-:-, J,” in Eqs. 
(4-2) and (4-3) give finite contributions to the observable quantities. 

We consider the meson current contribution to the vertex operator, but the following 
considerations are also applicable to the nucleon current contribution. 

Fourier transform of the meson current part of the vertex operator, /’{”(p’, p), is 


given by 
(27) ~*0(p'—p—k) IT” (Pp, p) 
ei dxdydé en ipl at tpy + ike Eee ( 3 ) [™ (x—y) ; (B: 9) 
where 
ETE A TOR 


2 Pee % 
=> sR Lis Za stale n(x—y) I5 Ay(x—€) Oy, Ap (E y) (B- 2) 


is the second order perturbation result and I” (x—y) is given by eq. (4-2). 
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Using an identity 


I" (x—y) = [™(y—x) = a ; | d*Leéletem™ [* (a) : 


one obtains from eq. (B-1), 


Po (p!, p)= = oe Ce | aU Oren 2D(pl—L,p—L). —_ (B-3) 


Expanding I” (a) to the powers of Ay @) asta), 


P= p In Geral (B-4) 


(odd powers of 4,*(0)4,(a) vanish by the integration over u and p in eq. (4-2)), 


one may write eq. (B-3) as 


P(e, —) =D [te (O)T* XC, pds (B-5) 
where 
Xi(p’, = coil Chel se (a) F* | he »(p/— —L, ae iL; (B-6) 


We first calculate X,, 1 to obtain 


; ; ae _ a 1 1 | (a—2y)?(1—2) —ildz el . 
X, > 3 &3 dx d . — —_— - wdu > . 
MWe rey cess ate (4). ie Pats ee 


where 


aaa (aig (aa, OPE pace) GBM) 27a PD 
ad J0 (4, —p2)® ; 


[Af eee 2XZ | 2(1—x)z 
(1=2y)*(1—2) = 1 


(1—2y)? (102) r= 1 


22 
=2 << — 4 
c + (ieaae hee > paul (1 x)p+xp’], 


A= cy — C14 — Ody — GAs, 


M? - 2 
u(MPB—ze) up + (1—w) M2, 


~ 4y(1—y)z<(1—2)-y(1—»)z 


u u 


5 “Tie kt ee 
y(l1—y)z A 


C= 


i AW Mace Pi - 
[(i—2y)?(1—z) —1’ 
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B= (1—2x+ 22?) (1—2y)?2(1—z) 
+[1-2+22¢(1—a) ][(1—2y)*1—2) -1] (B-8) 


Intergrations over 4,, 4, and /, may be easily performed to give 
7 A2 AD ERAS arth 
\, di, \, dA, | d4,(4,—p.] * 


= J? In E AF (c, +c) lie= fa Te (c, rath | 


Co Cy 


Q—=(G+G+G) —p2 Gaset =, 
ee G_ of pe | Se (c+) dey 25 
Syetes eo! =p, 
= (G+G)a—),; in a— +e +c) Sees 
COGS G— (Co tg) & spt 
= rey ie. jen 8% Ce) Het 


2 2 
C1 6a Ce Copel ps 


ae eat, Sat pe 


= 
C1 a C3 Co — Pu 


(B-9) 


and 
A2 AZ z A2 " a 
|, da[ > dée[ da ld. pele 
0 0 


xe eh sipaetCa: —atatal a ae aE Mike (ger ul Spi a 


& Sal cpt ap, O16 
Sin swtieies MEP eet seh in 0 (ara) i= “Pe 
g—(G--6) 4 “Pa C6 © —p; 
# beg (cy. 3 G) 4 sg may ae = (eer Co ;) Ae ai 2 
Co Cy o— (cj + 3) =p 
G— (@+c) Yeast (cape In © =e + Cs + ¢,). w= Pu 
zy Cc Co C3 te, -P ¢3) IP aa 
st eye gl ps)? In Ce eh) At == eh 
2 lnloaks —o,/—p? 
ta G A — pao) ia errand 
Co Cx a, alee Oe 


ie ree ow Pu eh pep: >) In co— (c, te) # od 


D cst Ona —py 
2 2 : 
sok ad 9 ae Fr 2 ed 2 (B- 10) 
DiCnGs GC. car pa 


5 fs ; 
Similarly one may write J;(0) as a function of cut-off momentum, A’, as follows : 
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7 AZ e k 
ee ec eee aust fae af 
0a) R42 (27)*Jo ‘) | (P+ pO +4, +45)" 


e | ett m BEL, — wAtE), (B-11) 
ce zy} P+ 
Writing 
i 2 ee 3 i 2 
itp’) Xu (p" p)u(p) =a(p)| 74°C) ni =, ) | u(p), (B- 12) 


we see from eqs. (B-7)~(B-10) that the worst singularities of x'(k), y'(k) and 
[4,(0) P are 
x1 () ~ At In(#27/M?), (RP) ~ # InC#’/M*), 


[4 4,(0) oe LE 


Accordingly, in the limit of A—>00, contributions of X, to 1’ 2” (p’, p) vanishes except 
x'(0) /[4,(0) P. Similar situation holds also for the case of X{(i>1) as will be seen 
from the form of eq. (B-6). 


So, if we write 


wp) Te? P's p)e (=>) Ti" 4-(0) |-*u(p) Xi (p’, p)u(p) 


wit F,° >) | wp), (B-13) 


=i(p) | 7A) +7 
Xji(i21) contribute only to F, (0). Therefore we need only to evaluate the zeroth 
order term in I” to obtain observable quantities. 


Above arguments hold also for the evaluation of Jj", ---, J,". 


Appendix C 


In eq. (4:2) 4y(x—y) appear in the exponential form, and one might think that 
this represents the propagation of many pions from x to y. That this is not the case 
may be seen from the results obtained in Appendix B. 

In this Appendix, we intend to investigate the above-mentioned situation more in 
detail and make some comments on the renormalization procedure proposed by several 
authors") to the case of usually unrenormalizable theories. 

Consider the neutral scalar meson field coupled with the two static sources at %i 
and x», the interaction Hamiltonian being given by 


H, (x) =4¢(x)[0(x—x,) +0 (x—x,) ]. (C-1) 


Then nuclear potential is easily obtained and the result coincides exactly with the lowest 
order perturbation result. 


Next we consider the case, where the interaction Hamiltonian is given by 
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H, (x) =i{ due £(y) [3 (x —x,) +0(x—x,) ], (Ge) 


to compare the situation with eq. (4:2). Nuclear potential in this case is given by 


rr 


V (x—x') = —4hi | \ dudu! f(u)f (u!) (Lem, e~#ONUY € (xx) 


a 
d 


azall dudu’ f(u)f(u')exp [ —4 {wd (0) +u?d™ (0)} ] 


XIm [exp {—}un! de(x—x')} J, (C-3) 


up to the second order in /. Eq. (C-3) is very similar to eq. (4-2), and indeed represents 


the exchange between two sources of may mesons in general. However, if one put 


fu) = (22) ~"| dy exp ize), (C-4) 


then H,(x) is equal to H,(x) and the number of mesons propagating from x, to x, 
should be one. Inserting (C-4) into eq. (C-3) and integrating over u and y, one 
obtains 


V (x—*') = —P A (x—2') ‘ (C-5) 


which coincides exactly with the nuclear potential derived from H,(x), as it should be. One 
may say about the number of mesons propagating from x to y only after performing the 
integrations over parameters. The situation is exactly the same for eq. (4:2), the 
number of propagating mesons being same with the second order perturbation calculation. 


Glauber’? has shown for the coupling 
FG 10 (xe Oe (C-6) 


which is obtained from eq. (C-2) by replacing f(u) =0(u—i/?) symbolically and replacing 
O(x—x,) +0(x—x,) by o(x), that S-matrix contains no infinities after the renormalization 


HM = exp [4 f? J™ (0) ]. (Gaz) 


This type of renormalization has also been applied to the case of usually unrenormalizable 
theories by Arnowitt and Deser, by Cooper and by Green.” If one intends to follow 


their prescriptions, one would renormalize /° in eq. (C-3) by 
A =}? exp [—3 {w'd (0) + ud? (0)} }. (C:8) 


However, if one adopt f(u) given by (C-4), the integrations in the right side of eq. 
(C-3) may be easily performed to give eq. (C:5), where apparently no renormalization 
is needed. So we may not renormalize /” by (C-8) in our case. 

Moreover, if exp| —4{u’d (0) +u?d(0)}] is absorbed into 2”, the absolute con- 
vergency of the integrals in eq. (C-3), which does hold originally, ceases to hold, and 


one encounters some difficulties in the integrations. 
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On a Non-Adiabatic Approximate 
Method in the Non-Linear 
Interaction Theory 


Shoji Ozaki 
Physics Institute, Kyushu University, Fukuoka 


May 27, 1958 


It has been shown by many authors 
that the perturbation method in quantum 
field theory leads to divergent series. B. 
L. Ioffe” showed that e?-powers expansions 
of the total Lagrangian in a uniform electric 
field diverge in the special case of quantum 
electrodynamics, which suggests the diver- 
gency even in the general case of Q. E. D. 
Such a divergence of power series ex- 
pansion is due to the non-linearity of the 
interaction Hamiltonian which results from 
the requirement of the conservation law of 
momentum. For the sake of simplicity, 
we consider the scalar field interacting with 
the self-field. 
divergency due to the power series ex- 
pansion and the divergency due to the 
infinite degrees of freedom of fields, we 
will consider without loss of generality an 
anharmonic oscillater with the interaction 
term /q‘, where 4 is a coupling constant. 


Even in the case of an anharmonic oscil- 


In order to separate the 


later, the power series expansion leads to 
divergence, whose proof will be shown 


later. Therefore it is necessary to seek for 


a non-adiabatic method. Here we develope 
a non-adiabatic procedure by Feynman’s 
functional integral method. The kernel 
for an anharmonic oscillater is given by 


the following : 


K(q’, aye ce he) 
ur 


oe \ oe | exp E | (q?— w?q? — 2/q") dt| Dq : 
' (a) 


The term which contains q' is expressed 
as follows by using the functional Fourier’s 


transformation : 


exp[ —i/ | q‘dt| 


ce {' ~ 


Then 
K(q’, q", t’, t’’) 


(x°—21q*) dé] Dx. (2) 


exp|i/ 


(- 


= | eis | exp |= (g?— wg? — 442g?) at| Dq 


x expLi | x*d Dx, (3) 


which means the weighted superposition of 
harmonic oscillators with a time-dependent 
frequency. The functional integrations with 
respect to'q are exactly performed, which 


will be shown soon later ; 
K(q’, ie hae £) 


{ Ca) gan 


1 1 
Se Ree pa ey, 
| | V 2miy,(t”) aco 
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” 


+y,(t!’) 4’? —2q'q"} | exp| id\ xdt] Dx, 


(4) 
where y,(t) and y(t) are the special 
solutions of the differential equation 


¥+ (w?+44x)y=0 and satisfy the follow- 
ing initial conditions : y,(t’) =0, y(t’) =1, 
y(t!) =1, 9, (7) =0. Hereafter we put 
w=0 for the sake of simplicity. 
to expand the formula (3) as the series 
of preducts of Hermite polynomials, we put 
y,(t) and y(t) as follows : 


t 


Oat cos | f°(=) dz/f(t), 


In order 


yl) =y(8) | de/y (5) 


t 
= sin | f°(2)d=/f() f®, (3) 
where f(t) satisfies the differential equation 


fit+f/f-2fh/Pr=4ix (6) 
and the initial condition f(t’) =0. Now 
we transform the variables x(t) to f(t) 
in order to perform the functional integra- 
tion of the formula (4). Then 


K(q’, gt a 


X exp | ey -cot |fde fea” +f2(t!")q!/) 


er 


—if (t’) f(t") q’q” cosec | fede 


‘ 
t 


meet ot (f"). 98 
2 fie”)! | 


. Lp 
Xexp| 1 | LF + ff 2f (pry de] 


a 


4 exp| : [rae| I 


ee 
412 f (t,)'* i (%) 


(7) 
If we apply the generating function of 


bilinear Hermite polynomials to (7), 


_ EG +y') —2by | 
20—Ff) =) 


= >)-H,@)H.(y), 


n ni: 


then we get 


Reet. te ||, FEE) 


rhea 


cl 5" | 
xere| (n+ 1) fra 


xexp| — 2 (F009? +979 | 


X ACF) 9) Af) 9") 
; Ad ; : . a 
x f+ fif—2f2 ‘} de 
exp| 16d \ {f +f/f f f 
ope | 4 ‘de | ont St 
ad le \f uefa I 
(8) 
It should be noted that as (7) and (8) 
contain the second derivative with respect 
to ¢, the multiplicative law of a kernel 
does not hold in this approximation. As 
the most contribution of the functional 
integral (8) comes from the part f*, we 
take only the part of f* as the zeroth 
approximation. As f which gives the largest 
contribution to the integral is 


f=VA(Qn+1) 4V 2Qn+ 1)?— (8/3)5 
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by using the saddle point method, the 
energy eigenvalues are expanded as 47>” 
for J>1: 


E,=— (2n+1/4)*? 7% 
+27 (2n+ tye dared ost 
La —8/3X Dd? (Inte AS bot - 2. 
(9) 


and the eigenfunction in this zeroth ap- 


proximation is given by 


2Qn+1 \'5, ea 
n(q’) a) ah Ges Mar 


y n! 


xexp| —(28 77 1)"9"| 


— 2n+1 . 1/6 \ 
Jeb fe Sey hans 
x (72 | ; | 7). (20) 


If we tend n->0co, 4-+0 by fixing the 
energy E=(2n+1/4)*°2'” in this zeroth 


approximate eigenfunction, ~,,(q’) tends to 
$:(q’) aw J3l16 exp [E : Baias ane iy/ 2Eq’] F 


whose normalization constant contains an 
essential singular point of 4 and does not 
tend to the usual normalized plane wave. 
This means the convergency of this ap- 
proximate method, whose proof will be 


shown later. 
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An Example of One-Dimensional 
Oscillatory Motion 
in Magnetohydrodynamics 


Nobumichi Mugibayashi 


Department of Physics, Kobe University, Kobe 
May 28, 1958 


It is the aim of this note to analyse an 
example of one-dimensional oscillations of 
a highly conductive medium caused by 
some external force. Such a problem has 
already been treated by Kulikovsky in 
connection with the pulsation of a cylindrical 
plasma in the hign current discharge.” 
Employing a procedure analogous to his 
method, we can find a set of particular 
solutions of the basic equations in one 
dimension. Then we get a peculiar feature 
of magnetohydrodynamics that the moving 
region of the medium can be adjacent to 


This affords us, at 


the same time, an example which does not 


the motionless part. 


conform to the fundamental theorem on 
the simple wave zone in the flow problem. 

The one-dimensional motion of the 
highly conductive medium is described by 
the magnetohydrodynamic equations in 


Lagrange’s form 
ee pels. Op paece 0H? | 
or p Ox 87 


pdx=const, Hdx=const, pf '=const, 


(2) 


where x is the particle coordinate, H the 
magnetic field which is assumed to be 
perpendicular to the direction of motion 
of the particle. The last three equations 
represent the conservation of mass, magnetic 


flux, and entropy, respectively, along the 
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path of the particle. 
We seek a particular solution of the 


equations by assuming 


x=%)+ (a—%) S(t), 
for 0 <x, 4 (region I), (3a) 
for a<x, (region II). (3b) 


X=X), 


For the boundary x,=0 we specify later. 
Here and in what follows the quantities 
with sufhx 0 pertain to the initial instant 
of time. In (3a), €(t) is independent of 
X), ¢(0) =0 and, as will be seen below, 
Sal, = (OO. 

In the region I, the conservation laws 


become 


p= Lo Hath p= Pe, (2') 

7) q] 7) 

with 7(f)=1—¢(). Eg. (1) then 
reduces to 


t A 

Oya dibedh SG, cif) 

ee 7 
in which 
2 a) 
A=— Po 5 

(49 (a—%) Ox, ; ( 

_ B= —- , 0H. , he) 


4z(a—x°) Ax 


and C are constants independent of x,, 
because the left-hand side of (4) is the 
function of ¢ alone. Conversely an ap- 
propriate choice of three constants A, B, 
and C not only characterizes the motion 
of the medium from (4), but also 
determines the initial distribution of pres- 
sure and magnetic field in terms of the 
density by integrating (5) and (6). 
Since 1<7<5/3 for usual media, eq. 
(4) has the solution only for positive 7 (t) 


(or VE.(2) <A)yic as expected. or the 


periodic solution, in which we are mainly 


interested, we have a set of conditions, 
A>0o, B<O and C<0. 


As to other possible choices of constants, 
it suffices to note that they lead ultimately 
to either the compression or the expansion 
of the medium without accompanying any 
discontinuity such as a shock. 

It is clear from (5) and (6) that if 
A> —B, the pressure gradient exceeds the 
magnetic pressure gradient at t=O and 
the periodic motion starts with the com- 
pression of the medium. If A< —B, on 
the other hand, the initial acceleration is 
negative, showing that the expansive motion 
occurs at first. 

In the region II, where the medium is 
at rest, the total pressure must be spatially 
constant. The only requirement at the 
boundary x,=a is the balance of the total 
pressure. No particle crosses the boundary 
and the particle velocity is the same on 
its both sides. These features are similar 
to those of the contact discontinuity” ; 
but here the density, the pressure and the 
magnetic field may change continuously 
across the boundary. Further it is worth- 
while to remark that the region I is not 
of simple wave, while its adjacent region 
II is of constant state. Thus, due to the 
presence of magnetic field, we have an 
exception of the fundamental theorem in 
flow problems that the flow in a region 
adjacent to a region of constant state is a 


— : 
8) We have also an interest- 


simple wave. 
ing case where there is no medium in the 
region II and the total pressure is supplied 
by an external magnetic field alone. On 
these points, we will discuss in more detail 
elsewhere. 


Whatever initial distribution of density 
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fy one may choose, we cannot make the 
total pressure at the boundary x,=0 equal 
to zero, because p— (A/B) 7-%-” (H?/87) 
is constant for a&<x<a, as seen from 
(2’), (5), and (6). That is, if the total 
pressure is zero at the boundary x=aé, 
then it vanishes everywhere in this interval. 
It might be suggested that by putting 
x=%,+ (4+%,)€(t), for —a<x,<0, we 
could balance the total pressure at x=at 
and obtain an oscillation of a limited region 
of the medium (—a<x<a). However, 
this is impossible, since ¢(t) has no neces- 
sary symmetry property. The motion given 
by (3a) seems to need, therefore, at the 
boundary x,=0, some external force like 
a piston or an applied magnetic field. 

In the isentropic case for which po-*=a@ 
is a constant independent of x,, the initial 
density distribution in the region I is not 
arbitrary but is determined by the equation 


a (a—x,). 


Ox, 2ar 


The author is indebted to Dr. T. 
Taniuti and Mr. Y. Kato for helpful 


discussions. 


1) A. G. Kulikovsky, Doklady Akad. Nauk 
S.S.S.R. 114 (1957), 984. See also, I. M. 
Javorskaya, Doklady Akad. Nauk S.S.S.R. 
114 (1957), 988. 

2) R. Courant and K. O. Friedrichs, Supersonic 
Flow and Shock Waves. (Interscience Publishers, 
Inc., New York, 1948). 

3) T. Taniuti, Prog. Theor. Phys. 19 (1957) 69. 


Absorption Effects in Antinucleon 
Phenomena 


Y. Kakudo, T. Kammuri and R. Nakasima* 


Department of Physics, Osaka University 
*Research Institute for Atomic Energy, 
Osaka City University, Osaka 


June 12, 1958 


It was shown experimentally” that the 
ratio of the antiproton production cross 
section for a nucleon in a carbon nucleus 
to that of a free proton is 0.77, and it 
was expected that the reabsorption of 
antiproton produced in a nucleus should be 


On the other hand, 


the cross sections of antineutron production 


taken into account. 


by charge exchange of antiproton for the 
cases of a free proton, carbon and lead 
nuclei were determined as 3, 4 and 3.8 
mb, respectively. The experimental fact 
that the antineutron production cross section 
is independent of mass number of proton 
number of the target nucleus, may possibly 
be interpreted by absorption effects of 
antinucleons in nuclear matter. 

The absorption effects in antinucleon 
phenomena will be examined in the present 
note in order to explain the above men- 
tioned experiments. Since the antinucleon 
production cross sections for elementary 
processes are known in both cases, only 
the consideration for the process of absorp- 
tion in a nucleus is necessary. 

(1) Antiproton production by 6.1 Gey 
proton beam. In the first place, the 
statistical calculation including only the 
effect of the intranuclear motion of nucleons 
indicates that the antiproton production 


cross section per nucleon o’ in carbon is 
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times as large as that of 
This result shows the 


about 8.5 
hydrogen o (1). 
discrepancy of about one-order in comparison 
with the experimental one. Hence one 
must take the absorption and reabsorption 
effects of incident proton beam and produced 
antiprotons. The antiproton production 
cross section (per nucleon) in which the 
absorption effects in the nuclear matter are 
taken into account based on a very simple 


model is expressed as 


Xexp[ —5,/42— (25—5,) /Aa], 


where /,, and 4, are the mean free paths 
of proton for the multiple z-production 
process and of antiproton for the annihila- 
tion process, respectively, and yp is the 
In this 


expression, the effects of transmission and 


nucleon density in the nucleus. 
refraction at the nuclear surface are 
neglected, and it is assumed that the anti- 
nucleon phenomena occur on a line along 
the direction of the incident particle. 

In the evaluation of this expression, the 
nuclear radius R is chosen as 1.3 X A" 
10-""cm and for the mean free paths, 
Ap=2.9X107% cm and 4,=1.1X107" em 
(for a carbon nucleus) are used convenient- 
ly. Inserting these values into the above 
expression, one finds o (for carbon) =0.10 
Xo’ (for carbon) and for the ratio of the 
antiproton production cross section for a 
nucleon in a carbon nucleus to that of a 
free proton, o (for carbon) /o(H) ~ 0.88, 
which has to be compared with the ex- 
perimental ratio, 0.77. The result seems 
to be in good agreement with the ex- 
perimental one since the latter is not so 
accurate. 


(II)  Antineutron production by charge 
exchange of 440 Mey antiproton. In this 
case the effects of absorption are treated 
as in the case of antiproton production. 
For the heavy nucleus, two different nuclear 
models are taken: first as model-A, we 
assume that R, and R,, the radius of 
proton ond neutron distribution, respective- 
ly, to be the same. Secondly, making use 
of the indication that the protons are 
confined to a region of a somewhat smaller 
radius compared with the nucleus as a 
whole due to the Coulomb potential, we 
assume R,>R,, as model-B. 
antineutrons may be produced only by 


Since the 


antiproton-proton collisions at this energy, 
the antineutron production cross section 
(not per nucleon but for nucleus) is written 
by 
Rp 
o=4r0../, | (S—s,) 
0 
X exp| —2(S—s,) /47—2s,/4,"] bdb, 


where o,, is the elementary cross section 
for charge process, /, is the proton density, 
s, is the width of the outer neutron shell 
along the particle path at the distance of 
6 from the center of the nucleus, and /% 
and /* are the annihilation mean free 
paths of antinucleon in the inner region 
with radius R, and those in the outer 
neutron shell region, respectively. We 
estimate this expression for each case below. 

i) For the case of a carbon nucleus. 
It may be chosen as R,=R,=R (nuclear 
radius) and A*=1.1X107-"cm, then we 
find o=0.22 Xo,,. 

ii) For the case of a lead nucleus- 
model-A. We may choose R,=R,=R 
and 4 =1.2X107" cm, and so this results 
in =0.20Xo,,. 
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iii) For the case of a lead nucleus- 
model-B. The choice of parameters in 
this case is as follows, R,=R=1.3X A’! 
2 LOTR cin RAG LB OGF 56107 8.cm: 14 
sel 1g. ems. anderdt— 1.056107" cm. 
By the use of these values, one obtains 
o=0.023 Xo... 

According to the above calculations, 
model-A for lead seems to be consistent as 
to mass number (or proton number) depen- 
dence, but it gives about 1/6 of the ex- 
perimental cross section just as in the case 
On the other hand, 


model-B seems to be excluded since the 


of a carbon nucleus. 


calculated value is smaller than the one 
for carbon by one order of magnitude. 
However, one must consider other effects 
of nuclear matter on the antineutron 
production mechanism such as the diffuse- 
ness of the nuclear surface or other intra- 
nuclear effects. As another important fact, 
it may be noted that the very large 
interaction cross sections of antiprotons on 
complex nuclei are reported experimentally. 

The authors would like to express their 
gratitude to Professor K. Husimi for his 
guidance. They are also indebted to 
Messrs. S. Yamasaki, Y. P. Kwon and K. 


Kobayakawa for many helpful discussions. 


1) Agnew, Chamberlain, Keller, Mermod, Rogers, 
Steiner and Wiegand, Phys. Rev. 108 (1957), 
1545. 

2) Button, Elioff, Segré, Steiner, Weingart, Wiegand 
and Ypsilantis, Phys. Rev. 108 (1957), 1557. 


Equivalent Hamiltonian of an 
Electron Gas at High Density 


Toshiyuki Nishiyama 


Department of Physics, Osaka University, 
Osaka 


June 17, 1958 


Recently Wentzel’ has presented a 
method of the so-called equivalent Hamil- 
tonian of an electron gas based upon 
Sawada’s theory.” In this theory the 
Bloch interaction has been adopted as the 
interaction between electrons and phonons 
together with the mutual Coulomb interac- 
tion of electrons in the explicit form. 
From our point of view, if the Bloch 
interaction is assumed, the mutual interac- 
tion of electrons V® should be of a 
screened form. Alternatively, if the mutual 
interaction V“ is assumed in the explicit 
Coulombic form we should adopt the attrac- 
tive Coulombic potential” as the interaction 
between electrons and ions in place of the 
Bloch interaction. Therefore some modifi- 
cation becomes necessary to derive physical- 
ly meaningful results. _Wentzel’s equivalent 
Hamiltonian is very similar to our Hamil- 
tonian of coupled oscillators proposed in 
previous papers.” Both for an electron 
gas and for an electron-ion plasma it can 
be shown that results derived from the one 
are in fair agreement with those from the 
other. 

We consider that the ionic motion is 
represented by a harmonic oscillation with 
the ionic plasma frequency 2,= 
(472 2N/M)', where N, z, and M are 
the ionic density, charge, and mass. Then 


Wentzel’s dispersion equation is changed to” 
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2A, py? Ss) Ora 


el j 5 
6, (2*) + 2°-—v <P yr—w, q 
p,q (p:q), (1) 
iD — ? 


setting v, = 2," and A,= |/4,| VIG Ane ff 
in Wentzel’s result (7a). 


nuum limit the phonon frequency given 


In the conti- 


by this dispersion equation becomes imagi- 
naty.” 

Next, to gain a great familiarity with 
the current theories and to investigate the 
damping of phonons, we should like to 
diagonalize the Hamiltonian of electrons 
in the absence of the ionic oscillation and 
to see to what extent the frequency of the 
ionic oscillation is shifted by eliminating 
the interaction between the electronic 
plasma oscillation and the ionic oscillation. 
Thus, after separating the plasma mode 
from the Hamiltonian we consider the 
interaction between individual electrons and 
phonons, quanta of the ionic oscillation 
modified by the plasma oscillation. This 
interaction resembles Froelich’s interaction” 
familiar in the theories of superconductivity. 
The coupling parameter corresponding to 
Froelich’s becomes M(q) =7(2,/P,)'”, 
z=1 and 2(2,/P,)'” for the momentum 
transfers of individual electrons parallel 
and perpendicular to the Fermi surface, 
respectively, where P, stands for the Fermi 
momentum and 2,= (3z/5M)" P.g implies 
the phonon frequency. 

Finally, we should like to give the 
change in Gell-Mann’s correlation energy 
caused by the ionic oscillation.” It can be 
obtained from the dispersion equation by 
performing a contour integral enclosing the 
zero points and the poles of ¢,(2°). Tak- 


ing the contour along the imaginary axis 
we get 


ao 


1. p 
de = (3/4) | a°q | (du/27) 
0 “— 


corr, 


x In{1— {X,(u) Y,(u)/A+X,(u))}], 
(2) 

X, (u) = (4e7/7P.q*) R(u/Po 9), 

R(x) =1—xtan'x", (3) 

Y,(u) =27/ (25+). (4) 


This result has been confirmed by a field 
theoretical perturbation methed similar to 
Gell-Mann’s theory for our coupled oscil- 
lator model. A fuller presentation will be 
given elsewhere in the near future. The 
author thanks Professor N. Fukuda for his 


valuable comments and discussions. 


1) G. Wentzel, Phys. Rev. 108 (1957), 1593. 

2) K. Sawada, Phys. Rev. 106 (1957), 372. 

3) D. Pines, Solid State Physics, Vol. 1, (1955), 
Academic Press. 

4) T. Nishiyama, Prog. Theor. Phys. 6 (1951), 
366; 14 (1955), 37, Lecture at the 12th annual 
meeting of the Phys. Soc. Japan, Oct. 1957 
(with Y. Watanabe). 

5) H. Froelich, Phys. Rev. 79 (1950), 845. 

8) M. Gell-Mann, Phys. Rev. 106 (1957), 364. 


Coupling Constant of Kaon-Baryon 
Interaction, I 


Akira Komatsuzawa, Reiji Sugano 
Department of Physics, Kyoto University, Kyoto 
and 
Yukihisa Nogami 


Department of Physics, University of Osaka 
Prefecture, Sakai 


June 19, 1958 


Although there have been several attempts 
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to investigate the dynamical properties of 
kaons and baryons, even their coupling 
type and coupling strength are yet un- 
known. In this note, assuming a simple 
coupling type, we shall point out that a 
physical requirement on the low energy 
K—N scattering restricts the magnitude 
of the coupling constant. 

We assume that the spins of baryons 
are all 1/2 and that of kaon is 0, A and 
+ are of the same parity and (NAR)-, 
(N2'K)- interactions are of the scalar 
type. All effects of pions and anti-kaons 
are taken into account, but direct interac- 
tions of kaons with pions are not consider- 
ed. Despite the large kaon mass, recoil 
effects are found to be fairly small in the 
low energy K—WN reactions through the 
scalar interaction and do not alter the 
results of this note more than 20%. In 
the following we use the static theory” 
from the beginning. 

The T-matrix element for the reaction 
K*+p—K*-+p is written as — 25 (a, 
Wy) h(w,). Here the indexes q and p 
indicate the initial and final kaon momenta 
and w,=(m;+p)'”. 


section is o.=47\h(w,) |? and the total 


The elastic cross 


cross section is 0,,=47Imh(w,)/p. The 
total cross section for the K~ +p reaction 


is denoted by &,,. The Low equation for 


h(w,) is 
} seach 
BA @., == —— \ dw,4 
Pp 9 q 
W,+4m 47%? 
Mis 
ge Frnt (,) a F pop (@) 
greg te Wq+, 
is Fay (Wy) 
: Olay (W 
fof dog eel) 
2 
AT P O, +o, 
mz + 4m 


Here J=9x°+9y", Ja and gy, being the 


renormalized coupling constants of (NAK)- 


and (N.'K)-interactions, and dim (0:5 mx) 
is the mass difference between hyperon 
and nucleon; mass difference between /: 


Flap (Wp) (Op << mx) 


is a formal extension of the K-~ 


and ' is neglected. 
capture 
cross section ¢,,, into the unphysical 
region. 

With the presumption that h(w,) has 
no zeros on the real axis of the complex 


w-plane, we have’ ”) 


Reh(a,) “'=4"'-(w, + dm) {1—4-A(o,)}, 
(2) 


A(a,) = Set onl dw, fe (Wg, Mp) 


Mig 


x Sele) 4g fi (Wy, Wy) Selo) | 


Ce (4) F, (@,) 
mre eo 
2 | dry f_. (Wy, O) Fenelon), 
Oey (4) 
mz+4m 


(3) 


where 
fax (@,, @,) =P-g/ (GO, + 4m)? (ws @,)- 


P stands for the principal value. ©,, is 
the cross section for the pure elastic reac- 
tion cKe ++ p> K> +p. 


scattering under present consideration is 


The low energy 


elastic (0,,=0,,) and the real phase shift 
is given by pcot0(p) =Reh(w,)~ 

Each term of eq. (3) is decreasing for 
W, >mx, and positive in the low energy 
Then, if the value of A is as 


large as 4: A(mx) =1, resonounce occurs at 


region. 

W,=mx. Fora larger value than this, b(«) 
acquires a new pole at w<myx and cortre- 
spondingly an inhomogeneous term should be 
added to the Low eq. (1).” Then the phase 
shift starts from 180°. This implies that 
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a bound state of K and N comes into 
existence. Whereas such a particle with 


strangeness -++1 and mass less than My+ 


my has never been observed. Therefore 
we may put a restriction on A; 
A<A=A(my)". (4) 


Now we estimate the upper bound G. 
Though 0;,/o,>1 at high energies we 


substitute unity for this ratio in the 


max 


integrand of eq. (3). Write | dw, fzlOg Wy) 


mic 
=F ,~(w,) and introduce a function R(,) 


by 
Mmax 
\ diy f_ (7, Gp) Free / Fer 
may 
ek (On) ts. (On): (5) 
Then 


Aaa (mx+ dm) (F.- 
sy hal 6h Ve OY gion (3 
F,.(@,) and consequently A depend on 
C is the contribu- 


As for 
the order of magnitude of C, a crude 


the cutoff energy Max - 


tion from the unphysical region. 


estimate will be given as follows. The 
cross section for K~ capture is expected to 


obey the 1/v-law ; we assume that 


Ohap/ F g == const ‘Pao /' Mx? q (7) 


with p,=/ (w,—4dm)?—m,2. The factor 
px’ comes from that the pion is produced 
in p-wave state. Notice that the denomi- 
nator q of eq. (7) cancels the factor q in 


pe ts (4, Wy) » 


Fcap/Fq <0.3 for the kinetic energy 10~ 


The experiment” shows that 


50 Mev in the center-of-mass system. Eq. 
(7) with the const. < 0.1 roughly fits the 
experiment. Then eq. (7) yields C(m x) 


<F,_ (mx) /30, which is negligible as 
expected from the smallness of the integre- 
tion interval. 

As for R(m,) the main contribution to 
the integral of RF_. comes from the low 
energy region where ©,,, €5'< 0.3ii> The 
cross section for the charge exchang scatter- 
ing K- +p—>K+n is not well known, 
but is likely smaller than o,,.” Therefore 
R(mx) =1.5~2.0 would be an appropriate 
assumption. 

Taking R(mx) =1 and 2, and dropping 
C, A is plotted in Fig. 1 versus @,.,=2 


1.5 
1.4 

£2 
12 
1.1 
1.0 

0.9 
0.8 


0.7 


Fig. 1. Upper bound 2 


~5 myx. If, for instance R(m,x) <2 and 


Wyygx = 2.5 mx (+hyperon mass) is assum- 
ed, we have 4~1.0. A more accurate 
estimate, taking account of o,,,/o,,>1 and 


C>0, gives more severe restriction. 


1) G. F. Chew and F. E. Low, Phys. Rev. 101 
(1956), 1570. 

2) A. Klein, Phys. Rev. 104 (1956), 1136. 
D. Amati and B. Vitale, Nuovo Cimento 6 
(1957), 1273. 
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In cases of other coupling types, corresponding 
values of 4 are not definitely positive, e. g. for 
the pseudoscalar type the corresponding value 
of 4 is negative; there is no possibility for 
bound states to appear. 

4) W. Alles, N. N. Biswas, M. Ceccarelli and J. 
Crussard, Nuovo Cimento 6 (1957), 571. 


Coupling Constant of Kaon-Baryon 
Interaction, II 


Akira Komatsuzawa, Reiji Sugano 
Department of Physics, Kyoto University, Kyoto 
and 
Yukihisa Nogami 


Department of Physics, University of Osaka 
Prefecture, Sakai 


June 19, 1958 


We have estimated in the preceding 
letter an upper bound of the coupling 
constant of the kaon-baryon interactions, 
basing on the requirement that there exists 
no bound state of K+ N. This restriction 
is valid but too lenient. The reason is 
that, if the coupling constant / approaches 
to the upper bound A the cross section for 
K‘-+p scattering becomes close to 4zp™° 
in the low energy region, which is very 
large ‘compared with the experimental cross 
section. In the present note, for a more 
practical orientation about the order of 
magnitude of the coupling constant, we 
make an estimation of the cross section 
for the low energy K‘-+>p scattering and 
compare it with the experimental data. 


The cross section is expressed as 


o (,) =4ER + Pp? + (w, + Im)?” 
x (1 —4-A(w,) 2} (1) 


with 4=g,°+ Qs" and 
, Z 
A (wp) = LEA 1 F_ (w,) 
TC 


+ R(«,)F_, (w,) C(w))}. (2) 


Here we have put o,,=0,,==0, since we 
are concerned with the low energy scatter- 
ing. Eq. (1) includes all the possible 
virtual effects of K, K and pions.” As 
discussed in I, C(w,) will presumably be 
small compared with other terms in eq. 
(a 

For given /, a lower bound for o(w,) 
can be obtained in the following way. Let 
o(w), obtained by dropping C(w) and 
replacing R(w) by an energy independent 
constant r, be denoted by o(w, 1). Since 
R(w) <1 and C(w) >0, a lower bound 


of the cross section is given as 

a(w) >a(a, 1). (3) 
If, for instance R(w) ~2, then 

o(w) ~o(a, 2). (4) 


a (wine for w¥ land .2,.M 4.9 215 neg 
and various values of A are shown Figs. 1 
and 2. Experiments” show almost constant 
cross section of 10~20 mb for the kinetic 
energy up to about 100 Mev in the center- 
of-mass system. Comparing these data 


with Fig. 1 we see 
A017. (5) 
If R(w) ~2 is the case, Fig. 2 shows 
1~0.6. (6) 


The result is not very sensitive to the 
variation of R(w). If larger Way is 
adopted, the restriction becomes more severe. 


On the other hand the 


comparison 
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Fig. 1. o(@, 1) for omax=2.5 mz 

between the experimental cross section for 
the reaction 7+p—K*+A and the result 
of the lowest order perturbation calculation, 
based on the coupling type adopted in this 
note, shows /,°=3 -+ 1”’” which is obvious- 
ly incompatible with eq. (5). If this large 
value is a reliable one, this may be taken 
as an evidence against the scalar type 
interaction. At present, however, this 
should not be taken too seriously since the 
result of the calculation is sensitive to the 
anomalous 


precise choice of 


magnetic 
moments of hyperons and moreover the 
very applicability of the perturbation calcu- 
lation to the kaon reactions is quite doubt- 
ful. 

An investigation of approximate solutions 
of the Low equations for the K and K 
scattering or capture by the nucleon has 
been carried out and will soon be reported. 


1) A. Komatsuzawa, R. Sugano and Y. Nogami, 
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2) 


3) 


4) 
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Fig. 2. o(@, 2) for @max=2.5 mx 

Prog. Theor. Phys. 20 (1958), 246. 

This will be reffered to as I. Notations and 
assumptions on the interaction scheme adopted 
in the present note are the same with those of 
I. 

Previously some non-perturbative approaches are 
done disregarding effects of K and pions. They 
approximately correspond to dropping R(w) and 
C(@) in eq. (1). 

D. Amati and B. Vitale, Nuovo Cimento 5 
(1957), 1533. 

C. Ceolin and L. Taffara, Nuovo Cimento 6 
(1957), 425. 

E, g., J. E. Lannutti et al., Phys. Rev. 109 
(1958), 2121. 

P. L. Donoho and R. L. Walker, Phys. Rev. 
107 (1957), 1198. 

A. Silverman. R. R. Wilsom and W. W. 
Woodward, Phys. Rev. 108 (1958), 501. 
Miyachi estimated the magnitude of the coupling 
constants for various types of interactions apply- 
ing the lowest order perturbation calculation to 
K*—N scattering and K-production by z—N 
collisions. For the coupling type adopted in the 
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A Note on the Magnetic Relaxation 
in a Low Field 


Motoyuki Tanaka 


Department of Physics, Kyoto University, 
Kyoto 


June 19, 1958 


Manetic relaxation in a low field has 
been investigated by many authors” 
phenomenologically, and their descriptions 
can be classified into two types; (1) the 
so-called isotropic relaxation characterized 
by only one relaxation time as a result of 
simple modification of the usual Bloch 
equation, (2) the Codrington et al. relaxa- 
tion” (hereafter abbreviated as the C. O. T. 
relaxation) characterized by two relaxation 
times T,, T,, longitudinal and transversal 
to the instantaneous magnetic field, respec- 
tively. Bloch” and Tomita? have shown 
that the dynamical theory can be applied 
to a.wider range of magnetic field strength, 
and have obtained much more knowledge 
mainly in the case of (1), but they did 
not really attempt to understand the C. O. 
T. relaxation. 

In these circumstances, the writer tried, 
along the line “of Bloch’s formalism,’ to 
see under what conditions the C. O. T. 
relaxation can be specified. 

In. order to compare the relaxation terms 
of eq. (4:55) of reference 3) with the 
damping terms of the macroscopic C. O. T. 


equation” 


1 
pa (MX, 
T, | se 


1 


1 
tp EX (Mx H) | (1) 


where 1/T,==1/T,—1/T,, we expand the 


former in terms of sin@, cos, where @ is 
an angle between the instantaneous field 
and zaxis in the rotating system. » The 
results may be written for the expectation 
values of the spin components (I,*), ¢I,*) 


and (J,*» 
— GE) (CL) — C*) — {9 cos? CL*) 


—7 sin@ cos#<T,*)}, (2a) 
= (o'-EC) Cie yee {7 cosyUE ) 
+ (€—€) sin? 6 (1,*)}, (2b) 


OY er a ee ain 0 CT 
— (§—€) sin cos@ (I,*)}, (2c) 


each referring to the rotating coordinate 
system by using Bloch’s notation. Here 
€, &, 7 and € are also expressed in terms 
of the same @, as in (4-40) in the 


following way, 

€==@0,(w) coshe, ¢/=@,(w’) coshk’, 

7=, (0) —®, (w') cosh’, (3) 

6C=0,(w+0’) cosh(k +x’) 

+, (w—w’) cosh («—k’) 
respectively, where K=hw/2kT, «= 
bw! /2kT, and w!= {(a@—o)?+07}"”. 

It will readily be seen that, when 7~ 
€—¢ and €’¢ are satisfied in (2), the 
generalized expressions (2) are reduced to 
the C.O.T. relaxation (1) as rewritten 


referring to the rotating coordinate system 


using the same field quantity as in (2) 


1 et 
——— (M*—7,, H* 
$ ( ) 


1 


1 : f 
—_-__{H* x (M* x H*)], 
TTA z 
(1) 


where H* is the magnetic field in the 


252. Letters to the Editor 


rotating system. But the conditions under 
which 7~&—€ and ¢/~¢ are satisfied 
cannot be explained without our knowledge 
about the behaviors of @,. Such a 
generalized form %, as in (4-40) being 
of less advantage, we make use of the 
asymptotic behaviors of Y, assuming that 
the correlation functions of the local field 
spectrum may be written in a_ simple 
exponential form. The results may be 
written as P(w,) ~|I'|/w* (for o,<o*), 
or P(w,) ~o*|I"|/w2 (for o* <a,), using 
a characteristic frequency w* of the lattice 
system. We can then obtain a condition 
under which the C.O.T. relaxation may 


exist ; 
ow! (= {(M—)? +073") Sa*>0, (4) 


from the fact that 7 and ¢—€ do not 
vanish in order to be distinguished one 
relaxation time from the other. In this 
case, 77~,(0) {1— (w*/w’)*}, F—C% 
P,(@) {1— (w*/|w’ +|)*}, and then we 
have 7~¢—C, as expected, because (0) 
~D,(w) for oF >. 

By comparing (2) with (1’), we obtoin 


@, (0) —®, (w’) coshx’ = —1/T, 
SUT, =) T;, (5) 


under the condition (4). 
On the other hand, the coefficients ¢’ 
and € of (2) which are independent of 


field quantities become, in this case, 

Ge as26=29,(w') coshe’=1/T,, (6) 
where use has been made of the relation 
PD, (o! +0) ~P, (ow) =9, (w’) /2. Hence, 
from (5) and (6) two relaxation times 


T, and T, introduced phenomenologically 
by Codrington et al. are given by 


1/T,=29,(w') coshe’, 


1/T.= P,(0) +% (w’) coshx’. (7) 


Finally, the writer want to emphasize 
that the condition (4) is difficult to 
attain experimentally, because it means an 
extreme off-resonance or strong saturation 
for the field quantities w, ©, and w,. 
However, a possible way is to keep a 
sample in such a way that * may satisfy 
the condition (4) under a given external 
magnetic field. 

In the usual condition 
(w*>o’', w), it is found that 1/T, ~ 
1/T, ~ 2@,(0) coshe’ because of the rela- 
tions 9,(w) ~9,(o') + 20,(w+0') 
P,(0). 


case of isotropic relaxation, and should be 


experimental 


This is the previously mentioned 
observed as the main feature of the 


magnetic relaxation in a low field. 
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4) K. Tomita, Prog. Theor. Phys. 19 (1958), 541. 

5) R. K. Wangsness, Phys. Rev. 98 (1955), 927. 


A Remark on the Anomalous 
Magnetic Moment in the Static Model 


Yasuo Hara and Ken Kawarabayashi 
Department of Physics, University of Tokyo, 
Tokyo 


June 23, 1958 


In the static model of pion-nucleon 
theory, the interaction Hamiltonian Hyy 
is given by 
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Hp tude! is \ oor ds: cdr. (1) 


When an electromagnetic field is involved, 
an additional Hamiltonian Hyy, is te- 


quired by gauge invariance. 
Hyg —e be ‘ * | p@) oA (46.26) dr. 


(2) 

This interaction is usually neglected for 
the calculation of the anomalous magnetic 
moment since the relativistic perturbation 
theory shows that the contribution from 
the corresponding interaction is negligible. 
However, if the source function p(r) is 
extended, this Hamiltonian, which seems 
to give only a small contribution to the 
magnetic moment anomaly, gives rise to 


an anomalous magnetic moment 
M; .=— 1/ (4 V 67R) 


k dv (k) 
V 2m, dk 


ary: jth Pies 


x Lo. 7. (a1 + bF_1) 
+o_7_(b,,1+a*1) 

=, T_ (bia-Fai,-1) 

=O. $5. dj21-+ 61) | (3) 


where u(k) is a cutoff function. This 
contributes only to the vector part of the 
magnetic moment, so we need to calculate 
the expectation value of this operator only 
for proton. Making use of the technique 
developed by Miyazawa,’ we get the 
following result (in nuclear magneton). 


peel 
(Mint) = a 


By de®) op 
x | o5r® 5p 


OK 


dw lon. \\ Bog @Qvuh) 
27%* O,0,(O, +) [v(D P 
eg 
dk 
a Se af 4 M 
47 BT. [e 
x i a — =| (b)2dk 
a Oye Ox 


etl 
27 O, 0; (OZ + @,) 


(2@,+,) ki | 


oO, @, (@,+@,))? 


Hope at) 
X oa (1) 6 dk dl (4) 


The last integral is obtained by a partial 


integration assuming that the cutoff func- 
We also 


assume the contribution from the (3: 3)- 


tion vanishes at high momenta. 


state to be dominant. Numerical values 


are as follows. 
(| Mine\ = — 4.12 —0.138= 471527 
(K=5/) 
(|Mine|) = —1.34—0.20= —1.54 
(K= 6p). 
The first numbers in these equations re- 
present the first terms in (4) and the 


second numbers are the remaining terms. 


Adding these to Miyazawa’s results, we 


obtain 
y= 1.54—0.18—1.27=0.09 (K=52) 
4y=2.05 —0.18—1.54=0.33 (K= 6/2). 


Compared with the experimental value, 
[y= 1.85, 
It has been believed that the vector part 


of the magnetic moment anomaly belongs 


this result is extremely small. 
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to a low énergy phenomenon to which the 
static model has been successfully applied. 
If so, this result calculated by the static 
model seems very difficult to explain ; there 
seems to exist nothing to compensate this 
large discrepancy. It may be noticed, 
however, that the static model yields a 
fairly good result as was shown by Miya- 
The difficult 


situation we have encountered may be 


zawa, if we neglect M,,,. 


avoided, in our opinion, in the following 
way. 

In the static model, the definition of 
becomes 


charge and current densities 


essentially ambiguous. Indeed, in such a 
theory it is not possible to fulfill the 
continuity equation everywhere inside the 


It is fulfilled only outside the 


source, and inside the source it can be 


source. 


fulfilled only in the form averaged over a 
volume V in which the source is included, 


5 2 
1. e.” 


\ Ee +divS dv=0. 

J Ot 

Expressions for S, and S which fulfill this 
weaker condition and gauge invariance does 
An additional 


appropriate current may contribute to the 


not seem to be unique. 


magnetic moment anomaly which cancels 
the contribution from M,,,, although we 
have not succeeded in finding them. 

On the other hand, if we take the 
standpoint that the static model must be 


regarded as a low energy limit of the 


relativistic theory, we may have to neglect 
M,,, in our case. 

Finally, we remark on the anomalous 
magnetic moment of the nucleus.” In a 
nucleus the propagator is not o(x—y) but 
0(x—y) —p(x—y), according to the Pauli 
principle, so the interaction becomes effec- 
tively an extended one, and we cannot 
neglect M,,,. However, we should not 
use the same p(r) in Hyy and Hyyer 
since their origins are not the same. 

The authors would like to thank Dr. 
Miyazawa for suggesting this problem. 

After the completion of this work, we 
happened to know that R. H. Capps and 
W. G. Holladay had treated a problem 
which is essentially the same with ours.” 
But their conclusion is different from ours ; 
and it seems that their treatment is in- 
consistent when seen from the standpoint 
of the static extended source pion theory, 
since the probability P,° becomes very large 
if we use f°/4=z~0.08 and K=5u~6p 
though they assume P,° to be 10%. If 
we put P,° to be about 60%¢~100% in 
their paper, we find <M.) to be about 
+ 0.78~ +#1.3, which is consistent with 


our result. 
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4) H. Miyazawa, Prog. Theor. Phys. 6 (1951), 
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